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Eine geschlossene Lösung bei Torsion zylindrischer Stäbe 
mit Umlaufkerbe, Ermittlung von Drillsteifigkeiten 
Von Wolfgang Vocke in Dresden 


Im Vorliegenden wird neben der Aufstellung einer geschlossenen Lösung der Umlaufkerbe am zylindrischen 
„Stab und Bestimmung der Formzahlen (Kerbfaktoren) besonders auf das Formänderungsproblem eingegangen. 
Die Ermittlung der Torsionssteifigkeiten wird dadurch erleichtert, daß man längs einer ‚‚Randschale“ inte- 
griert. 5 


Apart from giving a closed solution for the problem of the circumferential groove on a cylindrical rod under 
torsion. as well as determining the form factors (groove factors), the present paper deals in some detail wilh the defor- 
mation problem. Computing of the torsional rigidities is facilitated by integrating along a ‚‚boundary shell“. 


Dans larticle present on traite surtout le probleme de la deformation a cöte d’etablir une solution complete 
de l’entaille rotatoire au bäton eylindrique et de determiner les nombres des facteurs d’entailles. La determination 
des rigidites de torsion est rendue facile par l’integration le long d’une ‚‚cuvelte marginale“. 


B HacToameü CTaTbe IIOMUMO YCTAHOBJIEHKAA PelIeHUA B KOHEUHOM BHUJE NJIA HaceyKku Ha 
UNNIAHNPHYUeCKOM CTOP’KHe HU OUPENEeJIEHUA TEeOMETPNHYeCKHUX TMAapaMeTpoOB, 0C0O60€ BHUMAHHE 
YAENAETCA BONPOCAM U3MEHEHHA POPHMEI. ONPeNejIeHHe FKÖCTKOCTU OTHOCHTENBHO KPyYeHnf 
O6JETyaeTcA TeM, YTO HHTETPAal bep£tcH BAOIB „„KpaeBoüi‘‘ O00A04YKM. 


Die Formzahl 
Zur Lösung der Differentialgleichungen bei Torsion rotationssymmetrischer Körper 


3 

ber en N ME EL RE ER ET (1) 
3 

ne A SE N ES ERERSERNG (2) 


sind verschiedene Methoden bekannt (f: Drillwinkel, y: Spannungsfunktion). Eine weitere ge- 

schlossene Lösung läßt sich mit Hilfe der Superposition finden. Sind y, und y, zwei Spannungs- 

funktionen, dann lautet die neue Funktion: : 
Dame De waa ar (8 


Für y, und y, werden die Lösungen des zylindrischen Vollstabes und des Drehhyperboloids 


x R: 5 
gewählt. Unter Verwendung von Ellipsenkoordinaten mie Gin u-cosv und ve Cof u - sin v 
bekommen wir: 


y=A-(Coju-sinv — B-.cosv(3—cos?vV) .........() 
(a: Wellenradius im Kerbgrund bei der Spitzkerbe.) Die Randbedingung lautet (s. Bild 1): 
Re) en) 


Hieraus: 


. 6). 


Bild 1 


18) 
-1 
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Aus Formänderungsbetrachtungen folgen die Schubspannungsgleichungen 


I 
DE er 
m. 
a a Sn a er 
Die maximale Spannung im Kerbgrund 
T u RL {N 68) 
Max (5 => ‘) url zw Te 
ergibt sich zu 
€ A 3 9 
Te Br A ee a TE SR RE 
en, ? 
a 3 
Mit der Gleichgewichtsbedingung 
2n 6 E A 
A ce c V c 
M= [era 2214-642 3/15 +( 
06 
sowie Einführung des Krümmungsradius o 
3 
' e2\2 
0 a? AI u 1 
a .e@ Se y c2 
@ ar 
Mae Ze ae 
erhält man schließlich, wenn 7, = Fr ist und in den Kerbfaktor darstellt: 
IT * f} 
sn c\ 
ek AS | 
5 . r | 
& = | 
= BEE. 
Ca = 0) c# C i 
en ae = | 
| c* = | 


5 Re 
Bei der tiefen Kerbe je >1) geht Gl. (10) in die bekannte Formel für das Drehhyperboloid 


über, die zuerst von E. Melan!) gefunden wurde, während H. Neuber?) die Abhängigkeit 


3 {i + y: a 


b 
a: X en 


{132 V&4ı) 


Nach Gl. (10) tritt noch der Parameter hinzu. In Bild 2 ist x in Abhängigkeit von ° und 


von zeigte: 
5 5 


aufgetragen. 2 X 
Drillsteifigkeiten 
Integriert man die Gln. (6) und (7), so ergibt sich 
1 RE 3 —— 
I=a|t As Since +3 B-(au 4 area Yan) RE TEEN 


') E. Melan, Ein Beitrag zur Torsion von Rotationskörpern, Tech. Blätter, Prag 1920, Nr. 48 u. 49/50 
°) H. Neuber, Kerbspannungslehre, Springer 1937, 8. 92. D g ‚Nr u. 49/50. 
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Mit o = 0 erhält man eine Spitzkerbe, deren Flanken rund verlaufen (Bild 3), und es gilt 


A 2 


A 
DB . SueW/e te ol La e hal erh . . (12). 


Damit bekommt man für die Länge I! unter der Voraussetzung || und|,|>b:. 


lH l):e] ERTEILEN (13). 


GEST, Gl. (13) sagt aus, daß infolge der Kerbe die 
2,5 „effektive“ Stablänge um die Zusatzlänge 


Sir bi 
lo ae ı) BA, 


b 
Dre 
über in die effektive Länge des Drehhyper- 
boloids: 


3 


. größer als list. Mit > 1 geht Gl. (14) 


b 3% b 
junge de a er (15). 


Bei der tiefen Kerbe kann man also Gl. (14) 
oder (15) für geradlinigen Flankenverlauf 
benützen. 


15 


Bild 3 


Flache Umlaufkerben und -vorsprünge mit geraden Flanken 


Das Problem wird dadurch vereinfacht, daß man längs einer beliebig dünnen Schale inte- 
griert, die von zwei Spannungslinien y — konst. begrenzt wird. 


Bild 4 


Das Teilmoment dM liefert die Hauptschubspannung 
dM 


tee Br ER ERIO)E 
Tau 2.1721: Ah-cos« > 


Mit 7,, = Tun’ 605%, Tyy = Ty„ Sina und dr —=dx:tga bekommt der Drillwinkel die Form 


[s. Gin. (6) u. (7)] 


le une 
I=y.n6) 13: Ah: cos? «x 


7 


er init 


[89] 
1 
* 


2 B ae ö 3 h. 
412 Vocke, Eine geschlossene Lösung bei Torsion zylindrischer Stäbe „a. 37 a oe 


H. Neuber?) und F. A. Willers?) "wiesen nach, daß bei flachen Umdrehungsaußenkerben das 
Orthogonalsystem des Torsionsproblems in die der Randform entsprechende konforme Ab- 
bildung übergeht. Im vorliegenden Fall verwenden wir 


; w2 n 
!= Be BE EIER: SO (18) 
mit & =41 —2n) rn. 
%) @ = 0 (n = 0,5): Für die Spitzkerbe bekommt man 
ze wa —1| su 0 Eee ee (19% 


Man wählt vorteilhafterweise eine Schale nahe dem Rand und unterteilt in die Abschnitte 
L-NundN -—Q. 
Im Bereich L — N lautet mit cos& =1,r=b und Ah= (y + Ay) — y Gl. (17): 


IS Gesal dx 

Din 2. — 

dM b (v En Av)? -- e2 — 12 v2 B= e? — x2 

3 ; 22 A er 

E Form) = 
Da Av<v,v<e,y, und x, sowie e< 1 sind, erhält man: 
f'G 1 
an Ne a Ne en Aare ee ee (20). 


Im zweiten Bereich N - Q sind r und cos« variabel. Es entsteht folgendes Integral: 


2 
ne a 
dM b-Av v ER. 
} FI Rn 


Über die Substitution x = v - Sins und nach Addition mit Gl. (20) bekommt man die Zusatz- 
länge für die flache Spitzkerbe mit geraden Flanken: 


aan Bel 

2 - 2 b 
Zus e \? 1— arctg 

ee > ea je 

ö b b 


e 
+ ul sare tg { 
1— = a 1% = 
4- e\? b b 
Men 
En 
e 
Ey az 
de —e|ir 2 SZ VE A. 
"Gp-o P) ©. # 0: Bei der Kerbe gilt 0,5>n> 0 und für den Vor- 
REG 7) sprung Bild5) 0 >n>—0,. 
_ m _ Im Bereich L—- N können wir mit v < 1 schreiben: 
Az = |z'|- Au und Ah = |z'|- Av 
WW; ZU und im Abschnitt N -- Q: 


Bild 5 
Ax = Ah: cos? 


°) F.A. Wi Aispe ih 
3 Be Willers, Die Torsion eines Rotationskörpers um seine Achse, Z. Math. Phys. Bd. 55 (1907), 


a 1 ES 


Bi 5 th, H » e = : & = 
an nase [Dez.ı097 Vocke, Eine geschlossene Lösung bei Torsion zylindrischer Stäbe 413 


Damit ergibt sich die Zusatzlänge 


eg 
Ax “ Ax 
ED Pe rg 


zu 


Man entnimmt r«), einer Reihenentwicklung und ermittelt die Z,,, in Abhängigkeit von a. 


Ns 


Umlaufende Halbkreiskerbe 


Tiefe Kerbe: Wenn Dee 


in der Umgebung von x = 0 zur Verdrillung bei, so daß die elementare Formel gilt: 


<1 ist, dann trägt im wesentlichen nur die dünne Stelle 


Xp X 
dx dx dx 
us > red — na: = u ı (23) 
Te) [b — y@— a2] eg 
z i Vene 
% 
„b—e P 
und mit < x, erhält man 


b 


Flache Kerbe: Die Abbildung 


e? 
ne 
liefert für den Bereich L- N mite<!l: 
I 
a 35 
le ee Ey a er 22) 
und für den Bereich N — 0: 
Yyı Yı 
Pe dy re ZH U 
20° m |m-Al-sina be.Av Keule 


Yo 


Yo 


3 


_ bzw. nach der Substitution — = = „sowie Addition mit Gl. (25) 
e 
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= 
5 


N N 2. 


Ta IT er 


Die flache Halbkreiskerbe hat bei gleicher Kerbtiefe die doppelte Zusätzlänge wie die flache 
Spitzkerbe [vgl. Gl. (21)]. 


Die Auftragung von = über —- für die 


flachen Kerben in Bild 7 zeigt, da sich die 
Werte auf die glatte Welle beziehen, bei den 
Vorsprüngen negative Zusatzlängen. Bezieht 


Kerben 


Vorsprünge 


0,05 0,70 0,15 


Bild 7 Bild 8 
man dagegen auf die effektiven Längen nach der elementaren Theori i 
Te « a ‚me e, dann ergebe 
für die Vorsprünge positive Werte (Bild 8). seen en Biel zanich 


Zur Darstellung der Zusatzlängen für den ge t i 5 AI TORR : : a 
nolationsfpimel 8 gesamten Bereich der Kerbtiefe wird die Inter- 


e 


Be ze 
ne De e 1 ) | 
Zus a . i De flach as (5) ; Izus, Hei Tr SOSE (27) 
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es . . I, = x 
benützt. In Bild 9 ist ui > a für die Spitzkerbe (0, — 0) und in Bild 10 = 5 m 


für die Halbkreiskerbe über rn aufgetragen. 


Hyperbelkerbe allein (@=0 


2 


Bild 10 


Zusammenfassung 


Für das Torsionsproblem von Stäben mit kreisrundem Querschnitt und veränderlichem 
Längsprofil wird mit Hilfe der Superposition zweier Spannungsfunktionen eine Lösung für den 
Kreiszylinder mit Umlaufkerbe gefunden. 


Des weiteren werden genauere Drillsteifigkeiten ermittelt und in Form von Zusatz- 
längen dargestellt. Im einzelnen werden die Spitzkerbe, der Spitzvorsprung und die Halb- 
kreiskerbe behandelt. 


Eingegangen am 28, Dezember 1956. 
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Über die instationäre Überschallströmung durch 
Schaufelgitter mit Rückwirkung 
Von I. Ryhming, SAAB, Linköping/Schweden 
(Aus dem Institut für Theoretische Gasdynamik der Deutschen Versuchsanstalt für Luftfahrt e. V., Aachen) 


; jonäre Strömung durch ein Überschallschaufelgitter mit Rückwirkung in das Leitrad wurde in 
a [1 Tohanden. Wenn man die Relativgeschwindigkeit zwischen Lauf- und er be- 
rücksichtigt, treten instationäre Erscheinungen auf. Dabei kann man quasistationäre und En u ge 
Erscheinungen unterscheiden. Die Quasistationären haben praktisch die größere Bedeutung A 1, ni 
ganz ähmlich wie die Stationären behandeln. Der Hauptteil der Untersuchungen ist daher den typisch in- 
stationären Erscheinungen gewidmet. 


The stationary flow through a supersonic cascade with reachion into the guide vanes has been treated in a 
previous paper [. 1). n the dee velocity between rotor and guide vanes is taken into account, phe- 
nomena emerge which may be classified as quasi-stationary and non-stationary proper, respectwely. 
stationary phenomena, which are the more important ones. from a practical point of view, may be treated in 
much the same way as the stationary ones. Hence the main portion of the present inwestigation is devoted to 
the non-stationary phenomena proper. 


Le courant stationnaire & travers ume grille & paleties supersonique avec reaction dans la roue fixe a ete 
trait& dans-un article anterieur [1]. Si l’on tient compte de la vitesse relative entre la roue mobile et la roue fixe 
des phenomenes intra-stationnaires se presentent. Om y peut distinguer des phenomenes quasi-stationnaires 
et iypiquement intra-stationnaires. Ceux du caraciere quasi-stationnaire sont pratiquement plus importants 
que les autres ei peuvent tre trait&s de fagon similaire aus phenom£nes stationnaires, par consequent la section 
principale des examens est consacree aux phenomenes typiquement intra-stationnavres. 


VCTaHoBHBIIeecH TeyYeHHE C CBEPX3BYKOBOH CKOPOCTBIO Yepe3 NONACTHyF PeINöTKy, COIPO- 
BOKIAeMOE OÖPATHBIM Bo3NelicTBueM Ha HAaNpaBıIAmmmi anmapaT, paccMaTpuBaloch paHee 
B paÖorte [1]. Ecım mpmHunMmaerca BO BHHMAHHE OTHOCHTEIbBHAA CKOPOCTL MEKAy HUCKOM U 
HANPABAIAMINHM anmapaToM TYyPOHHbI, TO IPMHXONAHTCH YYHUTbIBATb H BO3HHKAMWIIMEe HecTa- 
IMOHapHpIe ABıTeHum. Ilpm 3TOM MO}KHO Pa3ıImyaTb ABJIeHUA KBA3HCTANMOHApHEbIeE U THLNHY- 
HO HecTammoHapHbIe. KBa3ucTaumoHapHbIe ABJIEHMHA HMEIOT IPAKTHYeCKH Ö6ONBbIIee 3HAYEHHE 
HU NONaMmTCcH O00Ppa0oTkKe NPHöMaMHu, BIIOJIHe AHAJIOTMYHBIMH IIPHMeHAeMbIM B CiIyyae YcTa- 
HOBUBIIETOCH TeuyeHum. IIoaTOoMy OCHOBHAA YACTb MHCcJIeNoBAHHÜ MOCBAMAECTCA THUIHYHO 
HecTauUMOHAPHbIM ABJIEHMFIM. 


Vorwort 


Der Zweck dieser Arbeit soll sein, ein Bild davon zu geben, wie instationäre Vorgänge 
sich in der Strömung durch ein Schaufelgitter abspielen. Um die instationären Effekte besser 
erkennen zu können und das Problem nicht unnötig schwierig zu machen, sind dabei mehrere 
Vereinfachungen vorgenommen. Besonders interessant ist die Strömung beim Eintritt in das 
Laufrad bei hohen Umlaufgeschwindigkeiten. Die Schaufeln werden als gerade Linien ange- 
nommen. Die Strömung wird als wirbelfrei und reibungslos angesehen und Grenzschichteffekte 
werden nicht berücksichtigt. Die mathematische Behandlung des Problems wird durch das 
Auftreten mehrerer unabhängiger Veränderlichen und komplizierter Randbedingungen er- 


schwert. Wir werden uns daher vor allem bemühen, die physikalischen Einzelheiten der Strö- 
mung darzustellen. 


I. Einleitung 


Für die folgenden Betrachtungen wird ein in üblicher Weise aufgerolltes Flügelgitter ge- 
wählt, welches zunächst aus zwei voneinander getrennten Flügelreihen besteht, dem Leitrad 
und dem Laufrad. Es wird ein rechtwinkeliges Koordinatensystem eingeführt, welches fest 
mit dem Leitrad verbunden ist. Die x’-Achse dieses Koordinatensystems weist in Haupt- 
strömungsrichtung und das Laufrad hat relativ zum Leitrad eine konstante Laufgeschwindig- 
keitL. Um die Verhältnisse nicht unnötig zu komplizieren, sollen die Schaufeln als gerade 
Linien angenommen werden, welche im Leitrad parallel mit der x’-Achse liegen. Die Lauf- 
schaufeln haben gegen die x’-Achse einen Anstellwinkel e, der Staffelungswinkel sei ß und 
Machwinkel«. Die Teilung wird mit ! und der Abstand von einer Laufschaufelvorderkante 
zu dem nächsten oben liegenden Leitschaufelende mit b=b(t), einer linearer Funktion der 
Zeit, bezeichnet (Bild 1). Um die Übersichtlichkeit der Arbeit nicht zu stören, wird dieselbe 
Teilung im Laufrad wie im Leitrad vorausgesetzt, was bezüglich der allgemeinen Gültigkeit 
keine Beschränkung der Darstellung bedeutet. Weiter soll die Behandlung auf solche Fälle 
beschränkt werden, bei denen eine „Rückwirkung‘ im Gitter vorhanden ist. Wie in [1] sind 
eine „Grenz-Machlinie‘“ und ein „Rückwirkungsgebiet‘“ am Ende der Leitschaufeln definiert, 


siehe Bild 1. Es ist angenommen, daß die Schaufeln genügend dicht liegen, so daß keine Mach- 
linie vom Schaufelende vor das Gitter gelangen kann. 


9 
2 
3 
f 
Bi 
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‚ Zunächst werden wir die Randbedingung auf den mit der konstanten Laufgeschwindig- 
keit‘L bewegten Laufschaufeln im x’ y’-System aufstellen. Diese Randbedingung_ besagt, daß 
die zur Plattenormale Komponente 


der Geschwindigkeit W gleich ist a „hukeltraa) 
mit der Schaufelgeschwindigkeit in > 
derselben Richtung. Sind die Kom- Mi ee 
ponenten der Störgeschwindigkei- Ux \ 
ten u und v in den beiden Achsen- NE 
richtungen und die ungestörte Ge- r KDD 
schwindigkeit im Leitrad U,, so ER A Re ee 
lautet die Randbedingung auf der Kagel N ur 
Schaufel: RN} KL Laufrad 
MN 
v+(Uo+ u) -tge RC Abe rg 
—=L-(sinß + cosß-tge) (1). Aa Er: 
. NV Nu 
. Füru=v=0iste=e,und x ” 
dann gilt folgende Beziehung: Bild 1. Koordinatensystem und Bezeichnungen 
L- sin ß 
t ee ee a 5 er iS ; 
5 © Us —L:cosß @) 


In diesem Fall gibt es eine ungestörte Parallelströmung mit der Geschwindigkeit U, durch 
das Gitter. Hier begegnen wir der klassischen Turbinentheorie, welche für das Geschwindigkeits- 
dreieck am Einlauf des Schaufelkranzes eine solche Stellung der Schaufeln fordert, daß die 
Laufschaufelvorderkanten mit der Relativgeschwindigkeit w parallel angeströmt werden (Bild 2). 
Bei unserer Behandlung werden wir immer 
bei vorgegebener LaufgeschwindigkeitL e 
so wählen, daß wir eine wenig gestörte 
Parallelströmungerzeugen, d.h.,wirsetzen 

BI LE, OA), 
wobei ö immer klein, positiv oder negativ, 
angenommen wird. Positives öd entspricht 
einer zu niedrigen Geschwindigkeit des 
Laufrades, negative ö seiner zu hohen 
Umlaufgeschwindigkeit. Wesentlich ist 
nun für unsere Behandlung, daß die 
Gl. (2) nicht genau erfüllt ist. 


Bild 2. Geschwindigkeitsdreieck beim Einlauf eines Schaufelkranzes 


2. Charakteristiken-Theorie für die quasistationäre Gitterströmung 


In der früheren Arbeit [1] wurde das Strömungsproblem im Gitter bei verschiedenen 
feststehenden Stellungen vom Laufrad zum Leitrad mit Hilfe der Charakteristiken-Theorie linear 
und exakt diskutiert. In diesem Abschnitt werden wir die früheren Ergebnisse auf ein Gitter 
mit einem bewegten Laufrad übertragen. 

Zuerst wollen wir jedoch zeigen, daß es im Gitter zu instationären Ausbreitungsvorgängen 
von Wellen kommt, wenn man die Laufgeschwindigkeit des Laufrades berücksichtigt. In 
Bild 13, [1] wurde die exakte Lösung des stationären Problems für b/l = 0,218 dargestellt. 
Zunächst nehmen wir an, daß das Laufrad sich seit einiger Zeit mit der konstanten Geschwin- 
digkeit L bewegt. Weiter nehmen wir an, daß L im Verhältnis zur Schallgeschwindigkeit c«, 
der ungestörten Parallelströmung, in der Düse klein ist. Damit können wir sagen, daß früher 
in der Strömung entstandene Wellen sich schon so weit fortgepflanzt haben, daß sie nicht in 
das Strömungsbild hineinkommen. 

Wir studieren nun die Strömung für Werte z < 0,218. Da das Verhältnis b/! im Laufe 
der Zeit kleiner und kleiner wird, ändert sich die Konfiguration des Strömungsbildes. Dabei 
bewegt sich der schiefe Stoß in Bild 13, [1] hauptsächlich gegen die Strömung, weil der Re- 
flexionspunkt des Stoßes an der Oberseite der Laufschaufel seine Lage von Zeitpunkt zu Zeit- 
punkt ändert. Nach einer gewissen Zeit ist b/l so klein geworden, daß der Stoß unmittelbar 
hinter der Vorderkante an der Unterseite der Laufschaufel auftrifft (siehe Bild 3a). Im nächsten 
Augenblick tritt der Stoß vor die Vorderkante der Laufschaufel und damit breitet sich eine 
instationäre Welle aus. In den folgenden Kapiteln wollen wir dieses Problem behandeln, welches 
wir ein „Ausgleichsproblem‘ nennen können. Bei der Behandlung dieses Problemes ist es nicht 
notwendig, L als klein anzunehmen. 


Yırn)ı Pütz m r En a 7 
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inblick auf das oben Gesagte kann man den Strömungsvorgang im Gitter in zwei ver- 
ea a nämlich: 1 die Laufschaufeln bewegen sich innerhalb solcher Gren- 
zen b/l, daß kein Vorgang auftritt, welcher eine Welle auslöst, und 2. ee 
sieren einen solchen Wert b/l, bei dem gerade eine Welle ausgelöst wird. Es ist offenbar, dal 
die stationäre Charakteristiken-Methode auf den letzten Fall nicht anwendbar ist, weil dort 
eine instationäre Welle sich in der Strömung ausbreitet. Für den ersten Fall jedoch ist die 
‚Grundströmung, d. h. die von instationären Wellen nicht beeinflußte Strömung, quasistationär 
und hier können wir, wie es sich zeigen wird, die Charakteristiken-Methode unter gewissen 
Bedingungen verwenden. ne BR 
Das Ziel ist dann, die Störungen zu berechnen, welche von den wesentlich instationären 
Wellen stammen, um durch eine Superposition der Grundströmung und der Wellen die Gesamt- 
strömung diskutieren zu können. Eine solche Superposition ist nur bei Beschränkung auf eine 
lineare Theorie der Störungen gestattet. Da wir in [1] gesehen haben, daß die lineare stationäre | 
Theorie nicht überall eine gültige Lösung gibt, wählen wir als eine Grundströmung die exakte | 
quasistationäre Lösung. \ A 
Ehe wir uns mit der Berechnung einer einzelnen Welle befassen, ist es jedoch notwendig 
zu untersuchen, wie die Grundströmung von der Laufgeschwindigkeit L beeinflußt wird. Dieses 
' Problem kann man auffassen als ein quasistationäres Problem, und es entsteht die Aufgabe, 
den Unterschied zwischen der stationären Lösung mit feststehenden Schaufeln und der quasi- 
stationären Lösung zu untersuchen, welche dadurch entsteht, daß sich die Schaufeln mit der 


Geschwindigkeit L bewegen. 


Die Strömungskonfiguration, welche wir im x, y'-System mit bewegten Laufschaufeln 
beobachten, wenn wir zunächst von instationären Wellen absehen, ist folgende. Infolge der 
Bewegung der Schaufeln ändert sich die Lage der Reflexionspunkte der Stöße und derMach- 
linien auf den Schaufeln. Es entsteht dadurch ein System von laufenden schiefen Mach linien 
und Stößen. Betrachten wir zunächst irgend ein kleines Feld konstanten Zustandes begrenzt 
von vier Machlinien, so ändert sich wegen der Bewegung die Lage dieses Feldes im x, y’-System. 
Die ursprüngliche Form des Feldes kann auch während der Bewegung deformiert werden, weil 
die vier verschiedenen Machlinien verschiedene Laufgeschwindigkeiten haben können. Um 
den Zustand in diesem Feld berechnen zu können, ist es notwendig, eine solche Variablen- 
transformation zu machen, daß das Abhängigkeitsgebiet des Feldes in dem neuen System ruht. 
Eine solche Transformation ist immer möglich, denn wählen wir z. B. als neuen Ursprung den 
Kreuzungspunkt zwischen jenen Machlinien, welche das Abhängigkeitsgebiet des Feldes tragen, 
dann berechnen sich die Zustände im Feld aus den stationären Gleichungen. Die erwähnte 
Transformation setzt sich im allgemeinen aus je einer Galileitransformation in x’- bzw. y’-Rich- 
tung zusammen. Führt man diese Transformation in die stationäre Verträglichkeitsbedingung 
ein, so. findet man, daß der Typus der Verträglichkeitsbedingung erhalten bleibt. Wegen der 
Transformation erhält man jedoch eine Änderung der Machzahl und der kritischen Schall- 
geschwindigkeit c*. Vergleiche etwa die verwandte Betrachtungsweise bei der Ableitung des 
stationären schiefen Stoßes (siehe etwa [2] Seite 272). 


Es ist also möglich, die quasistationäre Strömung feldweise durch eine geeignete Trans- 
formation der Variablen stationär zu machen. Dies ist als eine Definition der quasistationären 
Strömung aufzufassen. 


Im x, y'-System der quasistationären Strömung gibt es aber auch Machlinien, welche sich 
nicht bewegen. Diese Machlinien sind diejenigen, welche von der Endkante der Leitschaufel 
in das Rückwirkungsgebiet hineinlaufen. Dadurch ist die Anfangsbedingung für die quasi- 
stationäre Strömung festgelegt, weil der Zustand auf diesen Machlinien bekannt ist. Mit neuen 
Verträglichkeitsbedingungen im neuen System kann die quasistationäre Lösung feldweise ge- 
funden werden. Dabei braucht man noch die Randbedingung auf der Schaufel, welche im 


x, y'-System die Form Gl. (1) hat. (Am Ende des Abschnittes 7 ist ein Beispiel für eine solche 
Rechnung gegeben.) 


Mit zunehmender Laufgeschwindigkeit L wächst allerdings der Winkel. Bei einer Lauf- 
geschwindigkeit L c» ist tgem1, d.h., bei dieser Laufgeschwindigkeit und M& = y2 
liegen die rechtsläufigen Machlinien nahezu parallel mit der Schaufel. Dabei ändert sich der 
Charakter der Strömung, indem eine Schaufel mit einer höheren Laufgeschwindigkeit ein anderes 
Einflußgebiet besitzt. Da aber die Strömung bei kleineren Geschwindigkeiten von großem 
Interesse ist, wollen wir uns hier auf niedrigere Geschwindigkeiten L< c» beschränken. Bei 
solchen Geschwindigkeiten ist zu ersehen, daß die quasistationäre Lösung sehr wenig von der 


stationären Lösung in [1] abweicht. Es ist also mit Rücksicht auf die frühere Arbeit unnötig, 
neue quasistationäre Lösungen zu berechnen. 


Pe 
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3. Randbedingungen für das Ausgleichsproblem 


Im vorhergehenden Abschnitt wurde gezeigt, daß im Gitter, wenn eine zunächst be- 
trachtete Laufschaufel eine gewisse Stellung b(f) erreicht, ein schiefer Stoß plötzlich vor die 
Vorderkante der Schaufel tritt, (Bild 12—14, [1]). Damit entsteht ein reines instationäres Pro- 
blem, weil der Stoß sich mit einer gewissen Geschwindigkeit ausbreiten muß. 


‚ Um das Problem behandeln zu können, müssen wir zuerst die Aufgabe formulieren. Im 
Zeitpunkt {= f, haben wir in der Nähe der Schaufel die quasistationäre Konfiguration des 
Strömungsfeldes, welche in Bild 3a dargestellt ist. Zu dem Zeitpunkt {= 1, + At ist der Stoß 
vor die Schaufel gekommen und der höhere, in Feld 3 (Bild 3a) herrschende Druck muß sich 
ausgleichen mit dem auf der Oberseite der Schaufel herrschenden niedrigeren Druck. Es ent- 
steht folglich eine v-Komponente auf der in Bild 3b eingezeichneten gestrichelten Linie vor 
der Vorderkante der Schaufel und auf den beiden Seiten dieser Linie herrscht derselbe Druck. 
Infolge des Druckausgleiches auf der Oberseite und der Unterseite der Schaufel pflanzt sich 
eine Welle in der Strömung fort. 


Pı<Pa<P; 
1. Schaufe/ 


Bild 3. Anfangs- und Randbedingungen für das Ausgleichsproblem 


Da wir uns, wie bereits früher gesagt, auf eine lineare Theorie dieser Störungen beschränken, 
sind zunächst die Strömungszustände in den verschiedenen Feldern in Bild 3a als konstant 
anzusehen. Dies ist eine Annahme, welche sehr gut in der Nähe der Schaufel erfüllt wird. Damit 
können die Randwerte des Problems formuliert werden, indem die Zustände 14% Do Co und 
Ugg; Dos Cap auf den beiden stromaufwärts liegenden Berandungslinien des Ausbreitungsgebietes 
konstant und bekannt sind, siehe Bild 3b. 


4. Aufstellung der Potentialgleichung des Ausgleichsproblems und Herleitung 
einer Lösung des Störpotentials 


Mit der Einführung eines Störpotentials ® lautet die linearisierte Potentialgleichung in 
einem System, in welchem dieMachlinien der Anströmungruhen («', y’, t-System), folgender- 
maßen: (siehe etwa [2], Seite 408, Gl. (81)). 


1 1 
(1 — Mi) Bae + Br a du—2Me — Dn=0......0). 


oo 


Dabei ist ® das Störpotential. Eine Lösung der Gl. (5) ist gegeben durch: 


Br te Dit). 2 de) de (6) 
- \} VI—ME) SEEN F2: Mare kN) NE — Ey? | 
was unmittelbar durch Einsetzen einer Galileitransformation in die Lösung der gewöhnlichen 
Wellengleichung folgt (siehe [2] Seite 409, 
Gl. (85)). P(xyst) 
Für die Differentialgleichung (5) 
und ihre Lösung (6) liegt die Auffassung 
zugrunde, daß es eine Singularitätenbe- 
legung in der ganzen x’, t-Ebene gibt, 
deren Quellstärke auf der Stelle r, &’ 
durch v,(r, &') gegeben ist. Um das 
Potential in einem beliebigen Aufpunkt 
P(«',y',t) zu bekommen, müssen wir 
über das Abhängigkeitsgebiet des Punk- 
tes P in r’, t-Ebene integrieren. Der 
Abhängigkeitsbereich des Punktes P ist 
wieder durch den Schnitt des charak- 
teristischen Kegels für den Punkt P mit 
der x’, -Ebene gegeben, siehe Bild 4. Bild 4. Charakteristische Kegel im «', y', t-Raum 


B En , 2. . Math. Mech. 


Für eine Behandlung des Ausgleichsproblems ist es zweckmäßig, ein anderes Koordinaten- 
system als x’, y'-System zu wählen. Dieses neue Koordinatensystem ist in Bild 3b angedeutet. 
Die x-Achse fällt mit der Laufschaufel zusammen und der Ursprungspunkt liegt mit dem be- 
wegten Stoß fest. Wir suchen zunächst eine Transformation zwischen dem alten x’, y'-System 
und dem neuen x, y-System, mit deren Hilfe die Gl. (5) und (6) in dem neuen System aufge- 
schrieben werden können. Dabei rühen die Machlinien im neuen wie im alten System, so daß 
| wir die Möglichkeit haben, Gl. (5) und (6) zu 
transformieren. 

In Bild 5 ist die Geometrie der Bewegung 
für die Schaufel und den Stoß im x’, y’-System 
dargestellt. Wird die Laufgeschwindigkeit L 
normal und tangential zur Schaufel mit ZL, 
bzw. L,, die Normalgeschwindigkeit des Sto- 
Bes mit a, und die Projektion von a, in der 
‘ Tangentialrichtung der Laufschaufel mit a, 
bezeichnet, so erhält man folgende Linear- 
transformation zwischen dem x’, y'-System und 
dem x, y-System: 


z=2-cose—y'-sine+a,-t, (7a) 
fan 
— ! ” 1 I ” u . 
Bild 5. Geometrie der Bewegung: Schaufel und Stoß Yy Eu sine-+ A Aush L, { i 
oder ihre Umkehrung 
X = zcose+tysine—t-(a,-cose—L,-sine), (7b) 
y=—xsine+y-cose+t-(a,-sine+L,:'cose) 


Diese Transformation ist durch eine Drehung um den Winkel e und eine Galilei-Transformation 
mit L, in y- und mit a, in x-Richtung gebildet. 

Führen wir die Transformation (7) in Gl. (5) ein, so erhalten wir folgende Differential- 
gleichung für das Störpotential ® im z, y, t-System: 


1 M M, 
9,1 -M) +9, —-M)—— 9,20, 220,7 —20,,°M,M2=0 @), 
mit: 
M.-005e + =M;; Dos sa: 
| N . (9). 
M«:sine—— = M,; U'sne—L,=U, 


Dal,„=L sin (? + e) ist (Bild 5), finden wir, falls wir Gl. (2) in dem Ausdruck für M, Gl. (9) 
einsetzen, daß M,„= 0 wird. Durch die Beschränkung auf kleine Störungen Gl. (3) folgt jedoch, 
daß M, im Verhältnis zu M, klein, aber immer von Null verschieden ist. Ferner ist M, immer 
größer als eins. 
en = x, y, t-System bewegt sich die Laufschaufel der x-Achse entlang mit einer Geschwin- 
digkeit: 

+ DL, = 4, sin (&.+ &), Hr L cos (Bee) Re 
Falls wir den schiefen Stoß mit einer geraden Linie vom Machwinkela approximieren, hängt 
a, von L wie folgt ab: 
sin (# + &) 
Mo nee Spa Me a 


Eine Lösung von Gl. (8) bekommen wir durch Einführ der T i i 
(7) in GI. (6): ) ührung der Transformationsgleichungen 


a„=Lsin2o- 


ur ZN |&, 7): dE- dr 
7) Ya—M—MI) es dm H2M, co (EA) RE) +2 My: 0a dm) ya 
Br; 1 
Shi 169 (2): 48 dr SEE A AUTOR 


mit 


e= U —M— Mn). %(—D?+2M,- col— (2 —-9+2M, cl Dy-@- De, 
Im Integrai Gl. (12) steht im Zähler eine unbekannte Funktion /(&, 7), und wir werden zunächst 
die Bedeutung von f(£, r) untersuchen. Da wir ein neues Koordinatensystem genommen haben 
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welches im Verhältnis zum alten x’, y'-System eine Geschwindigkeit in einer gewissen Richtung 
hat, können wir nicht erwarten, daß v, auf dieselbe Weise in Gl. (12) wie in Gl. (6) eingeht. 


. Um die Bedeutung von f(£,r) bestimmen zu können, ist auf eine richtige Abgrenzung 
bei der Integration zu achten. Das Abhängigkeitsgebiet eines Punktes P(x, y, t) ist wieder 
durch einen Kegelschnitt gegeben und die Randkurve des Abhängigkeitsgebietes in der &, r- 


Pfxyt) 


Bild 6. Abhängigkeitsbereich eines Punktes P in der £&, z-Ebene 


Ebene ist durch das Nullsetzen des Radikanden in Gl. (12) gegeben (Bild 6). Bei der Aus- 


führung der Integration integrieren wir zuerst über 7 bei £ = konst., wobei die Grenzen für das 
Integral sind: 


1 RS the 
et = Ma y+M, @—H + Y@—H%- (1 —M}) +2-@—8)-y-M,-M,+y®-(1—M3)] 


mit (13) 
GM Me ee EL 
Für die Integration über & ist der Punkt E in Bild 6 durch: 
1 n 
&&=2-172769C- MM, +0 Kae SER (15) 


gekennzeichnet. Für y > 0 ist dabei das obere Vorzeichen vor o in &, zu nehmen. 
Für ® haben wir dann: 


&E Ts 


9-2 | IMea Sea 2.22.2202...) 


on 


Im allgemeinen bekommt man Auskunft über die Bedeutung einer unbekannten Funktion /(£, 7) 
in einem Integral über eine Quellbelegung ähnlich Gl. (6), wenn man die Ableitung ®, bildet 
und dann mit y gegen Null geht. Dabei steht im Zähler des Integranden ein y, was zeigt, daß 
der Integrand, wenn y gegen Null geht, über all Null wird, außer an der Stellex=&,1=[r. 
Diese Stelle kann daher einen Beitrag zum Integral liefern. 

Berechnet man nun die Ableitungen ®,, ®, und 9, von Gl. (16), so stehen im Zähler des 
Integranden lineare Bindungen zwischen x, y und f. Zunächst bilden wir eine solche Linear- 
kombination zwischen ®,, ®, und ®, daß im Zähler des Integranden gerade y als Faktor auf- 
tritt. Danach wollen wir untersuchen, wie diese Linearkombination sich verhält, wenn y gegen 
Null geht. 


Eine Linearkombination, welche die oben erwähnte Bedingung erfüllt, ist gegeben durch: 


EEE 
© i 13 Hi 
m :-®8,+M,;: 9) —M,)®,=3 _ | | y-(&.r)- (- Sde dEars LIIR 
Um in diesem Ausdruck den Grenzübergang y gegen Null machen zu können, ist es notwendig, 
erst Gl. (16) zu integrieren und danach die oben gegebene Linearkombination durch Differen- 
tiation von ® zu bilden. Wir entwickeln daher /(&, r) an der Stelle &=x und =tin eine 


Taylorreihe: 
END-ICH+LEIE-NHIEDE—Hr...... (8 


und setzen diese Entwicklung in Gl. (16) ein. Die Integration von Gl. (16) kann nun durch- 
geführt werden und es ist danach einfach, die Linearkombinationen Gl. (17) zu bilden. Gehen 
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wir dann in diesem Ausdruck mit y gegen Null, so erhalten wir: 


lim Dans B,-+ m, —M)9,\ = FIed. run Me (19) 
y>0 Co 
für 40. BR i 
Das erste folgende Glied enthält nur die Ableitung /,,, auf welche es natürlich ankommt, 
wenn die Linearkombination Gl. (17) in einiger Entfernung von der Schaufel bestimmt werden 
soll. Hier interessiert aber der Wert unmittelbar auf y = 0, wofür nur der Wert von Gl. (17) 
der unmittelbaren Umgebung verantwortlich sein kann. i 
Die physikalische Bedeutung der beiden ersten Glieder in Gl. (19) 


(Lose) 
ist in dem Potentialansatz von Gl. (8) zu suchen. Es gilt (nach [2], Seite 407) 
en 1 a 90 
— DS ee c) RAINER ar er Ir 


wo der Index 1 irgend einen bekannten Zustand bedeutet. Mit Beschränkung auf kleine Stö- 
rungen gilt dann: 


2 
— Bd, = (W— W)Wı en a ART BR 
oder 3 
a A Fe ir ET De RR (21b), 
d.h. 
z —- MD en (21ec) 
eh u 6% 908 In elle Gl a ; 


Die ersten beiden Glieder in Gl. (19) stellen also die mit M,„ multiplizierten Druckstörung an der 
Schaufel dar. 

Damit ist eine Lösung der Potentialgleichung (8) des Ausgleichsproblems durch Gl. (12) 
gegeben, wobei f(x, £) die durch Gl. (19) und Gl. (21) gegebene Bedeutung hat. 


5. Einige Eigenschaften der Potentialgleichung des Ausgleichs-Problems 


Die Randkurve des Abhängigkeitsgebietes eines Punktes P in der &, r-Ebene wurde im 
vorhergehenden Abschnitt durch Nullsetzen des Radikanden in Gl. (12) berechnet, Bild 6. 
Mit der Auflösung des Radikanden nach £ statt 7 erhält man: 


@— Hd = ct t-—TM)M + Yy& t—-WA—M)+2c t=d)yM, a. 
Ein Schnitt mit 7 = 0 gibt zwei verschiedene Werte von £. Suchen wir zunächst eine Kurve 
in der x, y-Ebene bei ?! = konst, für deren Punkte gilt, daß &,. = 0 ist, so erhalten wir: 

| — U, Php lV,.: dd (le Di a 
d. h., es entsteht bei Festhalten von ! ein Kreis in der x, y-Ebene mit dem Radius cz t und 
mit Zentrum in U,-t und U„-t. Wir wollen nun zeigen, daß die Potentialgleichung (8) sich 
innerhalb dieses Kreises auf die Laplacegleichung und außerhalb auf die Wellengleichung 
für die Störgrößen u und v zurückführen läßt. 

Durch Rücktransformation der Gl. (8) auf das x’, y'-System bekommen wir Gl. (5). Gl. (5) 


läßt sich durch eine Galileitransformation auf das zur Anströmung ruhende Koordinatensystem 
in die Wellengleichung zurückführen: 


mit 2er Als, Vz ce, 1 We 2176) 
ergibt sich 

Dry + Dyy — Dir = 0. 
Durch Differentiation nach x’’ oder y’ erhalten wir: 


Ur + yy ur —0, ee tyy—e 0.2.22... (8). 
All I 


. ı : 
en 2 u v nur Funktionen von Er und rn innerhalb eines Gebietes der Strömungsebene 
sind, d.h. 


x! %,) z' y' 
(ir). , (7). 


u A a un nn u 
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nennt man dieses Gebiet einen in der Zeit kegeligen Strömungsbereich. Führt man in die Glei- 
chungen (23) die Transformation 


2" =r.tf:.cosp, | y =r"t.sno 
ein, so bekommt man: 
U EP U Se EN (25), 
wobei wir der Einfachheit halber nur die erste Gleichung aufschreiben. Gl. (25) ist äquivalent 
mit der Tschapliginschen Gleichung (siehe etwa [2], S. 204 und 388). 


Die Gl. (25) ist für r > 1 hyperbolisch und für r< 1 elliptisch. Um für r > 1 die Wellen- 
gleichung zu bekommen, setzt man (siehe etwa [3], S. 136): 


& — 0050, 
T, 
wobei Gl. (25) in 
Deu) Lou 
Tu a 2 al ee er (26) 
übergeht. Mit 
1 


ya 
a 1 
15) 
für r< 1 geht Gl. (25) in die Laplace-Gleichung über: 


er N) . . . . . . . . . . . . . . . (27). 


Dieselben Gleichungen wie (26) und (27) gelten auch für die y-Komponente der Störungen. 

Jetzt muß nur gezeigt werden, daß r< 1 dem Gebiet innerhalb und r > 1 dem Gebiet 
außerhalb des Kreises Gl. (22) entspricht. Mit der Einführung von Gl. (7c) und Gl. (7b) in 
Gl. (24) bekommt man: 


ct. co8e + y-sine—t-(a-cose—L, sne&+ U.)=Ca'T'rT.cCosp, 
— Tess, U coset:ta, sine. L we)—=c.tL’r-sneo. 
Nach Elimination von @ entsteht mit Hilfe von Gl. (9): 
Bel STR len ler 
wodurch sich die Behauptung unmittelbar als richtig ergibt, nämlich daß dieser Kreis das 


elliptische vom hyperbolischen Gebiet trennt. Gehen wir nun über zu den „kegeligen Koordi- 
naten““ 


IE EEE IM D) 
j’ y TR ee Ba 


” T 
L =) — 
Co" Co 


wird die Kreisgleichung (22): 5 


Berge Mal een RERECH 


in der X, y-Ebene. Die Charakteristiken innerhalb eines Gebietes in der x, y-Ebene, wo die 
Zustände nur von £ und y abhängen, tangieren den Kreis, gegeben durch (29). Diese Tatsache 
ist sofort mit Hilfe von Gl. (26) und der Transformation zwischen ©, g- Ebene und x, y-Ebene 
zu sehen. 


6. Das Ausgleichsproblem 

Mit den in Abschnitt 3 gegebenen Anfangsbedingungen stromaufwärts und mit den im 
vorhergehenden Abschnitt diskutierten Eigenschaften der Potentialgleichung (8) ist es nun 
möglich, das Aussehen des ganzen Ausbreitungsgebietes in der %, y-Ebene darzustellen. Zuerst 
müssen wir jedoch zeigen, daß wir die Koordinaten X, y verwenden dürfen, d. h., daß der Aus- 
breitungsvorgang der Welle von kegeligem Typus ist. 

Wenn das Problem kegelig ist, muß das Integral Gl. (12) von kegeligem Typus sein. 
Daß dies der Fall ist, kann man auf folgende Weise einsehen. Mit u und v nur als Funktionen 
von X und y läßt sich ® (x, y, c» : ) umschreiben: 


DS. a Rlz.n), 
denn für ®,, ©, und O, entstehen nur Funktionen von & und y: 
i %,=-F&y); 9%, = F,&y), 
una ra 7 Fey. 


oo 


Z. angew. Math. Mech. 


424 Ryhming, Über die instationäre Überschallströmung durch Schaufelgitter ja. 37 Nr. 11/12 Nov./Dez. 1957 
en se 


\ ie Li naar ji Ableitungen, so 
i ir zunächst die Linearkombination Gl. (19) mit den oben gegebenen | Ki 
Rh ir. daß f(x, t) für y = 0 eine Funktion G(&) ist. Damit ist es möglich, kegelige ne 
naten in das Integral Gl. (12) einzuführen, wodurch gezeigt wird, daß das Problem kegelig ist. 
Es ist nun mit I“ ; ; 


ae 7 


und caT; 


EZ 


ls neuen Integrationsvariablen: 
FR ä dE-dr=rT:d(» Tr): dy. 


In Hinblick auf das Vorzeichen gemäß Gl. (19) ist das Störpotential füry > 0 gegeben durch: 
Co 
Dım “= | [&W x 


Dii 


T-d(c„?T): dy 
KT ME Pre A et MM) + uM et p) 

_ Ce > et ELITE 

-= | Ja rien 2) dp (30) 


Störpotential y>0:$;; 


* (Ms, Mn) 


kai] 


Ungs = O2 On; 


Störpotential Y>0 : ®,, 


Bild 7. Konfiguration der Wellenausbreitung 


Hier bedeutet ®,, das Störpotential und G,,(y) die zunächst unbekannte Quellfunktion in den 
in Bild 7 dargestellten verschiedenen Feldern. D,, ist weiter das zu einem beliebigen Aufpunkt 
in Feld 1i gehörige Integrationsgebiet in der y, c„ t-Ebene. Für y< 0 ist ®,, mit der unbekann- 
ten Funktion Gs,(y) gegeben durch: 


(as 1 
= | | Guy) Zr: den): dp en en 
Dei 

Zunächst studieren wir die Anströmbedingung für das %, y-System. Es ist mit Hinsicht auf 
die gegebenen Randbedingungen vorteilhaft, U, so zu wählen, daß u,, = 0 wird. Wir wählen 
also als Linearisierungs-Machzahl die Machzahl, welche in Feld 1, Bild 3a, in der Nähe der 
Schaufel vorherrschend ist. Damit berechnen wir die durch die Welle entstandenen Störungen 
gegenüber der quasistationären Strömung, welche in Abb. 3a gegeben ist. Da im &, y-System v,, 
verschwinden muß, können wir die Randbedingungen up = dp =, Cu = c„ des Ausgleich- 
problems bei y > 0 formulieren. Damit wird mit Gl. (21c) G, = 0, und Ö,; liefert nur einen 
Beitrag in den Gebieten 11, 12 und 13. Infolgedessen ist das Ausbreitungsgebiet der Welle für 
y > O0 völlig bestimmt. Für y< O gilt etwa dasselbe, nur sind die Zustände zz, Ds, Und C,, natür- 

lich von Null verschieden. 
In Bild 7 ist die dimensionslose Geschwindigkeit der Laufschaufel längs der X-Achse, welche 


wir mit 9 bezeichnen wollen, so gewählt, daß die Vorderkante der Schaufel im hyperbolischen 
Gebiet der Welle liegt, d.h. 


r% 1 — 
= @+W)<M—Y1l—-M a ee 
Dies entspricht mit der Randbedingung u, =, = (0, (oc, einer Laufgeschwindigkeit L, 


welche bei den am meisten vorkommenden praktischen Fällen genügend unter c,„ liegt. Auf 
diesen Fall wollen wir uns hier, wie früher bereits erwähnt, beschränken. 


une nn 
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7. Berechnung der Zustände im hyperbolischen Teil der Welle 


Zunächst sollen die Zustände ı,,, Ug1, Di = Va Und Cıı = C,, längs der gestrichelten Linie 
in Abb. 7 berechnet werden. Die Bedingungen v,, = d,, und Cı, = C,, liefern dabei zwei einfache 
Integralgleichungen für die unbekannten Funktionen G,.(y) und G,,(y). Um diese Gleichungen 
aufstellen zu können, müssen wir zuerst die Potentiale ®,, und &,, herleiten. Dabei ist nun eine 
erste Integration von Gl. (30) und (31) über 7 bei konstantem y sofort ausführbar. 


Y 


Bild 8. Integrationsbereich, hyperbolischer Teil der Welle 


Das Integrationsgebiet irgend eines Punktes P (x, y), y > 0 im hyperbolischen Gebiet der 
Welle ist in Bild 8 dargestellt. Die Grenze 7, und , bedeutet wieder das Verschwinden der 
Radikanden in Gl. (30) und (31) und y,,, ist gegeben durch: 


:— (M5—M,)G—M) + W@-M+G MS K 
E In | re 


für yS0. 
Die Integration von Gl. (30) und (31) über 7 liefert dann, wenn man die Grenze 7, und r, 
berücksichtigt: 
Ka ae =. (M,—p):2+ M2y+l11—- ;M—M,M,—p)l:ca1 
| rI8 SEHörr on IM, —y)? + M„—1]?” 


CooTı 


=. Ay) 2 + By) y+ C@y) 1]; 


wo A(py), B(y) und C(y) leicht identifizierbare Funktionen sind. Damit ist unter Berücksichti- 
gung der Grenze y, und y, Gl. (33) 


für y>0: Pi LO) - JAly) -2F-B(p) y -FCp) SH dpr ee A); 
2 
für y< 0: Pdy = F Gay) [Ay 2 + Bo) -y+ CC) I-dp... 2... 5). 


Wenn. wir zunächst die Komponenten v,, und v,, aufschreiben wollen, müssen wir die Rand- 
bedingungen erfüllen, d. h. v,, muß den Wert v,, = O0 und v,, den Wert v,, auf den beiden strom- 
aufwärts liegenden Berandungslinien annehmen: 


dv = dp + (Din), = (did > = U + (Ba)y- 
Mit Gl. (21e) gilt für c,, und c,.: 
1 2 1 1 2 1 
T, (Did = (Pd + a ar Co: M, D B, (Ba), = (Pad + TEN (Co — Co) ° M,' 


Für y = 0 gilt nun v,, = dz, und Cy; = Cgı- Berechnet man die verschiedenen Ableitungen von 
Gl. (34) und (35) und sucht ihre Werte für y gegen Null, so kann man die Bedingungsgleichungen 
Dj = Lg, Und Cıı = Ca, Pilden und bekommt: 


© 


FIIR ve M„M,—i) Gl) + GR) Pr Gyıl®) — Gzıl®) 
= N we Darin om 


0) 


2 


| 


ae Co) "7 


1 [ion Ep) + A| dy 
+1 elıe 2 G®) SE GR) u 
Y 


(M,— ©: +M: 1 
28 


U, 
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Diese beiden Gleichungen haben eine einfache Lösung, indem sich die Funktionen G,,(y) und 
G,,(y) als Konstante erweisen. Führt man nämlich die Integrationen unter der Voraussetzung 
durch, daß G,, und G,, Konstante sind, so erhält man folgendes lineare Gleichungssystem: 


M„'M, 
Don (L— M;) = (Gy + Ga) °° = — Gı + 62; 


2 


1 1 
re mn 
8 


worin o durch Gl. (14) gegeben ist. Die Lösungen von G,, und G5, wird: 


1 2 oe+M,' M rt 
G1 = Do, - vo (1 — My) + Fe (Ce — Co) ° a ... .@6), 
1 2 0 MV 
Gin = 5 m MD + te 5 2. @N. 
Die Zustände sind nun einfach zu berechnen und man bekommt: 
a 1—Mı "oa (2 0 M,M, 
usa een ge Mu 
13 1—M, | Ge — ts MM, 
er Fi De ee EM 68) 
1 1 M„M,, #1 M,(i — M,) 
CH 5 ce eh): Do (Os ie z: Zen £ ee > 
| nn et 
en 


Die Größen mit Index 11 sind wieder im ganzen Gebiet 11 und die mit Index 21 im ganzen 
Gebiet 21 der Welle gültig, wovon man sich leicht überzeugen kann. 

Die Zustände c,ı = Czı geben unmittelbar eine Auffassung vom Druck. Da M,, klein gegen 
1 ist, sind die beiden mit demselben Vorzeichen versehenen Korrekturglieder in c,, klein und der 
Druck ist hauptsächlich durch den Mittelwert des Druckes in Gebiet 10 und 20 gegeben. Es geht 
also eine Verdichtungswelle nach oben in das Gebiet 10 und eine Verdünnungswelle nach unten 
in das Gebiet 20. Da es einen Unterschied zwischen dem Drucksprung cz, — Cı» y < O0 und C,,—Cıo 
y > 0 über den Stoß gibt, erfährt der Stoß im Ursprungspunkt einen Knick. Die nach unten 
vom Ursprungspunkt ausgehende Machlinie stellt dann eine Verdünnungslinie dar, d.h. es geht 
ein Expansionsfächer von dem Ursprung nach unten aus (siehe Bild 10). 

Zunächst setzen wir die Berechnung der Zustände der Welle fort und leiten nun Potential 
und Zustandsgrößen im Gebiet 12 her. Es gilt für irgend einen Punkt P im Gebiet 12: 


PD, = — Gun | 1A) 2 + Bo): y+ Cal: dy — [ G1(y) -[AYy) x + Biy)y + Cpl dp, 


wobei G,, die schon bekannte Quellfunktion über der Strecke 0 bis 9 und G,,(y) die gesuchte 
Funktion über der Schaufel darstellt, siehe Bild 9a und 9b. Mit der Randbedingung v,, = 0 für 


UyyYıp Ep 
6 Gy) 


EN, x 


a 


Bild 9. Integrationsbereich, hyperbolischer Teil der Welle 


y = O und v,,, bekannt auf der von dem Punkt y ausgehenden Ch 


arakteristik, erhält 
Gleichung für G,,(y) auf y = 0: a 


. =0=v, + (B)y; 
oder nach Ausführung der Ableitung: 


7 x 


. 


du —Gu By) dy — Gualy) = Bl) dy + ae, | _ MM nd) 1)=0. 
l—-M [Ym,—2%2ı{M:—1 


ö „ 
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Die Funktion G1(y) erweist sich wieder als Konstante und man bekommt nach Ausführung der 
Integration: 


ee Mm 9m) 


G2 = 
12 Pr me ne 


woi2?=(M, — y)’ 2. M—1l und G,, und v,, durch die Gl. (36) und (38) gegeben sind. Damit 
bekommt man für die Zustände im ganzen Gebiet 12, indem man die Randbedingungen beachtet 


50% j) 
Eee 
Ua = Una + (Da) = 2 en) | 
4 —L 1 
Brsee ME 1 Did + (B,). AO: 
63 
R Bere "Be MM) M,-y)-M4 
ars Den A—M,(M,— 9) 


Für das Gebiet 22 ist: 
7 Yı 
Pr = GC [IA + By + Ci) dp + [Galy) [Az + By) y + Cy)ildy. 
7] 


Mit Hilfe der Randbedingung gilt auf y = 0: 


Up = 0 = vu + (®z.),: 
Daraus ergibt sich: 


D (1 — M,) AcıGa m,[M, 187 PM) 


G pe 
= A+M,„(M,—») Bean 
wodurch die Zustände in 22 leicht berechenbar sind: 
1 —=0, ) 
2 Ir (9 —(C% 0°— M„M 
Se REN _ ge 2 \ | 
22 BL M,M,—») (U: — U) ( n%) N M, > (42) 
1 A—M,„(M,—») (—M,M)(M,—y) | 
ln —— (U — ug): M,— — —— 4 a: m _& 8 
a ae) I mn RT MM MN. 


Damit sind die Zustände im hyperbolischen Gebiet 
der Welle berechnet. Die Machlinien, welche Expan- 
sionsfächern und diejenigen, welche Stößen entspre- 
chen müssen, sind in Bild 10 angedeutet. 

Die verschiedenen Zustände im hyperbolischen 
Gebiet kann man auch mit der Charakteristiken- 
Methode berechnen. Dabei kann man die Verträg- 
lichkeitsbedingungen in der Ackeretschen Form 
schreiben: 

u vtga,. = könst. 2.2.8), 
wo &,, der Neigungswinkel der „Machfront” gegen 
die x-Achse ist, Bild 8. Gl. (43) liegt die Auf- Bild 10. Strömungskonfiguration der Welle. (£- 4») 
fassung zugrunde, daß man eine bewegte Machfront * 
in einem System betrachtet, worin die Mach- 
front ruht. Dabei ändert sich nicht die Form der Ackeretschen Formeln, aber die Größe der 
Machwinkel wird infolge der Bewegung geändert. Um z.B.die vier Zustände 4,1, Ugı, d, und Cyı 
in den Gebieten 11 und 21 zu berechnen, braucht man vier Gleichungen. Zwei Gleichungen, 
welche die Zustände zwischen 10 und 11 einerseits und zwischen 20 und 21 andererseits auf 
rechts- und linksläufigen Machlinien verknüpfen, bekommt man durch die Gl. (43). Dabei ist 
&, gegeben durch: 


n. + M,„ M, + [07 
u Tin 
Zwei weitere Gleichungen erhält man aus der Energiegleichung in linearer Form, welche gegeben 
ist durch Gl. (21b) mit ©, = konst. Man benutzt diese Gleichung längs einer Stromlinie, 


weiche die Zustände 10 mit 11 und die Zustände 20 mit 21 verbindet. 


28* 
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ee RER NE I en 
8. Das elliptische Gebiet der Welle 5 h 
Die Berechnung der Zustände im hyperbolischen Teil der a im an: er 
i ößer ierigkeiten durchgeführt werden. Im elliptischen - 
Abschnitt ohne größere Schwierigkeiten rchgefü N 
j hebliche mathematische Schwierigkeiten, die Zustä i 
er Ba Schwierigkeiten hinweisen und danach nur die Zustände auf der Berandung 
des elliptischen Teils der Welle diskutieren. 


Bild 11. Integrationsbereich, elliptischer Teil der Welle. 


j jet i j Icher im elliptischen Teil der 
Integrationsgebiet irgend eines Aufpunktes (y > 0), we 1 
he Sn Bild 11 dargestellt. Die Grenzen des zuerst auszuführenden Integrals Gl. (30) 
sind nun anders und es ist: 


Coo Ta 
| en = H«(z, y, t,y), 

. . 
wobei der Einfachheit halber das Integrationsergebnis nur als eine Funktion H dargestellt ist 
Für das Potential ®,, bekommen wir dann folgenden Ausdruck: 


e a ME 
Gel [Haptan: | Mat | Cu Hapl.. . 
i 0 a 


IT 


[2 M,-Yı—M2 
worin G,, durch Gl. (36) und G,, durch Gl. (39) gegeben sind. Um eine Gleichung für G,, zu bekom- 
men, hat man die Randbedingungen zu erfüllen. Es gilt auf y = 0: 


VG =Vlart (D,;), =, oder (B,3), =. 


Diese letzte Bedingung führt zu einer Integralgleichung für die unbekannte Kernfunktion G,s(y). 
Wir bekommen aber hier ein sehr schwieriges Problem, und es ist nicht einfach, diese Integral- 
gleichung zu lösen. Be 

Wir wollen nun einige Aussagen über die Zustandsgrößen auf der Berandung des elliptischen 
Gebietes machen. Wir studieren das Problem zunächst in den Koordinaten r und p mit Gl. (24). 
Es ist u, =D, = Ugg = Vyg — 0, weil die Zustände stückweise konstant sind. In einer Umgebung 
der Strömungsebene, wo u,, — 0 ist, läßt sich Gl. (25) wie folgt umschreiben: 


((—r)- (u, +u]—r-u,=0 

oder RE NE 
(1 —r)- (ru), —- I (ru)=0 

Die Lösung dieser Differentialgleichung ist: 


Eos m a ee ee 


1—r 


wo C, (9), die Integrationskonstante, eine Funktion von p sein kann. Die Gl. (46) kann noch 
einmal integriert werden mit dem Ergebnis 
1+ Ir 
r 


a On 


Auf dem Kreisbogen r = 1 hat also u einen unendlichen Differentialkoeffizienten, aber u selbst 
ist über den Kreis hinweg stetig. Das ist ein ganz allgemeines Ergebnis bei kegeligen Strömungen, 
siehe etwa [4], Seite 55, Bild 30. 

Eine Ausnahme sind die T angierungspunkte der Charakteristiken, welche Punkte durch 


einen Sprung in den Randwerten selbst stark singulär werden. In diesen Punkten nehmen daher 
die Zustände mehrfache Werte an. 


Um die Geschwindigkeitsverteilun 
könnte man die elliptische Randwertau 


u = — C1(p) -log 


g innerhalb des elliptischen Gebietes zu berechnen, 
fgabe mit einer der Gleichungen (25) oder (27) angreifen. 


4 
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Dies gelingt allerdings nur für die v-Komponente, weil nur die v-Komponente auf der ganzen 
Berandung, Kreisbogen und Schaufel, bekannt ist (vgl. Dirichlets-Problem). 


9. Einige Gesichtspunkte auf der Ausbreitung der Wellen in einem Expansionsgitter 


A Durch die Auffassung des Ausbreitungsvorganges eines Stoßes in der Strömung, welche 
wir in den vorhergehenden Abschnitten bekommen haben, wollen wir nun die Folgerungen der 
Wellenausbreitung auf die Gesamtströmung in einem Expansionsgitter diskutieren. Weil die 
Strömung mit einer Rückwirkung verbunden ist, kommt es allerdings schnell zu komplizierten 
Reflexionsvorgängen des Stoßes und des durch das Vorrücken eines Stoßes vor die Schaufel 
hervorgerufenen elliptischen Strömungsgebietes in dem zunächst betrachteten Rückwirkungs- 
gebiet des Gitters. Eine allgemeine Behandlung dieser verschiedenen Erscheinungen ist mit 
unserer Theorie schwierig zu geben, weil die Randbedingungen auf mehreren, teils bewegten, teils 
stillstehenden Konturen zu erfüllen sind. Wir beschränken uns daher hier auf eine einfache 
en Diskussion einiger Einzelheiten, welche für das Strömungsbild im ganzen wichtig 
sind. i 

In dem betrachteten Expansionsgitter können nur Wellen ausgelöst werden bei kleinen 
Verhältnissen b/l, was aus einem Studium der Bildserie 11—14 [1] hervorgeht. Innerhalb des 
größeren Teiles einer Teilungsstrecke ist die Grundströmung quasistationär. Von der Stellung, 
welche etwa durch Bild 12, [1] gegeben ist, schieben sich allerdings Störungen von der Ober- und 
Unterseite der Laufschaufel vor deren Vorderkante. Dadurch werden wohl wesentlich instationäre 
Vorgänge ausgelöst von dem Typus, wie er in den vorangegangenen Abschnitten dieser Arbeit 
besprochen wurde. Die Störungen aber, welche in diesem ersten Stadium vor die Vorderkante 


Bild 12a. Schematischer Strömungsvorgang im Gitter. Zeitpunkt i = t, + 4t 


der Laufschaufel reichen, sind sehr schwach, so daß diese wesentlich instationären Vorgänge 
noch nicht in Betracht gezogen zu werden brauchen. Aus den zwei letzten Bildern 13—14, 
[1] ist jedoch zu sehen, daß es zwei Stellungen gibt, bei denen ein Stoß erstens von der Unterseite 
der Laufschaufel her und zweitens von der Oberseite der Laufschaufel her vor die Vorderkante 
der Laufschaufel treten kann. Mit den Laufgeschwindigkeiten, auf die wir uns hier beschränken, 
werden zwei Wellen bald hintereinander während einer Teilungsstrecke ausgelöst. 

Der schiefe Stoß erreicht im Gitter bei b/! = 0,16 die Position, bei der eine erste neue Welle 
in der betrachteten Teilungsperiode ausgelöst wird. Im nächsten Moment ergeben sich Reflexions- 
vorgänge des Stoßes und des elliptischen Gebietes mit der Unterseite der Leitschaufel Bild 12a. 
Durch die Reflexion des Stoßes an der Unterseite verliert der Stoß die erste Beziehung zum 
elliptischen Gebiet, welche er im ersten Augenblick bei Auslösung der Welle hätte. Die Aus- 
breitung des Stoßes durch das elliptische Gebiet geschieht daher in der Fortsetzung anders als 
vorher. Man erkennt durch einen Vergleich der Strömung in Bild 10 mit der der quasistationären 
Konfiguration in Bild 14 [1], daß der instationär-kegelige Ausbreitungsvorgang im hyper- 
bolischen Gebiet mit der quasistationären Lösung praktisch übereinstimmt. Der quasistationäre 
Zustand breitet sich also mit der Geschwindigkeit des linken Randes des elliptischen Gebietes aus. 
Die Zeit, die der Stoß braucht, um die quasistationäre Konfiguration zu erreichen, ist sehr kurz. 
Die Laufschaufel hat während dieser Zeit nur eine kleine Strecke zurückgelegt. 


Die Geschwindigkeitsschwankungen im elliptischen Gebiet behalten jedoch ihre Größen- 
ordnung, solange eine dafür entsprechende Energiemenge der Welle zugeführt wird. Die Energie 
erhält die Welle von den Quellen, welche auf der gestrichelten Linie vor der Laufschaufelvorder- 
kante angenommen sind (Bild 3b). Die Quellen haben hier, wie früher gezeigt ist, in der kegeligen 
Darstellung eine konstante Stärke. In der gewöhnlichen physikalischen x, y, t-Darstellung 
wachsen sie proportional mit der Zeit, d. h., die Störungen bleiben mit zunehmender Ausdehnung 
des Störfeldes unverändert, solange nur 9 Gl. (32) konstant ist. Diese Bedingung ist bei der 
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ü i ö j ährend einer gewissen Zeit 
lle anfangs etwa erfüllt. Die Quellen können jedoch nur währe 
Se a ausblasen, und zwar solange die erste BER ne 
i ö i ö i enn nac 
ält. Die erste kegelige Strömung wird abgelöst durch eine zweite, l ; 
= Hähtschaufel Ene a; Strecke zurückgelegt, daß der an der Oberseite der Schaufel reflek 


tierte Stoß eine zweite Welle auslösen kann. Um die zweite Welle auslösen zu können, muß die 


indiokeit einen solchen Wert haben, daß die quasistationäre Konfiguration, gegeben 
Fer ner Zeit hat, sich etwas auszubilden, ehe die Tale er 
aktuellen b/l-Wert für die zweite Auslösung einer Welle passiert hat. Bei den Lau Br win i8- 
keiten, welche wir hier betrachten, ist diese Bedingung erfüllt. Die zweite Welle wird dann a 
gelöst und damit fängt eine zweite kegelige Strömung an sich auszubilden. Die BEnNERe 2 
sind nun nur verantwortlich für die neue Welle. Da aber die zweite Welle im Gegensa zZ fü er 
ersten von dem Stoß an der Oberseite her ansBoe* wird, haben die Geschwindigkeitsschwan- 

in dieser Welle entgegengesetztes Vorzeichen. h 

San = Zeit es die erste Welle mit der Strömung weg und der Radius des Stör- 
gebietes wächst proportional zur Zeit. Die Welle bekommt aber keine Energie mehr und sie 
breitet sich daher mit gleichbleibendem Gesamtimpuls aus. Die Geschwindigkeitsschwankungen 
gleichen sich schnell aus (vgl. [2] Seite 86—89). In Abb. 12b ist ein Beispiel gezeigt, bei welchem 


Bild 12b. Schematischer Strömungsvorgang im Gitter, Zeitpunkti = i, + 54. Erste Welle nicht kegelig, zweite noch kegelig 


die erste Welle nicht, die zweite aber noch ihren kegeligen Charakter hat. Auch die zweite kegelige 
Strömung verliert aber bald diesen Charakter und das Resultat des Ausbreitungsvorganges werden 
zwei Wellen mit gleichbleibendem Gesamtimpuls, welche mit der Strömung wegschwimmen. Die 
Störgebiete der beiden Wellen überdecken sich dabei teilweise, und wegen entgegengesetzter 
Wirkung wird die Gesamtstörung der Wellen in diesem Gebiet außerordentlich gering. 

An der Wellenfront fangen aber die Störgrößen an zu wachsen mit einem unendlichen 
Differentialquotienten, Gl. (46). Eine neue Welle wird ausgelöst, sobald die Wellenfront die 
Vorderkante einer Laufschaufel einholt. Die Störgrößen in dieser letzten ausgelösten Welle sind 
nun sehr klein, weil die Störungen an der Wellenfront sehr klein sind. 

Die Strömung im angenommenen Gitter ist aber mit Rückwirkung verbunden und daher 
können die ausgelösten Wellen auf den Strömungszustand in den nächsten Teilungsperioden 
einwirken. Diese Einwirkung wird jedoch sehr gering wegen des Abklingens der Zustände in den 
Wellen und der zahlreichen verschiedenen Reflexionsvorgänge und vor allem, weil die Störung 
an der Wellenfront klein ist. 

Die reflektierten Teile der Wellen können im Gitter nur eine beschränkte Anzahl von 


Teilungsperioden vorwärtslaufen. Der Expansionswinkel des Expansionsfächers am Ende der: 


Leitschaufel nimmt kurz bevor neue Wellen im Gitter ausgelöst werden kräftig zu, so daß das 
Rückwirkungsgebiet ganz überdeckt ist. Da eine Störung nicht über die letzte Machlinie dieses 
Expansionsfächers treten kann, werden im ganzen Gitter sämtliche Störungen in das Laufrad 
„hineingeschoben“, ehe neue Störungen ausgelöst werden. 

Das Ergebnis dieser Überlegungen ist, daß die quasistationäre Konfiguration des Strömungs- 
feldes nur während einer kurzen Zeitdauer von diesen Wellen beeinflußt wird. Die instationären 
Effekte des Ausbreitungsvorganges der Wellen haben zwar Einfluß auf die ganze Strömung im Lauf- 
rad, aber durch das Abklingen und durch die Reflexionen wird dieser Einfluß im ganzen gering. 
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10. Zusammenfassung 


Die Lösung des Strömungsproblems durch ein Schaufelgitter mit Rückwirkung wurde im 
stationären Fall in [1] gegeben. Wenn man die Laufgeschwindigkeit L des Laufrades berück- 


_ sichtigt, zeigt sich für ziemlich kleine L-Werte im Vergleich zur Schallgeschwindigkeit, daß die 


stationäre Lösung im wesentlichen richtig ist. 

Es gibt eine typisch instationäre Abweichung von der quasistationären Lösung, welche in 
dieser Arbeit hauptsächlich behandelt ist. Die Ursache dieser Erscheinung ist, daß ein Stoß vor 
die Laufschaufelvorderkante gelangt. Der Einfluß dieser instationären Erscheinung auf die 
quasistationäre Lösung ist nicht wesentlich. i 

Diese Arbeit wurde unter der Leitung von Herrn Professor Dr. Klaus Oswatitsch ab- 
gefaßt, der sie immer mit großem Interesse verfolgt und viele wertvolle Vorschläge sowohl bei 
der Planung als auch bei der Durchführung gemacht hat. Ich möchte Herrn Professor Oswatitsch 
für alle-seine Hilfe an dieser Stelle herzlichst danken. 
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Eingegangen vom 28. Januar 1957. 


Einfluß einer unendlichen Reihe gleicher Kreislöcher 


auf die Durchbiegung einer dünnen Platte 
Von O. Tamate in Sendai (Japan) 


Die Aufgabe der Abschätzung des Einflusses einer unendlichen Reihe gleicher Kreislöcher auf die Spannung 
in einer unbegrenzten, biegebeanspruchten Platte wird auf Grund der Kirchhoffschen Biegungstheorie und 
der komplexen Darstellungsmethode von Muskhelishvili behandelt. Zur Bestimmung der Koeffizienten, die in 
der Durchbiegungsfunktion der Platte enthalten sind, wird die Störungsrechnung angewandt. Numerische 
Ergebnisse für einige speziellen Fälle werden mitgeteilt. 

The Kirchhoff theory of bending and the method of complex representation due to Muskhelishvili are 
used to obtain an estimate of the influence which an infinite series of equal circular holes will have on the stress 
distribution in a plate subjected to bending moments. In order to determine the coefficients included in the de- 
flection function a perturbation theory is used. For some special cases numerical results are given. 


Sur la base de la theorie de fleche d’apres Kirchhoff et de la methode complexe de description selon 
Muskhelishvili l’&valuation de Vinfluence d’une serie infinie de trous circulaires et egaux sur la tension 
dans une plaque non-limitee et flechie est traitee. Afin de determiner les coefficvents compris dans la fonction de 
fleche de la plaque la methode de perturbation est appliquee. Des resultats numeriques pour quelques cas speciaux 
sont communiques. 


3anaya OleHKU BJIUAHNA ÖeCKOHEYHOTO PANA ONHHAKOBEIX KPyTJIbIX OTBEPCTHÄ Ha HanprRe- 
HUE B HEOTPaAHNYeHHON, HATPy}KeHHOH Ha H3TUO MNACTHHKE, HCCJIeNYeTCA HA OCHOBAHHH TeOPHH 
usruba Kupxroha u MerTona N306pa;keHnA B KOMIWIEKCHOM Bune Mycxennmsnam. as 
oNpeeleHuUA KOIPPHILMEHTOB, CONeP>RAMHXCH B DYHRIIMM MH3TU0a INTACTUHKU, IIDHMEHAECTCH 
BbIYMCJIeHMEe BO3MYINeHUÜ. „II HEeKOTOPbIX YACTHbIX CIIyyaeB MaloTcA YUMCJIEeHHBIe Pe3syJIb- 


TATbI. 
I. Einleitung 

Bekanntlich hat das Vorhandensein eines Loches in einer biegebelasteten Platte eine große 
Wirkung auf die maximale Spannung, und ein große Anzahl theoretischer Untersuchungen!) 
wurde durchgeführt, um den Einfluß eines Loches in verschiedenen Formen zahlenmäßig zu 
erfassen. In den meisten Fällen dieser Arbeiten wurde aber die Platte als unbegrenzt in zwei 
Dimensionen und als geschwächt durch ein einzelnes Loch angenommen. Es ist klar, daß die so 
gewonnenen theoretischen Lösungen ein um so genaueres Bild von den Spannungen in einer 
Platte von endlicher Ausdehnung liefern, die durch das eine oder andere Loch geschwächt ist, 
je kleiner die Abmessungen des Loches im Vergleich zu denjenigen der Platte sind. 

Folglich ist es wichtig in Ingenieur-Praxis, daß die Einwirkung der benachbarten Löcher 
oder Grenzen auf die größte Spannung numerisch abgeschätzt wird. In den vorhergehenden 
Arbeiten [1], [2]?) hat der Verfasser die Aufgabe der Bestimmung der Spannungserhöhung neben 
einem kreisförmigen Loch für den unendlich langen durch Randmomente verbogenen Streifen 
endlicher Breite betrachtet. 

Wir suchen nun zahlenmäßig den Einfluß von einer unendlichen Reihe gleicher Kreislöcher 
auf die Spannung in der unbegrenzten biegebeanspruchten Platte zu erfassen. Hierbei wird die 


1) Eine kurze Übersicht der vorhergehenden Arbeit ist inO. Tamate [2] gegeben. Siehe auch G.N.Sawin[3]. 
2) Die Zahlen in eckigen Klammern verweisen auf das Literaturverzeichnis am Schluß der Arbeit. 
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be auf Grund der Kirchhoffschen Biegungstheorie elastischer Platten [4] und der kom- 
En Darstellungsmethode von Muskhelishvili [5] behandelt. Zur Bestimmung der Koeffi- 
zienten, die in der Durchbiegungsfunktion enthalten sind, wollen wir die Störungsrechnung an- 
wenden. 
II. Problemstellung 
Wir nehmen eine dünne Platte aus elastischem isotropem Material mit einer unbegrenzten 
Reihe kongruenter kreisförmiger Löcher vom Radius a an, die so verteilt sind, daß jedes folgende 
Loch durch Verschieben des vorhergehenden um die Strecke b längs einer Gerade erhalten wird. 
Wir legen den Ursprung der Koordinaten in den Mittelpunkt 0 irgendeines (als Bezugsloch 
angenommen) der Löcher und die x-Achse längs der durch die Mittelpunkte der Löcher laufenden 
Gerade (Bild 1). In bezug auf diese Achsen können wir die komplexe Koordinate z in folgender 
Weise definieren: 
z=r+iy =TIo HA a ae (2.3): 
Ferner werden wir den zusätzlichen Koordinatenachsen durch den Mittelpunkt jedes aller übrigen 
Löcher parallel zur x- und y-Achse legen. Wenn wir z. B. den Mittelpunkt des m-ten Loches als 


Bild 1. Koordinatensystem. 


sammenfällt und die y„-Achse nach unten gerichtet ist (siehe Bild 1). Damit werden wir auch die 
komplexe Köordinate z,, gebrauchen: 


zu, = Hiyn=7n (6020 SSLSn GO, ae ee 
Wir beziehen alle Maße auf den Lochabstand b, führen also die dimensionslosen Koordinaten ein: 


=zib, E=vlb, n=yib, oe=r/b, N. (2.3) 
Im Zuld, Em End, Nm Ymlds Om = mid; = alb 
Beide Koordinatensysteme werden nach der Formel 
Ze zz enbi cl, ein mel) 


verbunden. 

Nun nehmen wir an, daß die Platte der Beschränkung der Kirchhoffschen Biegungstheorie 
elastischer Platten folgt. Danach wird die Dicke der Platte als klein gegenüber den linearen 
Oberflächenabmessungen vorausgesetzt und auch die Durchbiegung der Platte als gering im 
Verhältnis zu ihrer Dicke. Dann ist die Differentialgleichung für die Durchbiegung w(z, y) einer 
auf ihrer Oberfläche kraftfreien, nur durch Randmomente und Randkräfte gebogenen Platte: 


Aa (2.5) 
Dx3. ..Qy2 bu 

Es ist bekannt, daß für die biharmonische Funktion w stets folgende Darstellung gilt: 
oz, y) = Re [Zo@) -Aala)lo2 2 Ser RT 


Hierbei ist Re[- - -] der Realteil des Ausdruckes, der in der eckigen Klammer steht; p(z) und x(z) 
sind analytische Funktionen der komplexen Variablen z=x + iy; ferner it z= x — iy. Die 
Spannungen werden durch die in Richtung der Dickenerstreckung zusammengefaßten Kräfte 
und Momente dargestellt. Für diese Größen gelten die folgenden Beziehungen: 


M,+M,=—4D(l +»)Re[p'@)], 
M,—M,+2iH,,=—2DA—)Ep'@&)+y@], ya) =dyaldz,t.. &7) 
N,—iN,=—4Dg"(), 


wobei M,, M, die Biegemomente, H,, das Torsionsmoment und N,und N, die Randkräfte sind. 
Jede der Größen wird auf die Längeneinheit bezogen und die x- und y-Suffixe bedeuten die 
Wirkung auf den Rand senkrecht zur x- bzw. zur y-Achse. Unter D und » sind hierbei die Platten- 
Steifigkeit bzw. die Querkontraktionszahl zu verstehen. 

Mit Hilfe der Funktion z = w(f) = w(rei®) gehen wir von den kartesischen Koordinaten x, y 
zu den orthogonalen krummlinigen Koordinaten r — const, © = const über. Alsdann ergibt 


die Formel (2.7) für die Komponenten M,, Me, H,o, N, und N. der Randmomente bzw. der 
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Randkräfte in bezug auf die Achsen (r), (©) 
M,+Mo=—4D(l +»)Re [O0], 
M,— Me+2iH,e = Na _ 


ZOO +7. 98 


iD Dt), 


r| 20 RxXo] 
wobei zur Abkürzung 


oe) = Pod] = PM, ya) = ylod] = yı ld, 


A) AU) (2.9) 
on -99 eg 
a. > un 
gesetzt ist. 
Bekanntlich lauten die Randbedingungen für freien Plattenrand, wie folgt: 
u® 
Mt, = N eo N (2.10), 


wobei M,„, H,„ und N, das Biege- bzw. das Torsionsmoment und die Randkraft je Längeneinheit 
des Plattenrandes sind, dessen äußere Normale n ist und dessen Bogenlänge s. ‘Die Randbe- 
dingungen (2.10) kann man leicht durch die Funktionen o(z) und y(z) ausdrücken: 


1 Sanlorznune 
— pe) rzy ed) try =iCz+G, TERN! 
Dabei ist C eine reelle, C, eine komplexe Konstante mit 
1— 
en era ar: roh Ai 
(2.12) 


Mit Hilfe von diesen Beziehungen und Bezeichnungen bestimmen wir den Spannungszu- 
stand in der Umgebung der Löcher unter der Bedingung, daß einerseits auf den Rand der Löcher 
keine äußeren Kräfte wirken, andererseits die Spannungen in großer Entfernung von den Löchern 
den folgenden Werten entsprechen: 

a) Reine zylindrische Durchbiegung um y-Achse 


Me9=M,, M@=,M, HD9=ND=NM=0..... (21a). 
Folglich gilt für die Durchbiegung der Mittelfläche der durch keine Löcher geschwächten Platte 
M, b? 
W, rn, nm N 5 ee er 2.13b): 
Die Funktionen 9 und », die der Funktion w, (2.13b) entsprechen, haben die Form: 
oz) = ee yNz) = Ei Bull Bebsactl Se Hola ch 


b) Reine zylindrische Durchbiegung um x-Achse 
M»=M, M®=»,M, H9=-ND=-NM=0 ;) 


M,b? 2, 
Ta re A: 
Mb, Mb, 
c) Reine einachsige Biegung um y-Achse 
M®=M,, Ma»=-HY)=N®=-N®=0, ) 
M, b: 
Bar one: eo! 
M,b M,b - | 
(0) Se ee (Oz) — 0 Be 
ee man ae Ban 


d) Reine einachsige Biegung um x-Achse 
Mo'=M, MM=H%=- NO=NM=0, 


M, b? 


| 
= apa Bere); N 12.10). 
Mb ee | 
( 0 E (EYE ge ee 
ee ATIDU—n° } 


Z. angew. Math. Mech. . 


434 'Tamate, Einfluß einer unendlichen Reihe gleicher Kreislöcher 134.37 Nr. 11/12 Nov./Dez. 1957 


Ferner nehmen wir an, daß die Platte keine Querlast auf ihrer Oberfläche trägt. Es ist klar, daß 
jede der Durchbiegungen (2.13), (2.14), (2.15) und (2.16) den gleichmäßigen Spannungszustand an 
den Rand eines jeden Loches liefert. 

Zur Erfüllung der den Rändern der Löcher vorgeschriebenen Randbedingungen erfordert es 
noch die Zusatzdurchbiegung. Die Zusatzfunktion hat dann der biharmonischen Gleichung (2.5) 
zu genügen. Mit anderen Worten, können wir die Zusatzdurchbiegung nach Gleichung (2.6) 
durch zwei komplexe Durchbiegungsfunktioneır darstellen. 

Wegen der periodischen Anordnung der Löcher und der Gleichmäßigkeit des Grundspan- 
nungszustandes, müssen die komplexe zusätzliche Durchbiegungsfunktionen auch gleiche Formen 
für alle Löcher haben. Nun drücken wir die Zusatzfunktionen, die den Einfluß des Vorhanden- 
seins des m-ten Loches charakterisieren, in der Form i 


NZ 


aus. Wenn wir diese Zusatzfunktionen insgesamt addieren, in denen wir z,, durch z— mb er- 
setzen, so können wir die dem ergänzenden Spannungszustande entsprechende Durchbiegung in 
der Umgebung des als Bezugsloch angenommenen Loches in folgender Form darstellen: 


w = Reg) +) + I E-mb)g® @—mb)+ I W@—mb)] (2.18), 
M= — © m= —o 
wobei x®’(z) = [ y®(z) dz, und $’ die Summe ohne m = 0 zeigt. Dieser Ansatz ist fast analog dem 
vonH. Saito [6] zur Lösung der speziellen Aufgaben der ebenen Elastizitätstheorie angewandten 
Verfahren. 
Die gesuchten Funktionen o® und y® sind nun so dadurch zu bestimmen, daß die Rand- 
bedingungen an den Rändern aller Löcher erfüllt sind. Wegen der Symmetrie des Problems 
benötigen wir jetzt nur die folgende Beziehung an Stelle von der Randbedingung (2.11): 


tür © doll =A0,0 e&), 


1 = ea ee = en. —. 
|#9o + yo) + I 9 @—m ++ + I @—-mb)gW @—mb) 


m=— co mMm=— © 


+ yo) Hyde + > yOE—mD=0 ee Ben 


m=—&© 


III. Reine zylindrische Durchbiegung um y-Achse 


Nun betrachten wir die Aufgabe der wirklichen Bestimmung der unbekannten Funktionen 
p(z) und g®(z) für den Fall der reinen zylindrischen Durchbiegung um die y-Achse. 

Es ist klar, daß die Funktionen 9 und y® in den unendlich entfernten Teilen der Platte 
einen Nullspannungszustand ergeben muß und auch die eindeutigen Komponenten der Ver- 
schiebung im ganzen Bereich der Platte. Beachten wir dies sowie die Symmetriebedingung des 
Problems, so stellen wir die Funktionen in Form 


Mb e, M,b|a < d 
(1) Ben 0 2; v28 a 0 ) SS 2 - 
Y (z) D — 2e=1’ yD(z) — D E + aan Et (3.1) 
dar, wobei dy,d,,...; 6 &,... reelle Konstante sind. Setzen wir die Funktionen y%(z); y%(z) 


(2.13c) und g(z); y®(z) (3.1) in die Bedingung (2.19) ein und führen die Beziehungen 
, = er ER ann, 62 n-+1 { ) 
n=0 


Er l I 2 
> (25 een = Dan ET a 


n=0 


EIERN IN 9 zn, 2 —1 had 
= Den Lu 
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mit 
1 Pe | 
2 er ge 
er lanı = al Ban Jana 
FERH 9 a s (3:3): 
2 A n FA — 
es | On Jun: SR yanrs >| In+i Jane 
0 a a re | 
m=1 e i 


ein, so erhalten wir 


1 1x — = 
eh Ar RE: Da @s—1) g-(2s—ı1) Die N arg, ABuHı orn+l 
de & z sl n=0 


De 7} Hi > (2 Sn 1) &,4 @:-D o2s+1 en Se. 3,2 n u o 2r—)) 
s=1 s=1 n=1 


neo n + 1) 2ny, , Aantl gen) — rel 
s=1 s=1 


ee > 25de,A%s+D g2s+1 + d, >, 2 n2n x, Jan o “n—)) 
s=1 


n=1 
+ > ds, Se er a ee, 
s=1 n=1 


Da diese Beziehung für jedes o gelten muß, können wir die Koeffizienten verschiedener Potenzen 
von o gleich Null setzen. Sodann bekommen wir folgendes Gleichungssystem zur Bestimmung 
der Koeffizienten der Funktionen g und y®. 


Is l 1 o 
—— u Haug‘ rt u 0 2 I. == 
ale ı) 7 % a 122 Y2s 025 0, 


1 oo 
= = Öl,n > N &n+2n Az De dz; | 
s 2 3.0), 
+ 2n Ant Damp, — (2 + 1) Arrtl Am, nl, je 1) 
s=1 s=l 
1 oo 
— 2nıertd) ,a„— (2n — erde, + — Amt Moe el | 
& s=1 
wobei ö,,„ Kroneckersches ö ist. Durch einfache Rechnung überzeugen wir uns, daß 
1RURT 1 „<< } 
=—(_—1)2?+[— —1)2 O4 65 > 
a 
BE Sn, * } 
dan mel d1,n — I, ""Ies eitn da; 
ri (EAN 
& Br | 
Ba“. ge (?% = ie ARKo, ) a 24% Con, (n > 1) 
s=1 
ge 3% ee N eiin—2d, + Dir, & — 2npy, E22] At" 2 e,, | 
2 s=0 s=1 
wobei zur Abkürzung folgende Bezeichnungen eingeführt sind: 
En 1 
2nn,, —_ (2 ne 1) ZNyn, a anj., = (2 no 1) En, 5 ANe — a | 
RHRRIR 


92 [3 1 6) 
a  — (an a) : 


9 ; p c 29 
anf, =9n anyo; 2 2N]]9% —9Pn?n er Anne, &z (2 n + 1) nn 
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i ? BE 2 den 

Zur Lösung des Gleichungssystems (3.7) können wir die Störungsrechnung anwenden, 
in der A und & die Störungsparameter sind. Da das Verhältnis a/b =) als klein angenommen 
ist, so mögen wir die Ansätze er 


0 © fi 
' 2r+2 
de 23 a) EP 72 + 2 a W ep Z2rt+2, 


p——1 


r-+1 
den = VO En 2 N 
ED 1) 
Be VUN. z j Wen, ep artan—2 
r=1 9= 


machen. Setzen wir diesen Ausdruck in die Gl. (3.7) ein und vergleichen die Koeffizienten der 
gleichen Potenzen von A und e, so erhalten wir 


ii 1 
a0) — 7% (Val —_ — zn x 
' 1 ie er 
VA) — Z DEN — =; 
m na nn; 
; s=1 
( ee e 
de) — — 2m], @—1)e(0) d,r — Pr, e-D d® — 2j, P+De | >) 
: t Tz 
En 2nK, @—1e(0)] Ö2, a 2nn, u 07, u) —y 2ni,, we = +1) 
s=1 
W 
ar [2rh,, PVae—2 —2nj,, (BHDER—2) — Ink, P—1) er—29)], (3.10b), 
s= 
ni 
LY sn 2317, P—de®) Ö1,r u [2° N Fer zu P—1el] O8, r ( = ) 
t 
+ ur ed) a 33 20, Wi Zar 
s=1 
a 
FL er ee, ) nl 2ny,, eDeeze] ) 
wobei die LimiteZ und !’ der Summe folgende Werte haben: 
r/2 r2) — 1 erade 
— / ja wenn r=J® 
(r— 1)/2, (r — 1)/2, ungerade. 
Also im allgemeinen lauten folgende Beziehungen: 
Ydar—ı) — 
% ner ren een 


gen = Veen = 0 ,J 


Mit diesen Gleichungen kann man die Koeffizienten @d(), @e) der Reihen (3.9) nacheinander 
bestimmen. 


Die vorhergehende Lösung ist nur formal und um die Gesetzmäßigkeit der Lösung aufzu- 
stellen ist es nötig, daß die Konvergenz der Reihe (3.9) bewiesen wird. Der direkte Beweis ist 
ziemlich schwer, aber das Gleichungssystem (3.7) kann auch durch sukzessive Approximation 
gelöst werden und dabei wird es leicht auf ganz ähnliche Weise wie von R. C. J. Howland [7] 
erwiesen, daß die Lösung in vorliegendem Fall für mindestens A = a/b < 0,35 konvergiert. 
Ferner zeigt das Ergebnis der numerischen Rechnungen nach den Formeln (3.10), daß die 
Reihe (3.9) bis zu dem achten Gliede einschließlich praktisch befriedigende Konvergenz für 
) = 0,3 ergibt. 

Die Werte der Koeffizienten, die in der Reihe (3.9) enthalten sind, berechneten wir nach 
Formeln (3.10) bis zu r = 8, aber der Platz erlaubt uns nicht die Werte in Tabellen beizu- 
ügen. 

Mit den Koeffizienten kann man leicht die Durchbiegungs- und Spannungsverteilung in 
der Umgebung eines jeden der Löcher erhalten. Bild 2 zeigt die Biegemomentenverteilung 


neben dem Loch für » = 0,3 und A = a/b = 0,25; 0,15. Der größte Wert von Biegemomente 
in der Platte wird beir=a, 9 = + 90° erreicht. 
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&) 
Z 
Io ex 
y & Bild 2. Biegemomentenverteilung neben dem Loch 
Ss für reine zylindrische Durchbiegung um y-Achse 
SZ” er (Querkontraktionszahl » = 0,3). 
| 06 (a) M, und M, längs z-Achse. 
ZN (b) M, und M, y längs Y-Achse, 
04 Ge (ec) Mg am Lochrand. 
02 My I: a/b = 0,25. II: a/b = 0,15. III: a/b —0. 
Fre 
Q 01 02 03 04 05 
— yyb 
(b) 


Da M,= 0 am Lochrand, so ergeben sich die Momente M, an demselben Rand: 
[Mel-=-a/Mo = —4D(1 +») Re [p'@)]--a/Mo 


ns 1— 4m, Hader DE @n+ men 23} os2n 9| (3.11), 


2n 
1 Hr n=1 k 
wobei 
x 
Na = 3.293} e25 - ae (3.12). 
s=1 
Mit diesen entwickeln wir wieder das größte Biegemoment in Reihe nach A und e und erhalten 
o r+1 
M. 1—e — 
mar =27 1+2e+ > > Mar, p e? A?" ne 
My, +e ee 
wobei 
r 
> — rn+1 
Mao =ALZ (1 (2n 1) (M td + Om), 
n=1 
r 
M. —n+1 1)me—n+1 
MR ALM On 1) Br + On 
n=1 
02 r—p+l .. . (8.14). 
—n u 
My» =4| 35 1" (@n — 1) @dt+D + Omi) 
n=1l 
— +2 (9 r = “r 
Tr IT ae DE U FE Ze 1) 
M;,;, 41 E — 4 +Deh) u 
Tafel1. Zahlenwerte der Koeffizienten My, „ 
3 pe 4 2 | hr: eh u | 8 
r 1 | 2 3 | 4 | 5 6 | 
ak ! | | | l | 
l l | u lMae N 
ie .10° | +1,7317 10% | —6,7145-101 | +2,4901 -10° | —1,0156 -10° | +4,0838 10° | —1,6280 -10° | +6,5816 - 10° 
1 Zy2800: 10 ae 3294 10° | +1,6786 -10° | —1,2464 10° | +7,1501-10° | —3,7715 + 10° | Filsoı.. 105 — 8,0522 - 10° 
2 | +19160-10: | —1,2122-12% | +4,9748.10° | 8,6403 .10° | 5,7685 -10° | +7,2058 10% | 5,0680 100 | +3,6268 710, 
eh, | 18/6586 10° | 1,1109 10° | +7,5466 -10° | —8,9232 10° | +1,5498 10° | —2,3430 10° | —2,0910 - 10° 
4 j | 15.697110 | —_9.3714-10° | +8,5228 10° | —5,9186 -10° | -+3,3840 10° | —1,5937 - 10 
5 : | 43.7486 - 108° | —7,5227 10° | +8,4400 -10° | —7,1572 10° | +5,0416 - 107 
| | | +2,4665 - 10° | —5,8423 - 10° | -+7,7396 10° | — 7,6734 + 10’ 
zu | | | +1,6229 10° | —4,4314 -10° | +6,7436 - 10? 
5 +1,0678 - 10° | —3,3022 + 10’ 
2 | | +7,0258 « 10° 
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In Tafel 1 sind die Werte von Ms,,, angegeben. Die Abhängigkeit des größten Momentes von & 
(oder v) und A zeigt Bild 3. Hieraus ersehen wir die „entlastende“ Wirkung einer unendlichen 
Zahl gleicher kreisförmiger Löcher für die Spannungserhöhung in diesem Fall. 


0,5 04 0,3 02 0 0 


Be 
BI 
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Bild 3. Maximales Biegemoment für reine zylindrische Durchbiegung um y-Achse 
Wenn wir in (3.13) A = 0 setzen, so ergeben sich die größten Biegemomente für ein ein- 


zelnes Kreisloch: e7 
u > KR ren. 


mar __ I 2m 
ae a 34» 


Dieser Wert stimmt mit dem Ergebnis von J. N. Goodier [8] überein. 


IV. Reine zylindrische Durehbiegung um x-Achse 


In gleicher Weise wie im vorhergehenden Abschnitt können wir die Funktionen o(z) und 
v(z) für den Fall der reinen zylindrischen Durchbiegungen um x-Achse erhalten, indem wir die 
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Bild 4. Biegemomentenverteilung neben dem 
Loch für reine zylindrische Durchbiegung um 


x-Achse (v = 0,3). 


Mo 


(a) M in und M y längs x-Achse. 
(b) M,„ und M, 
(ec) M& am Lochrand. 


längs y-Achse., 


— Mix 


1: a/b = 0,25. II: a/db = 0,15. III: a/b —0. 
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Grundfunktionen g%(z) und y%z) (2.14) und die Zusatzfunktionen gU(z) und y®(z) (3.1) an- 
nehmen. Die Koeffizienten, die in den Zusatzfunktionen enthalten sind, wurden auch bis zu 
r = 8 berechnet. 

Bild 4 stellt den Biegemomentenzustand neben dem Loch für » = 0,3 und } = 0,25; 0,15 
zusammen. Der maximale Wert von Biegemoment in der Platte wird bir=a und 9 =0° 
und 180° erreicht und kann man das maximale Biegemoment in Reihe nach A und g entwickeln: 


Men Se 1 9 Te kes \ 
nn r 
M, I: Pr e+ 2 2 Mann 1 te ae Menke 


r=1 p=0 


Die Werte von M;,,p sind in Tafel 2 angegeben. Die Abhängigkeit des maximalen Momentes 
von € (oder v) und A ist in Bild 5 ersichtlich. Das Bild gibt uns zahlenmäßig die „überlastende“ 
ara einer unendlichen Zahl gleicher kreisförmiger Löcher für die Spannungserhöhung 
wieder. 

Tafel 2. Zahlenwerte der Koeffizienten Mr,p 


a 


+3,2899 - 10° | +1,7317 -10: | +6,7145 102 | +2,4901 102 | +1,0156 10° | +4,0838-10° | +1,6280 10 | -+6,5816 - 10° 
+3,2899 - 10° | +3,8964 - 10: | +1,4548 102 | +4,8597 - 10° | +2,2136 -10° | +8,4119 10° | +3,1691 10° | -+1,3323 - 105 
| +1,3160:- 10% a . ig +1,7998 - I +4,5508 «10: | +2,5455 10° | +8,3248-10° | +3,4002-10% | +1,7657 - 10° 
‚6586 - —1,4243 - +1,9131 10° | —2,2132 10° | +2,8563 - 10° | —6,9332 - 10° 5,0484 - 10° 

i +5,6971 10° | —2,9989 -10° | +2,5105 - 10° | —1,1247 - 105° | +6,6137 - 105 F:6077 -108 

| +3,7486 - 10° | — 3,3297 - 10° | +2,9664 10° | —1,8995 - 10° | -+1,2007 - 107 
| | +2,4665 - 10° | —3,0834 - 10° | +3,1396 -10° | —2,4688 - 10? 
| -+1,6229 105° | —2,6161 10° | +3,0601 - 10? 
| +1,0678 -10° | — 2,1078 - 107 

| +7,0258 - 10° 


| 
| 
| 
| 
| 


SOSTOSOUPUNDHO 


30 


iS 


N 
Oo 


—— Mmax/Mo 


Fr 
m 
Rs 
u 

FF 


0,15 020 0,25 


Bild 5. Maximales Biegemoment für reine zylindrsiche Durchbiegung um &-Achse 


V. Reine einachsige Biegung um x- bzw. y-Achse 
Wie leicht gezeigt, durch geeignetes Überlagern beider vorhergehender Lösungen für die 
reine zylindrische Durchbiegungen können wir die Lösung für die reine einachsige Biegung 


um x- und y-Achse erhalten. 

Die größten Biegemomente an den Rändern der Löcher treten bei reiner Biegung um die 
y-Achse an den Stellen © = + 90° auf, während bei reiner Biegung um x-Achse an den Stellen 
0 = 0° und 180°. 

Bild 6 zeigt das Ergebnis für reine Biegung um y-Achse; Bild 7 für um x-Achse. Für beide 
Fälle bedeutet A» 0 die unbegrenzte Platte mit einem einzelnen Kreisloch. In diesem Fall 
beträgt das größte Moment 

Ne D+3» 
Be Een 2 ee ee Del). 
M, £ 3er (5.1) 


Dieser Wert stimmt mit dem Ergebnis von H. Neuber [9] überein. 


045 020 025, 00. 


Bild 6. Maximales Biegemoment für reine Biegung um y-Ache 


Gr an m 02 os je 


015 0,20 0,25 


Bild 7. Maximales Biegemoment für reine Biegung um x-Achse 


In Bild 8 sind diese Schaulinien dargestellt. Ist das Verhältnis von Lochradius zur Stegbreite 
groß, also A— 1/2, so kann man die nähernde Spannungserhöhung auf eine andere Weise be- 
stimmen. 

G. H. Lee [10] hat mittels passender Koordinaten für ein Blech mit Hyperbelkerben 
nach Bild 9 eine Lösung gefunden. Ist der Steg im Verhältnis zum Krümmungsradius klein, 
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kann man dieses Ergebnis ohne weiters auf den vorliegenden Fall der Platte mit periodischer 
Bohrung anwenden. Der von Lee gefundene Wert für die Spannungserhöhung ist 


BE ee ee 
(3-+») arg ie! eig el 


wobei ist x = 1/24. Dies stimmt auch mit dem Werte von Neuber 
[9] überein. In Bild 8 ist das Ergebnis mit gestrichelten Kurven 
aufgezeichnet. Hieraus dünkt mich der Näherungsgrad ist befrie- 
digend für A > 0,3. 


Meinem verehrten Lehrer, Herrn Prof. Emeritus Dr. S. Higuchi 
möchte ich an dieser Stelle für sein dauerndes Interesse an dieser 
Arbeit und für viele Ratschläge herzlich danken. 


. (53), 


10 
(0) 01 [4 0,3 
Bild 8. Koeffizient der Spannungserhöhung für reine Biegung um z-Achse gegen a/b (v = 0,5; 0) Bild 9. Blech mit Hyperbelkerben 
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Die axiale Strömung von Bingham Plastiken 
in konzentrischen Rohren 
Von A. Slibar, Wien und P. R. Paslay, Schenectady, New York 


Unter Benutzung dreidimensionaler Spannungs-Verformungsgeschwindigkeitsbeziehungen wird für einen 
als Bingham Plastik bezeichneten Stoff das Strömungs- und Spannungsfeld bestimmt. Die Rechnung zeigt, 
daß ein konzentrischer Teil des Materials sich zwischen zylindrischen Grenzflächen mit gleichförmiger @e- 
schwindigkeit bewegt. Die Ergebnisse einer numerischen Auswertung der erhaltenen Lösung werden für einige 
im Tiefbohrwesen auftretende Fälle zeichnerisch dargestellt. 


The flow and velocity field for a material known as Bingham plastic is found by using 3-dimensional stress- 
deformation-velocity relations. The calculation shows that between cylindrical boundaries a comcentric 
part of the material flows with uniform velocity. A graphical representation is given of the numerical results 
of some calculations pertaining to problems such as occur in deep-well drilling. 
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En employant des relations des vitesses de temsion et de deformation & trois dimensions on determine les 
lignes de courant et les lignes de tension pour une matiere nommee substance plastique de Bingham. La cal- 


culation montre qu’une partie concentrique de la matiere se trouve ‚entre des surfaces margimales avec une vitesse 
uniforme. Les resultats d’une Evaluation numerique de la solution obtenue sont illustres pour quelgues cas 


se presentant au domaine du sondage ü grande profondeur. 


HenonsayA COOTHOLIeHUA, CBABBIBAIIMA B TPÖXMepHOM IIDOCTPaHcTBe HAIPSIKEHNE, 
nedbopmanmio u CKOPOCTB, MIA MaTepuara, Ha3bIBaeMOTO mtacTMmaccoü BUHrxaMa, ONpeNesreTch 


noJIe TeyeHnfi uU HanpsBKeHHH. BpIyucJeHuA MOKAa3bIBAlIoT, YTO YaCTb MATEPHAJIA, OTpaHu4YeHHaA - 


KOHICHTPUYeCKUMH IMIIUHNPHUYeCKHUMU OTBEPCTUAMN, HBUSKETCH C IIOCTOAHHOÜ CKOPOCTbR. 
Ina pesyABTaToB, MONYYalINmxcA IIPH UMCHEHHOM NCHONIB30BAHMN HauJeHHOTO pelIeHun, 
B HEKOTOPBIX CIIYYAHX, BA7KHBIX B CBA3U C BOIPOCAMH TIIYOUHHOTO ÖypeHun namTcH Tpadukn. 


1. Einführung 


Versuche haben gezeigt, daß verschiedene, angenähert als viskos zu bezeichnende Stoffe 
ein physikalisches Verhalten zeigen, welches durch Benutzung der Verformungsbeziehungen 
viskoser Stoffe nicht vorausgesagt werden kann. Eine technisch sehr wichtige Gruppe solcher 
Materialien (Schmierfette, Farben und Lacke, Bohrschlamm) wird in der englischen Literatur 
als Bingham Plastik bezeichnet. Im eindimensionalen Fall tritt in einer solchen Plastik Ver- 
formung erst ein, wenn eine kritische Spannung überschritten wird. Ist Fließen eingetreten, 
so ist das Verhältnis der Spannungsänderung zur Deformationsgeschwindigkeitsänderung 
konstant. 

Die wesentlichen Kennzeichen der Fettschmierung können analytisch durch Annäherung 
der Schmierstoffe als Bingham Plastiken erhalten werden). Gleichfalls ist für eindimensionale 
Strömungsvorgänge das Verhalten des bei Tiefbohrungen verwendeten, als Bohrschlamm be- 
zeichneten Materials erfolgreich unter Benutzung des eindimensionalen Binghamschen Fließ- 
kriteriums (siehe Bild 1) 


<— =0,—% für ER 0;x 


beschrieben worden’). 


Bild1. Schema der übereinander gleitenden Platten. Schubspannungs-Verformungsgeschwindigkeits-Diagramm einer Bingham Plastik 


Eine der wesentlichen Aufgaben im Tiefbohrwesen ist die Vora ür di Ö 
er ussage des für die St 
einer sekundlichen vorgeschriebenen Bohrschlammenge erforderlichen ge Obwohl 
aus Versuchen lang bekannt ist, daß Bohrschlamm als Bingham Plastik angenähert werden 


kann, ist keine brauchbare Lösung für die Ström i i ; 
i ung dies € 
ind Verrohrung bekaeyene g dieses Mediums zwischen Bohrgestänge 


Für die Strömung einer Bingham Plastik durch ei i indri 
I ıng | : ein kreiszylindrisches Rohr i 
E. Buckingham?) die Lösung für laminare Axialströmung ai en worden. Die Ba. 
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ng and Completion Operations“, 
Production, 1954, Houston, Texas. 
ling Mud“, World Oil, Februrary 1, 


tn 


2. . Math. Mech. 5 : & B a 
Ba. 37 Nr. 12 ee 1957 Slibar und Paslay, Die axiale Strömung von Bingham Plastiken usw. 443 


Arbeit gibt die Lösung für die laminare axiale Strömung zwischen konzentrischen Zylindern 
unter Benutzung von dreidimensionalen Verformungskriterien an. Die Lösung wurde für den 
Bereich praktischen Interesses im Ölbohrwesen graphisch dargestellt. 

Als Beschränkung für die Anwendung der angegebenen Lösung ist das mögliche Auftreten 
turbulenter Strömung anzusehen. Der Übergang von laminarer zu turbulenter Strömung bei 
Bingham Plastiken wurde untersucht und tritt angenähert bei einem entsprechenden Wert 
der Reynoldsschen Zahl wie beim viskosen Medium ein. Für turbulente Strömung einer 
Bingham Plastik ist üblich, das Medium als viskos anzusehen. 


2. Bezeichnungen 
Die folgenden Bezeichnungen werden in der Arbeit verwendet: 
C,D,, D, Konstante 


d;, Komponente der Verzerrungsgeschwindigkeit, z.B. in kartesischen Koordinaten 
1 on ä 
1 
Nor, Di O2 
1 | 
J= >: 'S;;S;; Quadratische Invariante des reduzierten Spannungstensors, 
4,7 
K Druckgradient in axialer Richtung der konzentrischen Zylinder, 
Ds; Druck im Querschnitt z = (0, 
0 Durchflußvolumen je Zeiteinheit, 
T,2, © raumfestes zylindrisches Koordinatensystem, 
FrsuEo innerer und äußerer Halbmesser der begrenzenden Zylinder (Bild 2), 
aD su innerer und äußerer Halbmesser.des mit konstanter Axialgeschwindigkeit strömen- 


den Materialkerns (Bild 2), 
1 
= 9,9 & 0;,„ Komponente des reduzierten Spannungstensors, auch als Spannungs- 


deviator bezeichnet, 
DV, De, V, Komponenten des Geschwindigkeitsvektors in zylindrischen Koordinaten, 


V Axialgeschwindigkeit des fließenden Stoffes zwischen den Halbmessern r* und r** 
(Bild 2), 

4012 raumfestes kartesisches Koordinatensystem, 

6; Kroneckerdelta, 

m Neigung im Schubspannungs-Verformungsgeschwindigkeitsdiagramm für reinen 


Schub (Bild la, 1b), 
Komponente des Spannungstensors, 
To kritischer Wert von 5; 


3. Die Spannungs-Verformungsgeschwindigkeitsbeziehungen der Bingham Plastik 

In einem, dem Binghamschen Fließgesetz gehorchenden Stoff tritt Fließen ein, wenn 
eine Kombination der wirkenden Spannungen einen kritischen Wert überschreitet. Liegt der 
Wert dieser charakteristischen Größe unterhalb eines kritischen Betrages, so tritt im Werk- 
stoff keine Verformung ein. Die erfolgreiche analytische Behandlung des plastischen metal- 
lischen Werkstoffes mit Hilfe des Hencky-Huber-von Misesschen Fließgesetzes legt den 
Schluß nahe, daß auch im Falle Binghamscher Materialien. die zweite Invariante des reduzierten 
Spannungstensors für die Entscheidung, ob Fließen im jsotropen Medium eintreten wird, von 
maßgebender Bedeutung sein muß. In Zylinderkoordinaten schreibt sich die zweite Invariante 
des reduzierten Spannungstensors 


1 er ’ a 
h,=5 (+ Be SH ae en, PNORL EBr 1-7 BEE ERBEREE GE © 3° 


Wie bereits von K. Hohenemser und W. Prager‘) für die Anwendung an Bingham Pla- 
stiken vorgeschlagen wurde, soll auch hier angenommen werden, daß Fließen im Material auf- 
tritt, wenn die zweite Invariante J, einen kritischen Wert überschreitet. Weiter wird angenom- 


*) K.Hohenemser und W. Prager, „Über die Ansätze der Mechanik der Kontinua““, Z. angew. 
Math. Mech. 12, 1932, 8. 216—226. . J.G. Oldroyd, „A Rational Formulation of the Equations of Plastic 
Flow for a Bingham Solid“, Proc. Cambridge Phil. Soc., 43, 1947, p. 100. W. Prager, „On ‚Slow Visco- 
Plastic Flow“, Studies in Math. and Mech. presented to R. v. Mises, Acad. Press, 1954, N.Y. 
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i er ältni hwindigkeiten 

‚ daß der Werkstoff inkompressibel und das Verhältnis der Verformungsgesc 
lei beliebigen Punkt des strömenden Mediums gleich dem Verhältnis der dort auftreten- 
den reduzierten Spannungen sei. Diesem physikalischen Verhalten entsprechend haben die 
Verfasserb5) die von K.Hohenemser und W.Prager erstmals vorgeschlagenen Ver- 


formungsgesetze 


a a 7 ae esse 
1 J—% .. 

= — I SET a 3. 

“ N VI Js: a 


angewendet. In den Gin. (2), (3) bezeichnet z, den kritischen Wert von Js weleher zur 
Erzielung einer Verformung überschritten werden muß. Die Größe u stellt die zweite Material- 
konstante dar. Die Werte von z, und u müssen durch Versuch bestimmt werden. Gl. (3) stellt die 
Verbindung zwischen der Deformationsgeschwindigkeit®) d;, und der reduzierten Spannung S;; 
unter Benutzung der Stoffzahlen 7, und # und der positiven Quadratwurzel der zweiten 
Spannungsvarianten her. 

Für den eindimensionalen Verformungsfall einer Bingham Plastik zwischen parallelen 
Platten (Bild 1) unter Einwirkung einer in x-Richtung wirkenden Kraft kann leicht gezeigt 
werden, daß sich die Beziehungen (2) und (3) auf die von Bingham angegebene Lösung speziali- 

sieren. Die Lösung ist gegeben durch: 


= Syy = Sa =: Gy. 0 
ee re en, 
und 4 
Ba Amer 2 Dar: 
1 e 
Ge = 2,09) für Gr % 


Bild 1b zeigt die Schubspannung o,, über den Abszissenwerten Qv,/oöz. Diese Lösung zeigt, 
daß die Spannungs-Verformungsgeschwindigkeitbeziehungen nichtlinear sind und daß im 


allgemeinen im strömenden Medium Bereiche vorhanden sein werden (bestimmt durch A =) 
in denen das Material sich wie ein starrer Körper bewegt. 


4. Axiale Strömung zwischen konzentrischen Zylindern 


Bild 2 zeigt die betrachtete Anordnung. Die äußere Belastung ist durch den axialen 
Druckgradienten gegeben. Die Zylinder werden als konzentrisch angenommen und zur Ver- 


Z 


Materialkern mit Axialgeschwindigkeit % ? 
A 
\ 
> 


22, u [) 
N T 


ni 


| 

“* . 

| 
ih) 


Bild 2. Anordnung des betrachteten Systems konzentrischer Rohre, 


L 


7 


») P.R. Paslay und A. Slibar, Die Fließbedingung i i 
ee nn rn ne ie Fließbedingung und das Verformungsgesetz viskoser plastischer 


*) Da der Werkstoff als inko ibel ij i i 
ee „apressibel angenommen wurde, stellt d;; auch die reduzierte Verformungs- 
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nachlässigung der Ein- und Auslaufvorgänge wird die Länge I als groß im Vergleich zum Zy- 
linderradius r, angesehen. Die für den stationären Zustand gültige Lösung muß die sechs Span- 
nungs-Deformationsbeziehungen und die Gleichgewichtsbedingungen erfüllen. 

Im zylindrischen Koordinaten lauten die Gln. (3): 


0, 1 (here 2 
or 2u VI; 3 
ver Nee le 


gsi 
= dee — = (+ =] 


1 
Ir ss (000 Fr | 


0) ei va 
eu, wel In 2 nn an oo) 


92728 VI; 5 
6), 
1% ,@_ (hr), 
a) er > 
DEE 
oz! a or i= N VI, [6; -] 
We de u 1a ] 
or r r09 7: Si 
u\ yJ } 
und die Gleichgewichtsbedingungen’) 
00,r 17.09,68..0.00.5.,, Inne 080, 
a a er 
00,0 1 0006 006, | 20re Ey 
lese ek: > alla 0 ER ITE RNRE E)8 
00,, e 00,0 005, Orr a 
or + T-00 a 02 3 £ “r 


Da die Thyxotropie des Materials nicht berücksichtigt wird, kann angenommen werden, daß 
die Stoffkennzahlen 7, und u konstant seien und weiters der Werkstoff an den begrenzenden 
starren Zylinderflächen hafte. Die Geschwindigkeit wird daher an diesen Grenzflächen ver- 
schwinden und das Spannungsfeld im Gleichgewicht mit den an den Rohrenden aufgebrach- 
ten Drücken sein. 

Aus Symmetriegründen folgt, daß die Spannungen und Geschwindigkeiten von © unab- 


hängig sind und daß 


DEN Va Ve ee re 

Die GlIn. (5) und (6) werden durch den Ansatz 
ee Re ps oe De ee ee) 
v1, De UF) a NE ee) 


erfüllt. In Gin. (8) bedeuten K den axialen Druckgradienten und p, den Druck für z= 0. Da- 
mit folgt aus (5) und (6) 


OD 0, er 
u 2 =, 7) für I ee sei) 
r2 
RE für ee) 
or 
Dumala  .K a 2: 
dr ü 


?) Die in den Gln. (5) und (6) verwendeten Normalspannungen sind ohne den durch die Schwere her- 
rührenden Teil zu verstehen. Da J, unabhängig vom hydrostatischen Spannungszustand ist, wird das Fließ- 


kriterium von diesen Beiträgen unabhängig sein. 
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Bedeutet C eine Konstante, so wird die Lösung von Gin. (12) 
C 1 | 
gg Kr D has sl RB 


i% 


ür di i ichtli itä . (10) und (11) der Wert 
Für die Bestimmung von v, muß auf Grund der Nichtlinearität der GIn ( 
von o,, zuerst bekannt sein. Die Lösungen von (10) und (11) ergeben sich unter Verwendung 


von (13) zu 


u (Br #r4@mr+D,) für ee 

a für REP n (15) 
u 4 u M 

= Ve: für lo. Sr an 


wenn die Integrationskonstanten mit D,, D, und V, bezeichnet werden. 

Mit wachsendem Druckgefälle K werden sich von der inneren und äußeren starren Zylinder- 
wand aus zwei Fließgebiete entwickeln, während das verbleibende Material sich mit ansteigender 
gleichförmiger Geschwindigkeit als starrer Körper axial verschieben wird. Abb. 2 zeigt die zur 
Bezeichnung dieser Gebiete verwendeten Symbole. r* ist der Radius der zylindrischen Grenz- 
fläche zwischen innerem Fließgebiet und unverformten Kern, r** der Radius der zylindrischen 
Grenzfläche zwischen äußerem Fließgebiet und unverformtem Kern. 

Ist der Druckgradient K positiv, dann wird im Bereich 


hr 
eine positive Schubspannung o,, vom en. als der kritische Wert z,, im Gebiet 
rt <r< r** 
eine Schubspannung vom Betrag kleiner als z,, im Gebiet 
ge 
eine negative Schubspannung vom Betrag größer als z, vorhanden sein. 
Die beiden letzten zu erfüllenden Bedingungen sind: 


1 lo, = 1, an dersStelle re r*sundvr = 
2. D, muß eine stetige Funktion von r sein. 
Bei Benutzung der ersten Bedingung führt Gl. (15) auf 
C 1 
Heinz Arte are ee 
C 1 
er zKır* ar Base ee ee 


Zur Erfüllung der zweiten Bedingung muß die Axialgeschwindigkeit v, aus den Gln. (14) und (15) 
unter Einhaltung der Randbedingungen 


Va) fürn 
N) für Eur, 


bestimmt werden. Ist der Druckgradient K genügend groß, um Fließen zu verursachen, dann ist 
das Geschwindigkeitsfeld gegeben durch: 


To 


Bü er ng! Br RT R 
£ ee Bag zu u 
D, — ° für r* < = rE* A (20) 
= To IR n C r > 
ö = u Ver) 4 u (r n)+ 5 Se für ER a © (21). 


Aus der zweiten obigen Bedingung fol 


5 t weiter bei genügend Fe, FE 
geschwindigkeit des Kernes S g a großem Wert K für die Axial 


2 T LER ar 
V SEEN: a I: 
n u ( N) 4 u (7 ers rt) 
To un: IK C x 
V = ES 0: ER 2 ee. Ren D) 'E 
0 Pr (7 WG Hi (r 1) + a Se ee: 
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Die Gln. (17), (18), (22) und (23) sind die vier Bestimmungsgleichungen für die Unbe- 
kannten C, V,, r* und r** und wir erhalten daraus 


/ Bl: Yannzaa {eh 
N a ar N a 7 1,208 || 4 
ER aa Er ARE FG u ins 
r z +14 14 
Bllskaeelt ı s® 25 
Ee| 1+ 142 3 ER AD VAL EIER IE (25), 
u ee I VL Be a a 26), 
1 2|ı+ Yır> 7 = 
CK 
Po _ m, BUAURDIE nV Er er (27) 
ok r, ei er T, To To SE Were 
% 


Für einen vorgegebenen Wert von n und einen gewählten Parameterwert C K/r3 bestimmt 


Gl. (24) den dimensionslosen Druckgradienten Kr,/ty. Mit dem so bestimmten Wert können’ die 
Gin. (25), (26) und (27) zur Bestimmung der Größen r*/r, r**/r, und u Vo/r,T, herangezogen 
werden. 


Das kleinste dimensionslose Druckgefälle EZ ‚„ welches Fließen verursacht, kann aus 
0’ krit u 
den Beziehungen (25) und (26) erhalten werden, wenn *=[r;, und r**=r, gesetzt werden. Löst 
man diese beiden Gleichungen nach den zur Einleitung des Fließens notwendigen Parameter- 
werten auf, so erhält man 
age ee N en es 
krit 


T, 
0 ie; 
0 


und 


Überschreitet sch (28) gegebenen Wert, so erhält man für das Durchfluß- 
% 
volumen je Zeiteinheit 


r, ne 


er ER AR A Er 
0-27 | 1.1 dr= 2m | en up (?—r)-+ z I 
ur r; 
"fr DENE 
+ zn (r**2 — r*?) Vo+22| 1" (r— rn) - E (#3 Dee = r:dr, 


yar 


oder dimensionslos geschrieben 


RN a\ ONE Bee 
BE else], 2 MR ap Au) a) \ SL . (80). 
nen 3\Kn a; Ib 9%, rn 6 1 AT, Rn 


Ist also größer als der kritische Wert, so bestimmen die Gin. (24), (25), (26), 
To 


i A gr ER a Be 

(27) und (30), die den Strömungsvorgang kennzeichnenden Größen een 
#Q 

una („&R . 
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er nme 


j ie (X i Verhältnisses nach 
In Abb. 3 ist das kritische Druckgefälle (2®) as Funktion des Verhä! = 


ö i i Flüssigkeit tritt bei einer 
lt. Im Gegensatz zur Strömung einer viskosen Fl it 
Ge Die Fließen a ein, wenn das dimensionslose Druckgefälle einen kritischen Wert 


überschreitet. ’ - j (en A 
i i —- — unk- 
In Abb. 4 sind für den Parameterwert 7 —= 1,5 die re = 3% 1, als = 
tion von sth gezeichnet. Das Diagramm zeigt die Abhängigkeit der radialen Dicke des sic 
T 


mit der Geschwindigkeit V, verschiebenden konzentrischen Materialkernes. 


e | 
To/knit | 


12 


0 


er 6 8 / rg 


3 R 

i i vadi i i 2rr,r9, als Funktion des dimen- 
Bild 3. Der dimensionslose kritische Druckgradient Bild 5. Das dimensionslose Durchflußvolumen uQ! ro 
(Krofto)kyig als Funktion des Radienverhältnisses r,[ro sionslosen Druckgradienten Kr,/t, für verschiedene Werte des Parameters rolt; 


Bild 4. Die Radienverhältnisse r*/r, und r**/r, der Grenzflächen Bild 6. Die dimensionslose Kerngeschwindigkeit uVo/r,T, als Funktion des 
des Materialkernes mit Axialgeschwindigkeit V, als Funktion dimensionslosen Druckgradienten Kr,/t, für verschiedene Werte des Para- 
des dimensionslosen Druckradienten Kr,/r, für rolr; 155 meters rolr; 


Für verschiedene Parameterwerte (r,/r;) zeigt Abb. 5 die Abhängigkeit des dimensionslosen 


Durchflußvolumens en vom dimensionslosen Druckgefälle .s Entsprechend der Nicht- 
* . 0 0 0 

linearität der Verformungsgesetze ist das Durchflußvolumen bis zum Erreichen des zugehörigen 
kritischen Druckgradienten gleich Null. 


Abb. 6 zeigt die dimensionslose Kerngeschwindigkeit n als Funktion von er 
T 
gewählte Werte des Parameters Ka 5 Eu ; 


T; 


für 


Bi 
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5. Bemerkungen 


Die exakte Lösung für die laminare Strömung einer Bingham Plastik zwischen konzen- 
trischen Rohren wird angegeben. Die Lösungen werden für praktisch wichtige Zahlenwerte im 
Tiefbohrwesen graphisch dargestellt und die besonderen Kennzeichen der Strömung einer Plastik 
erörtert. Die resultierende Strömung unterscheidet sich im besonderen von jener einer viskosen 
Flüssigkeit durch: 

1. zur Erzielung einer Strömung muß der Druckgradient einen kritischen Wert über- 

schreiten und 

2. ein konzentrischer Materialkern bewegt sich mit gleicher Geschwindigkeit, die radiale 

Dicke dieses konzentrischen Materialkernes wird durch ansteigenden Druckgradienten 
verringert. 

Die Dicke und die Geschwindigkeit des sich als starrer Körper bewegenden Material- 
teiles sind für die Bestimmung des durch einen vorgegebenen Druckgradienten geschaffenen 
Durchflußvolumens von Bedeutung. 

Die in den Gleichungen (24) bis (30) angegebenen Lösungen gehen für verschwindende 
kritische Schubspannung z, in die für die viskose Flüssigkeit gültigen Lösungen über. 
Eingegangen am 8. Februar 1957. 


Eingegangen am 8. Februar 1957. 


Beiträge zur Diaz-Greenberg-Methode 
Von Rüdiger Nicolovius in Hamburg !) 


Das Verfahren von Diaz und Greenberg läßt sich, wie man durch Heranziehung der Hyperkreismethode 
von Synge zeigen kann, auf allgemeinere Randwertaufgaben anwenden. Der Beweis wird für Randwertaufgaben 
vierter Ordnung geführt, das Verfahren an einem Zahlenbeispiel erläutert und mit der Methode von Maple und 
Synge verglichen. Schließlich werden einige andere Aufgaben genannt, auf die die Methode anwendbar ist. 


By means of the hypercircle method due to Synge it can be shown that the Diaz-G'reenberg method may also 
be applied to boundary value problems of a more general character. T'he proof is given for 4th order boundary value 
problems. T'he method is illustrated by a numerical example and compared with the Maple-Synge method. 
Finally, some further problems are mentioned to which the method may be applied. 


On peut appliquer la methode de Diaz et Greenberg aux problömes plus generaux de valeurs marginales, 
comme on peut demontrer en considerant la methode du hypercercle de Synge. La demonstration est faite pour 
des problemes de valeurs marginales de la quatrieme ordre, la methode est explique a l’aide d’un exemple numerique 
et compare a la methode de Maple et Synge. Enfin on indique quelques autres problemes auxquels la methode est 
applicable. 


Merox Inana u ['pmnH6epra, RaAK MO5KHO HOKA3ATb BOCNONB3YACh METONOM TULEPOKPY’KHO- 
CTU CuHre, TaK;ke IIPUMeHNHM U K ÖoJlee OÖMMM KpaeBbIM 3alayaMm. JIORAa3aTelIbcTBo HaeTca 
AAIA KpaeBbIXx 3amay yeTBepToro HOPANKA, CAM IIPHeM HIHOCTPHPYeTcA Ha UNCHIEHHOM UPHMepe 
MH CONOCTABJAIAETCH C MeTonoM Marsa u CGuure. B mocHenHeä yacTu pa60TbI ykasbIBaeTca Ha 
HEeKOTOPbIe 3amayu, K KOTOPbIM ONNCBIBAeMblÜ METON TAK3KE IIPHMEHHM. 


I. Die Diaz-Greenberg-Methode für Randwertaufgaben vierter Ordnung 
1. Problemstellung, Hilfsfunktionen 
Der Aufgabe werde ein zusammenhängendes, beschränktes, offenes Gebiet B der euklidi- 
schen x-y-Ebene mit einem stückweise glatten Rand /’ zugrundegelegt, der in eine endliche An- 
zahl von zusammenhängenden Teilen so zerlegbar sei, daß diese sich zu punktfremden Mengen 
T,, I, I, I, zusammenfassen lassen (7, + [3 + 3 + T,=T). Mit den Abkürzungen 


2 
Ade==;+ = = Pt Pyyo . — Ableitung nach der äußeren Normale auf J', 


L[p] = AAp — (a 9), — (ap, + bp; NE Se 


sei die Randwertaufgabe 
L[w] =p in B 


| 
w — s M[w] = f Ml[w] = / 
o9w auf /7 ow “rauf 
ow a 
w + Da Ml[w] = f auf I, M[w] = f air, | 
Aw + qMl[w] =g dw=g 


Pa 1) Institut für Angewandte Mathematik der Universität Hamburg. 
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j etzt: a = a(t, >0inB+T,b=b@,y) > 0 in B,p = p(&y) 
Dt ar Fe en auf se I definierte Funktionen, s > 0,12 0. Die 
Randwertaufgabe () sei eindeutig lösbar. Alle Daten sollen so beschaffen sein, en ı 
Betrachtungen, insbesondere die beiden benutzten Sätze von Green, Gültigkeit e = S ). 
Das Ziel der Diaz-Greenberg-Methode für die Randwertaufgabe (1) ist es, den : = u w 
in einem Punkte von Bin Schranken einzuschließen. Dazu werden zunächst zwei Hilfsfunktionen 
u und v benötigt, die man sich so gewählt denke, daß die Bedingungen 


v—sM[ul =/ Pe £ 
dv auf I\ 0 —_ lauf; 
7 trtdu=g ra a) 
0 at 
ne 15, lu] = auf I, A auf I, 
Au -+qM[u] =y du=g 


erfüllt sind. Die Randwertaufgabe (1) sei so beschaffen, daß die Angabe eines u und eines in B 
zweimal stetig partiell differenzierbaren v möglich ist. 


2. Die Hyperkreismethode von Synge 


In diesem Abschnitt wird das Verfahren von Synge [8] — in Analogie zu seinen Betrach- 
tungen über die dritte Randwertaufgabe der Potentialtheorie — auf die Randbedingungen von (1) 
verallgemeinert. F sei der lineare Funktionenraum, der von allen reellwertigen, in B zweimal stück- 
weise stetig partiell differenzierbaren Funktionen gebildet wird. Definiert man als Skalarprodukt 
zweier Elemente o, y von F 


99-5 | dp Ay+ age tw +berv) daduy=(W9® -»:..:.0®) 
B 


und als Norm eines Elementes von F 


TE ESS oe (4), 


so sind Assoziativ- und Distributivgesetze erfüllt und die Metrik ist mindestens positiv semidefinit. 
Sie kann definit sein, was aber nicht erforderlich ist, da die für das folgende wichtige Schwarz- 
sche Ungleichung auch im semidefiniten Fall leicht zu beweisen ist. 

Der Greensche Satz lautet in der geeigneten Form 


1 
99-5 | tApay + (Pepe + 9 yı) + bp yY} dedy 
B 


küfı 1 (fo 
-5| 9 Liyldady +, | IE Ay Miy)) ar +.) 
B 7 


Mit seiner Hilfe erhält man 


(w — v, w— u) =; | w-»Lw— u] dx dy 


B 


1 dlw — 
+z || 2: > Aw) w—)Mw—u] EIS. ereln 


Durch Subtraktion der Gleichun 


& ur gen (2) von den entsprechenden Gleichungen (1) leitet man 
leicht einige Beziehungen her: 


In B: 
Mufr: L[w— u] = uw] — Zur DUPog a (7). 
w—))—sMw—u=0, Zy iiw—no, 
wegen s>0, >60 gilt hiernach 
— w—v)MfWw—u<0, 2 Alw-— u) nk ee (8,). 


°) Diese Voraussetzung schränkt die Allgemeinheit am wenigsten ein und ist in konkreten Fällen leicht 
kontrollierbar. | 


h 
’ 
4 
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Auf T5: 
° 
(w—v) +q = D) —sMw—u=0, A|w—u)+qgMWw—u =0; 
kombiniert man diese N linear so, daß die g-Glieder herausfallen, so wird 
o(w— v 
en ee, ee Kon UBER SWERERS SOREM (8,). 
Auf T},: | 
0 (w — v) 
also j 
8 (w — 
en ee — w—-v))Mw—u=0 ......... (85). 
Auf T;: 


Allzeit) ee FE el a (8): 


Während nach (7) das Gebietsintegral in (6) verschwindet, ist der Integrand des Rand- 
integrals wegen (8,) bis (8,) nirgends positiv, so daß 


DW — WU) SU. 0 ae le) 
gilt. Diese Beziehung läßt sich auch als Hyperkugel in F deuten: 


EN 


ee (u-+ v) 


Die Minimalisierung von ||u—v|| durch: einen Ansatz u=w+q,u,, v=w,+q4D, 
(Ü=1,...,n), in welchem u, und v, die inhomogenen Bedingungen (2) und die bzgl. B ortho- 
normalen?) u, und v; paarweise die entsprechenden homogenen Bedingungen erfüllen, ist wegen 
(17) (s. u.) zweckmäßig. Man denke sich diese Minimalisierung analog zu den bekannten Ver- 
fahren von Ritz und Trefftz durchgeführt. Für Einzelheiten werde auf die Arbeit [7] des Ver- 
fassers verwiesen, in welcher diese Minimalisierung für einen Spezialfall ausführlich behandelt ist. 


3. Die Formel von Diaz und Greenberg 


In der ‚‚dritten Greenschen Identität‘ 
we [F» de dy iz Aw Mal) ar | 9 wmm)ar. Sat 
B 18 on 


soll w» den Wert von w in einem festen Punkt P aus B (nicht auf T'!) bedeuten und Y die Green- 
sche Funktion, die zu diesem Punkte P und zur Randwertaufgabe d) gehört, d.h., daß y die 
homogenen Bedingungen (1) erfüllt, ausgenommen im Punkte P, wo y wie ?lnr (? =? nl 
singulär wird. Für das folgende wird vorausgesetzt, daß eine Elementarlösung h mit L[h] = 0 
(außer in P) bekannt ist, welche die gleiche Singularität besitzt, wie y, so daß w = Y — h in ganz 
B regulär ist. Dieser Regulärteil w der Greenschen Funktion ist dann Lösung einer Randwert- 
aufgabe vom Typ (1), welche sich nach der Hyperkreismethode des vorigen Abschnitts behandeln 
läßt. Zusammen mit der Randwertaufgabe für w sollen daher hier die an die entsprechenden 
Hilfsfunktionen u und v analog (2) zu stellenden Bedingungen angegeben werden: 


L[w] = L[a] =0 in B 
w— sM[o] = v—sMüül=/=—h+sMih] 
ow je. OD LEW = oh E auf /}, 
— +t4w = een 
= dw EEK ZERO en = oh n 
w+47,—5sMlu] mas ME =i> hg, +sMiA auf I, = 
Aw + qM[w] = Au+qMf[ül=9= — Ah—qMIh] ie 
M[w] = Mjül == -—Mih] 
dw @ ov DR: Oh = auf I’; 
M[w] = Miu] =/=-MlIh] er T, 
An = Au=g=- Ah j 


3) Orthonormal im Sinne von: (u v, w—-u)=6ir für ,k=1,...,n. 
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fund g sind dabei wie angegeben aus h zu berechnen. Analog zu (10) gilt die Ungleichung 
IERTERE 
| <z lä—nl ee en Zr en (13). 


Löst man die Randbedingungen für w (1) und y (12, homogen) nach dem jeweils ersten, 
nicht mit einem Faktor s, t oder q behafteten Glied auf, setzt man sie in (11) ein und spaltet man 


dann >in h + w.auf, so erhält man mit der aus bekannten Funktionen zu berechnenden Größe 

3 ie = oh 2. 
= | ir azay + m ar— | rar + | 028 ar— (gAkar 
f B Ts Tzı Ta Tja 

(dabei bedeutet I), = I; + I},) die Formel 

2 2 24 ow = 
Wp—2Z= | Braray 2 ma) ar— | war + | ar— | v4war 
B Tzı Ta Ts 


Tıa 


0,04 


er, 


e n 0 o 
- [Tammar— |Twar + [fr ar— [zAwar . 110 
Tja Tai Tza Tja 

Die zweite Zeile entsteht aus der ersten durch zweimalige Anwendung des Greenschen 
Satzes (5) auf das Gebietsintegral und entsprechende Ausnutzung der Randbedingungen für w 
und w. 

Durch einfache vektoralgebraische Umformungen, anschließende Anwendung des Green- 
schen Satzes (5) und durch Einsetzen der Randbedingungen erhält man nach Vergleich mit (1 4 


bene. a 7)+m- a4), u») (u ,0— 0) 2(u,7) 


+ (ww) + wu) —- ,w—)—-W,w—u)=...=w-—-N..... (15) 
mit 


N=2-- (uw ,ä—5)+2(u) + [ (Mful + fMfä] + sM[u] Mia) dr 
— f @ Au. 4 Au mt Ar A 


Durch Abschätzung von (15) nach der Schwarzschen Ungleichung sowie nach (10) und (13) 
ergibt sich die Formel von Diaz und Greenberg: 
w—N|< |u—ov||- |a—ovl|| |... :....2222.(2. 
Literatur: Diaz und Greenberg [4] behandeln die Aufgabe (18); Washizu [10] zerlegt Z’ in zwei Teil I, , 
7, und gibt auf I‘, die Randbedingungen von (18) vor, auf I‘, Biegemomente und Scherkräfte. 
4. Allgemeinere Skalarproduktsdefinition 
Bei der ersten Randwertaufgabe der Plattengleichung, 
AAw=p in B 


ee ne Br 


ist auch folgende Definition des Skalarproduktes brauchbar: 
1 
( y) = Do IE (1— 9) (Or: %rz + 2 Pay Yau t Pyu Yyy) +rv4Ap Ay} dry . 2. ulgr 
B 
Aus der Darstellung 


at + +2) dr 
und den Beispielen ; 
(zy,2y) = (1 —») J dx dy, @?+y,22+y)=4(l+v)[dedy 
geht hervor, daß die (19) entsprechende Metrik genau für —1<v <+Il E semidefinit ist. 


Daher ist — 1 <» < + 1 vorauszusetzen. Anwendung des 
ar g des Gaußschen bzw. Greenschen Satzes 


Auf IR 6) ° 
(9 y) = 5 | p AAy de dy + = ee | Ya ) 9p oAy 
2 R : ») 2 IN v) Pr On +9 an 169 an Ay (07) an af (20). 


. 
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Mit der „dritten Greenschen Identität“ 


r I oe die 0Y _9Aw 
WW) = AAwdxed | —_ | oe en y— — a b 
e IE nd (+) +9 Aw 7 al 
wa [are öy, 22 ESTER Zu 
IG v) (0,2%: an er AN Ale a a (21) 


kann das Diaz-Greenberg-Verfahren dann für die Aufgabe (18) analog zu den vorhergehenden 
Abschnitten abgeleitet werden. 


Die Abschätzung der Ableitungen von w bis zur dritten Ordnung einschließlich ist hier, wie auch bei Auf- 
gabe (1) durch Differentiation der Elementarlösung nach den Koordinaten des Aufpunktes P möglich. Washizu 
[10] schätzt Biegemomente und Krümmungen einer Platte ab, indem er für die Zahl v in (19) des Poissonsche 
Verhältnis des Materials der Platte einsetzt. Die vom Verfasser in [7] behandelte Methode von Maple und Synge, 
die der Diaz-Greenberg-Methode verwandt ist, läßt sich ebenfalls mit der Metrik (19) durchführen; beim folgen- 
den Beispiel wird davon Gebrauch gemacht. 


5. Beispiel 
Das Verfahren von Diaz und Greenberg werde an der Plattenaufgabe 
Adw=1 RE A a a a Ta a pn. = 
P) 2 Pad Q 
| re et QaNer 
On kl=1, W<ı 
erläutert (wel. [7]); es sei etwa die Durchbiegung wp = w (0,0) in Schranken 7 a 
einzuschließen. Dazu verfährt man zweckmäßig in folgender Weise: Ansatz 
geeigneter inhomogener und homogener Funktionen, Berechnung sämtlicher 
Skalarprodukte, Orthonormierung der homogenen Ansatzfunktionen, Trans- 
formation der Skalarprodukte mit den inhomogenen Ansatzfunktionen mittels 
der Orthonormierungskoeffizienten, Berechnung von Näherungswert und 
Schranke. I-Y2 
Hier werde ausgegangen von den Funktionen 
1 
= + Mi = 0 
ug = 0 De lan fürn 
= 3% 
© = 1 
N - 2 2 EN .& .. m 
H=H#+Yy 35, 4 se) ae 
U, = a — yi (v, ist auf r= 1 zweimal stetig differenzierbar) (23). 
= typ —5aryt + y h=0.rhr (r? = 22 + y?) 
y=a—y = + y? — 2) (0 — 1)2 
= r—6ey? + yi %, = (x? + y? — 2)2 (02 — 1)2 22 
nr Bry + l5rryt— 2, = (a? + y? — 2)2 (22 — 1)? y2 J 


Dabei bedeutet das ““ über den w; und v; (t, j# 0), daß diese Funktionen noch nicht orthonormiert sind. 
Mit den durch die Orthonormierung entstehenden Funktionen «; und v; macht man den Ansatz 


uw=w+t uw; a=W+ Su (summiert über? =]1,..., 6) 
v=u+bv, v’=Ww+ b; Y (summiert über j = 1, 2, 3) 
zur Minimalisierung der Faktoren der Schranke ||v — v|| - |@ — || in (17). Die folgenden Ergebnisse (für drei 


Werte des Parameters » in (19)) werden mit den nach der Maple-Synge-Methode erhaltenen Werten verglichen 
und anschließend diskutiert (vgl. [7]). 


v Diaz-Greenberg Maple -Synge 


lwp — 0,03153.046| < 0,00137.536 | 4,37% I vr — 0,03171.151| < 0,00271.548 | 8,57% 
hop — 0,03157.676| < 0,00141.633 | 4,49%, | |wp — 0,03175.780| < 0,00275.555 | 8,68% 
op — 0,03159.035| < 0,00142.254 | 4,51% | wr — 0,03177.139| < 0,00276.117 | 8,69% 


Dabei geben die Prozentzahlen den Bruch Bchranke 100 an. An diese Prozentzahlen mögen folgende Be- 


merkungen geknüpft werden: Näherungswert 

1. Die Diaz - Greenberg-Werte sind etwas besser, als die entsprechenden Maple - Synge-Werte. Hin- 
gegen ist die Diaz - Greenberg-Methode mit etwas mehr Arbeit verbunden. 

2. Die Ergebnisse hängen nur wenig von» ab. Daher ist es am einfachsten, v = + 1 zu setzen, denn es haben 


allgemein beim Beispiel 
fürv= +1 | harmonische Funktionen | ds und % 
für —1< vl) Ebenen I 
fürv = —1 ı Ebenen und Rotations- | % 
|  paraboloide | 


(und nur diese) im Funktionenraum F die Länge Null. 
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für r <1 auch anders definieren, als in (23), z. B. kann man Polynome in r benutzen, die 
ee bar anschließen, oder solche, die zwecks Minimalisierung der Schran- 
ichkeiten bringen aber mehr Arbeit, während die Schranken kaum 


3. Man kann v, I ; 
beir — 1 mehr als zweimal stetig differenzier 
ken einige Parameter enthalten. Beide Mögl 


den. r 
RZ ne an der Genauigkeit der Werte ist der Arbeitsaufwand sehr groß, zumal man nur Schranken für 


win einem einzigen Punkt P aus B erhält. Jedoch gibt es heute kaum andere Methoden zur Abschätzung der 
Lösung von derart allgemeinen Klassen linearer Randwertaufgaben. 


II. Die Methode von Diaz und Greenberg für andere Randwertaufigaben 


Der Kürze halber werden in den folgenden Fällen nur die Definition des Skalarproduktes, 
die den Formeln (5) und (11) entsprechenden Greenschen Sätze und die Abschätzung angegeben; 
die Beweise lassen sich mit diesen Hilfsmitteln analog zum 2. und 3. Abschnitt führen. 


6. Selbstadjungierte elliptische Differentialgleichungen zweiter Ordnung 


B sei ein zusammenhängendes, beschränktes, offenes Gebiet des N-dimensionalen, eukli- 
dischen Raumes mit kartesischen Koordinaten &,,..., &y. /} und 7, seien zwei punktfremde, 
aber den stückweise glatten Rand I’ von B ganz ausfüllende Teilmengen endlichen Zusammen- 
hangs. Über doppelt auftretende Indizes ist zu summieren, und Indizes nach einem Komma 


DR - 
bedeuten Differentiation nach den entsprechenden Koordinaten. 33 ist die Ableitung nach der 


inneren Konormalen (Definition siehe |1], S. 23). Es wird wieder eindeutige Lösbarkeit und Gül- 
tigkeit der folgenden Greenschen Sätze vorausgesetzt. Mit dem Differentialoperator 


A; we een! 5 me 
Lie) = (4; ); + : = 5 In 
[9] ( r Pır);i 19 el, ee 
schreiben sich die Bedingungen an die abzuschätzende Lösung w, an die Hilfsfunktionen u und v, 


an den Regulärteil » = y — h der Greenschen Funktion y zu einem Punkt P aus B (Die Ele- 
mentarlösung h sei bekannt.) und an die entsprechenden Hilfsfunktionen u und D folgendermaßen: 


L[w] = Lu] ur A] — L{u] = in B 

w N iur en -—/=-—h auf /, (26) 
ow Be! ER oh = en: 
ie A a RN auf /, (2 >0) 


Der Funktionenraum F bestehe aus allen reellwertigen, in B definierten und dort stetig 
differenzierbaren Funktionen. Das Skalarprodukt werde definiert und nach dem Gaußschen 
bzw. Greenschen Satz umgeformt: 


1 1 Le 
en =5 | Au tarwan-4 [pzwias—. [Ay Mar .. 27). 
B B g: 


Dabei ist A durch A,, cos w, 2) = A, = A cos (o, 2,) gegeben (v = innere Normale). Hiermit 
ist die Gültigkeit der Hyperkugelbeziehung (vgl. (10)) leicht zu beweisen. 

Literatur: [2], [8]. 

Mit Hilfe der ‚dritten Greenschen Identität“ 

erg _ dw 0y o% 0y = 
[5 an— [15% ol + [4 wa = |5ras + [41 Zar—|a [yarT (28) 
B I P B IS I% 

erhält man nach Durchrechnung des Beweises 


£ BENENNEN u ja 
0 (2 + 2(u, u) — (u— v, uU—)) + | A Va+rz ar+ | Asvv ar)<]a—o]-jü—) 
T; 


Ti 


m — nn 


| (29) 
mit der Abkürzung 


= Oh el 
«= |hrap+ [ar ®ar— [Ayhar 3 
ame. ji 


B A 


Literatur: [3], [6] und [9] (Dort wird als allgemeinster Fall s 2 0 — vgl. (26) — behandelt.). 
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7. Eine parabolische Randwertaufgabe vierter Ordnung 
Unter fast wörtlicher Übernahme der Voraussetzungen des 1. Abschnitts mögen die Aufgabe (w), die 
Hilfsaufgabe (w = 7 — h) und die vier Hilfsfunktionen (u, v, % und ») mit den Abkürzungen 
L[9] = Przzz — (a Py)y >» a=alk,y)>0, Ay stetig, 
M[g] = Przz Nr —@PyNy >» NR %y = Komponenten der äußeren Normale auf T', 
gegeben sein durch die Bedingungen 


Deal =p, 2er = Tran 0 in B 

w —u ee =) ieh auf 2, und I,, 

M[w) =M[u] =f, Mw] =Mia) =f=—-Mh) ufl,undl,\.... (31). 
Wh = U —9 N um =9g=— han: auf 7, und I’, 

UN: Ana Wrı nn =I=— Be N; auf], und T, 


Um einen Widerspruch zu vermeiden, sei g = 0 dort, wo n, = 0 ist. 
Die Definition des Skalarprodukts und die Umformung nach dem Greenschen Satz lauten: 


Ich: ua 1 
NZ, | (PrzYzz Tr ayyy)ddy=— | p Liy] dedy + 75 n (Px Yan, — pP Mly))dT . . (32). 
B 25, I. 


‚... Pie Greensche Funktion y und die als bekannt vorauszusetzende Elementarlösung h mögen der homogenen 
Differentialgleichung (ausgenommen im Aufpunkt P) genügen und den Limesbeziehungen 


lim /MipJldal=1, lim / Mihlat =1 
0— Do eo T, 

(Dabei ist über einen Kreis vom Radiuso um P zu integrieren.). y erfülle außerdem die homogenen Rand- 
bedingungen. Die ‚‚dritte Greensche Identität‘ lautet hier: 4 en 


wp EBAR +S 92 Wenns —y Mo) rl N, —w MY))dl ......88) 
Als Ergebnis der Diaz-Greenberg-Methode erhält man: 


op — (2 + 2 (u, ©) — (u — 0,0%) + SF Mlu)+f Mia) AT—f Gun t+gün)dN)| < |u—v|- | a — 5, 


13 13 


(34) 
mit I; = I; + T,. und 
z = Shpdady + SI Mh) AP — SFRAT 4 Sghadl— Sghaadl SE. Geh. 13t(39) 
8. Randwertaufgaben der Elastostatik 


- Unter Anpassung an die bisher verwendeten Bezeichnungen lautet die Aufgabe: Es sind in einem zusammen- 
hängenden, beschränkten, offenen Gebiet B des dreidimensionalen euklidischen x, — x, — x, — Raumes mit 
stückweise glattem Rand I’ drei Funktionen w; so zu bestimmen, daß die 


Kompatibilitätsbedingungen 
1 

w=- lm; wi): 
Spannungs-Dehnungs- Beziehungen in B 

wi; = Cijmn Wmn » er A RL) 
Gleichgewichtsbedingungen 

Wjj=ki, 
Randbedingungen 

zw; bij W; —h mit Wi=Wi; N; auf I' 


Nm 


erfüllt sind. Dabei bedeuten Indizes nach einem Komma Differentiation nach der entsprechenden Koordinate 
und es ist über doppelt auftretende Indizes von 1 bis 3 zu summieren. Ferner bedeuten w; die Verschiebungen, 
wij die Dehnungskomponenten, W;; die Spannungskomponenten, W;die Randspannungen, n; die Komponenten 
der äußeren Normale, c;jmn die Koeffizienten des Hookeschen Gesetzes, die so beschaffen sein mögen, daß auch 
die 21-Konstanten-Theorie der Kristallelastizität erfaßt wird. Weiterhin seien die k; (Gravitation, Eigenspannungen) 
in B und die a;j, bij, fi auf I’ so vorgegebene Funktionen, daß die Aufgabe eindeutig lösbar ist und die folgenden 
Betrachtungen, insbesondere die Greenschen Sätze, gültig sind. Überdies seien in jedem Randpunkt die durch 
die a;j und d;j gegebenen Matrizen A und B so beschaffen, daß A’ B (' = transponierte Matrix) positiv semidefinit 
ist und daß es zwei Zahlen ! > 0, m > 0 gibt, mit denen 


1A! A m BUB=E (Einheitsmatrix) . ... .. en 20... (87) 


gilt. Letzteres ist bei vielen praktischen Problemen erfüllt. Für diese Aufgabe hat Synge [8] die Hyperkugel- 
beziehung aufgestellt; die Randbedingungen sind in dieser Form jedoch zum Einsetzen in die „dritte Greensche 
Identität“ (40), und damit für die Anwendung der Diaz-Green berg- Methode nicht geeignet. 

Um die notwendige Transformation der Randbedingungen durchführen zu können, wird zunächst der 
Rand 7 nach dem Rang der Matrizen A und B gemäß folgendem Schema so aufgeteilt, daß jeder Punkt von I’ 
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genau einer der Punktmengen I’, 1, 1» I; angehört: 


Taoderf, | T; 


Die Fälle *) sind ausgeschlossen, da wegen (37) Rang A + Rang B> 3gilt. 

Auf T, werden nun die Randbedingungen (36) von links mit A1 multipliziert. 

Auf T', multipliziert man sie entsprechend von links mit B”. } day: k s 

Auf IT’, und I‘, existiert eine orthogonale Matrix F, die die (reelle, symmetrische) Matrix A 4 auf Diagonal- 
gestalt transformiert: FF’ A'’AF=D. Mit Hilfe von (37) und der Semidefinitheit von A’ B läßt sich zeigen, 
daß die Matrix K= F’(A' + B’) nichtsingulär ist. Die Randbedingungen von (36) gehen durch Linksmulti- 
plikation mit K (bis auf Faktoren # 0) in die Randbedingungen von (38) auf 7’, und IT‘, über. 

Als Vektor (u, v, w,...) des Funktionenraumes F sei eine Menge von 6 Dehnungs- und 6 Spannungs- 
komponenten bezeichnet, welche in B stückweise stetig differenzierbar sind und die Spannungs-Dehnungs-Be- 
ziehungen erfüllen. Diese werden daher nicht mehr hingeschrieben. Die transformierte Aufgabe und die Bedingun- 
gen an die Hilfsvektoren « und v lauten: 


Kompatibilitätsbedingungen 
1 l nt 
w=zlmitwmd, = 5 wi; +), \Wij=1 2, 3) 
in B 
Gleichgewichtsbedingungen 
Wi;=Uyi=hi; 
Randbedingungen m. we Eee (38) 
w+ WW; =-ut+n;l;=fhi auf I, 
W;+ GW = U; + GG; =; auf I, 
wa u — fa 
W =D; =h 
W.+sw =U.+sy =fk aufl, 


Punkt zu Punkt ändern. Nach der oben beschriebenen Transformation gilt s> 0. Es wird zusätzlich voraus- 
gesetzt, daß die I, aus je.endlich vielen zusammenhängenden Teilen bestehen. Das Skalarprodukt werde wie 
folgt definiert und nach dem Gaußschen bzw. Greenschen Satz umgeformt (p erfülle dabei die Kompatibilitäts- 
bedingungen): 


i ie) 1 
en = 4 mw} | mar | mwi,jaB EEE NEN) 
B T B 


(dB = Volumenelement in B, dI’= Oberflächenelement auf T). Setzt man die Reziprozitätsrelation 
pij Pij = Yij Pi; und die positive Definitheit der Dehnungsenergie 5 Pi ®;; voraus, so ist das Skalarprodukt 


(39) kommutativ und die entsprechende Metrik positiv definit. Beide Forderungen an die Cijmn aus den Span- 
nungs-Dehnungs-Beziehungen sind bei der 21-Konstanten-Theorie der Kristallelastizität erfüllt. 


Die ‚‚dritte Greensche Identität‘ lautet; 


ee - „äh = n PeB 
up) = [PP Wi,dB— [PP Widl+ SwTiwar ni en. 


_Die Regulärteile %® der drei Greenschen Vektoren 7% — 9% 1 und die zugehörigen Hilfsvektoren a® 
und 9% erfüllen die & 


Kompatibilitätsbedingungen 


e 1 1 
(ed) ) 7,(k) Ar) an, 
wK) — wer. u (k) — —__ (meh, (K), As 
ij D) ( er ei) ’ Y%; D} KO ) (i, ,k= 1, 2, 3) 
Gleichgewichtsbedingungen I 
ww». _ U. — 
Win; = U; = 0 
Randbedingungen 
716) .. Weck) _ 7566) en RT 7 12; le. 
a THWr = Ton Up == PM — 0,80 auf 7, 
(k) EEE TIE, ER en H 7 
WP+g u = UM +; vn == —Hb_—_ ih® aufT, ” 
ud — och) = fm EN hd 
17(k) Ru m. — auf J' 
w — Ub = jo en Hi 3 
ck) u ETC =, m 7a 7 
WVozren —ZUw, im — = H® _— sh® auf 7‘, 


uk) We — 5% U % 1 u 
wrsWP = +3U® == —hb—sH® zur, 


ee I: 


Z. angew. Math. Mech. , . . 
Be 11/12 Nov.'Dez.1957 Ludwigund Rollnik, Theorie von Regelsystemen usw. 457 


Die Durchrechnung des Beweises führt auf: 
(k=1, 2,3) | 


A| u e ers: (42) 
mit 
N=2— mn, ab 456) + Fu fRAT— fu FRAar 
In "Ds . 
+ SW AT — SuiPar+sUVTo ar 
N’zı T5 : T; 
nee 43 
= 4 —utva® ob) SP kaB+ SUP hal — SvR far nn 
> IB 7 b T, 
SUR IE MAT— ob far + y UP rar 
en T, r, 
und 


2 = NP dB + [HR AT — rw HAT + SHE SG NRA)AT—SK®LAT+SHRLAT (44). 
24 102 IE T, T; - 


Abschätzungen für die Ableitungen der w; und damit für die Dehnungs- und Spannungskomponenten 
werden durch Differenzieren der Elementarlösungen erhalten. 


Literatur: [3], [5] (Dort werden die Fälle ;;= 0, 7=T,und f=T,+ T, behandelt, vgl. (38).) 


Herrn Prof. Dr. L. Collatz, unter dessen Aufsicht dieser Aufsatz vorbereitet wurde, danke 
ich für sein reges Interesse, Herrn Dr. J. Albrecht für seinen Rat bei der Formgebung des 
Manuskripts. 
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Eingegangen am 22. März 1957. 


Theorie von Regelsystemen mit zeitlich variabler Regelstärke 
Von G. Ludwig und H. Rollnik!) in Berlin 


Die bisherige allgemeine Theorie von Regelsystemen wird dadurch verallgemeinert, daß für die Regelstärke 
eine zeitliche Variation zugelassen wird, die das System in einem endlichen Zeitintervall in den gewünschten End- 
zustand überführt. Die Lösung dieses Problems gelingt mit Hilfe einer allgemeinen La place-Transformation. 
Unter der Voraussetzung, daß die auftretenden Differential-Operatoren Polynome beliebigen, aber end- 
lichen Grades sind, werden die strengen Lösungen in Integralform aufgeschrieben und asymptotische Entwick- 
lungen sowie Potenzreihen daraus gefolgert. Außerdem werden die Bedingungen für die Stabilität des Regel- 
systems angegeben. Eine Anwendung der Formeln soll in einer späteren Arbeit erfolgen. 


The existing general theory of servomechanism is generalized so as to cover systems where the control strength 
varies withtime in such a way thatthe system is brougth into the prescribed final state within a finite period of time. 
The problem is solved by using a general Laplace transformation. For the case that the differential operators are 
polynoms of any finite degree rigorous solutions in integral form are given from which asymplotie expan- 
sion and power series are derived. T’'he stability conditions of the system are also derived. Applications will be dealt 
with in a forthcoming paper, 

La theorie generale de systemes de röglage existant jusqu’ici est generalisce par le fait que pour le degre de 
reglage une variante temporaire est permise transmettand le systeme au cours d’um intervalle de temps fini dans 
l’etat final desire. La solution de ces problems r&ussit a V’aide d’une transformation generale d’apres La place. 
En swpposant que les op£erateurs differentiels apparaissants soient des polynömes de degr& quelcongue 
mais finis, les solutions exactes sous forme d’integrales sont donnees, et des developpements asymptotiques ainsi 
que des progressions de puissances en sont conclues. Une application des formules devra Etre donnee dans un 
article futur. 

30 
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CyıectBylomar 10 CHX IIop Teopus peryJIupoBOUHBIX CHCTEM o6oÖmaeTca, NOIycRKaf MIA 
pery.InpoBo4UHOH MOINHOCTU MH3MeHeHNE BO BPeMeHN, KOTOPOe B TeyeHuM KOHeUHOTO IIPOMe- 
3KYTKA BPeMeHH IIPHBOAUT CHCTEMy B TpeöyeMmoe OKOHYATEIBHOE COCTOAHHE. Pemenne 3ToüA 
samayı ynaerca NpM MOMoIMm oÖmero npeo6pasosauun JIarıaca. B IpenmoJIorkeHNun, TO 
BXONAINMe IMHeÄHBIe IndhepeHnmansHsIe OfepaTopbI ABIIAIOTCH MHOTOWIEHAMU 10060H HO KO- 
HeyHOoH CTEMEHN, NAIOTCA TOYHLIE PeIIeHHA B UHTETPAaIBHOH dopMe, U U3 HUX BbIBONATCH AcHMII- 
TOTHYecKuUA PA3JI05KeHHA U CTeNeHHBIe PANBI. HpoMe TOTO BEIBONATCA YCJIOBUA YCTOÜYUBOCTU 
perysmuofoyHoä cucrempI. IlpmmeHeHnue hopMmyıl ÖyNeT NaHHO B Pa0oTe, IIyOJIuKyeMoA NO3sMe. 


1. Einleitung 


Zur Behandlung von linearen Regelungen, für die allein heute eine allgemeine Theorie 
besteht, und über die eine ausführliche Lehrbuchliteratur berichtet ([1]), ist die Laplace-Trans- 
formation das geeignete mathematische Hilfsmittel. Die bekannten Stabilitätskriterien von 
Nyquist und die für die praktische Beurteilung eines Regelkreises so wichtige Transferfunktion 
können mit ihrer Hilfe leicht gewonnen werden. Der formale Grund für die Anwendbarkeit dieser 
Methode liegt natürlich darin, daß die betrachteten Regelungen durch lineare Differential- 
gleichungen mit konstanten Koeffizienten beschrieben werden können, die im Laplace-schen 
Unterraum zu rein algebraischen Aufgaben führen. Aber auch der Einfluß einer endlichen 
Verzögerungszeit, der auf eine Differential-Differenzen-Gleichung führt, kann mit ihrer Hilfe 
dikutiert werden ([2], Anhang IV). In den folgenden Überlegungen werden wir Integraltrans- 
formationen mit Laplace-Kern aber noch offenem Integrationsweg in der komplexen Ebene 
ausgiebig benutzen, um das Problem einer zeitlich sich ändernden „Regelungsstärke‘‘ zu be- 
handeln. 

Ohne uns auf eine spezielle Anwendung festzulegen, wollen wir unser Regelproblem durch 
eine Amplitude & (f) beschreiben, die im stationären Endzustand den Wert Null haben soll. 
Die Aufgabe der Regelung: besteht darin, das System von einem beliebigen Anfangswert Z zu 
diesem Gleichgewichtszustand hin zu „steuern“, auf jeden Fall 3 (f) im Laufe einer gewissen 
Zeit so nahe an Null heranzubringen, wie es die Toleranzen der speziellen Aufgabe zulassen. 

Die bisher behandelten Regelungen lassen sich durch folgende Differentialgleichung 
beschreiben: 


P(G) 30 + 0(£) 30 = n@ ren) RE 


Dabei sollen P undQ Polynome vom Grade n bzw. m in £ bezeichnen. P = 38 = nit) 


möge das Verhalten ohne Regelung beschreiben, so daß das Polynom P die Trägheiten, Reibungen 
und „inneren“ Kräfte des ‚freien‘ Systems enthält, während n(f) eventuell auftretende, von 
außen „eingeprägte“ Kräfte beschreibt. Die Koeffizienten von Q hängen von den physikalischen 
Größen ab, die die Regelungsstärke bestimmen und durch deren geeignete Wahl die Stabilität 
und das praktisch richtige Arbeiten des Regelkreises erreicht werden soll. 

Wir wollen in der vorliegenden Arbeit folgende Probleme betrachten: Zu einer bestimmten 
endlichen Zeit, für die wir £ = 0 setzen, soll der Regelvorgang abgeschlossen sein, d. h. der Fehler 
gleich Null oder wenigstens klein genug sein. Daher soll die Regelstärke um so größer sein, je mehr 
man sich dem Zeitpunkt ? = 0 nähert. Genauer setzen wir voraus, daß die zen Variation 
der Regelstärke durch folgenden Ansatz beschrieben werden kann: 


en a\ 30 
ION SORT 
IMEOELU 2 
wobei das neue Polynom Q* konstante Koeffizienten haben soll. Physikalisch bedeutet dieser 


Ansatz, daß die Regelung eigentlich mit Hilfe einer Größe e — er erfolgt, die kurz vor{=0 


außerordentlich empfindlich auf eine Abweichung der Größe j 
n : E {) von Null anspricht. Besteht 
für {= 0 ein endlicher Fehler 3(0), so wird &(!) an der Stelle i ı singulär. n 


Es ist bequem, die Größe e(t) als abhängi i ä j i 
‚ die (r ge Variable zu wählen; wir schreiben daher 
Ausgangsgleichung (mit leichter Abänderung der bisherigen Bezeichnungen) in Keeee Form 


P (7) to + 0 (4) = ni. N a 3 


Dabei sollen also die Differentialoperatoren durch Polynome gegeben sein: 


Po)= pr; NOEPE: 0 A 


!) Aus dem Institut für theoretische Physik der Freien Universität Berlin, 


rare 
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n Über die Inhomogenität n(f) werden wir soweit wie möglich keine speziellen Voraussetzungen 
machen. 

Wie schon erwähnt, ist es zuweilen sinnvoll, die im Regelkreis auftretenden Verzögerungen 
durch eine endliche Verzögerungszeit zu berücksichtigen. In solchem Falle muß (2) verall- 
gemeinert werden zu j 


’ 


P (te +0 (A) td = ni. DE len 


Da die Laplacetransformierte von e(E'—r) mit der vone(f) durch &;_,(z) = e'? 2,(z) 
gewonnen werden kann, kann man (2') formal in der Gestalt (2) schreiben, wobei Q aber kein 
Polynom mehr ist, sondern eine ganze Funktion in z. 

Um die Lösungen von (2) zu bestimmen, müssen wir zunächst die allgemeine Lösung der 
homogenen Gleichung 


P (410 +04): =0 Area ol a Ru 


_ aufsuchen. In der vorliegenden Arbeit wird die allgemeine Methode dafür angegeben (Abschnitt 2) 
und mit ihrer Hilfe Eigenschaften der Lösungen von (4) abgeleitet, die für die Beurteilung der 
Stabilität des Regelkreises wichtig sind (Abschnitt 3). Im 2. Abschnitt werden wir von P(z) 
und Q(z) nur fordern, daß sie ganze transzendente Funktionen sind, während der 3. Abschnitt 
speziell Polynome behandelt. Verallgemeinerungen der Ergebnisse dieses 3. Abschnittes werden . 
in einer späteren Arbeit publiziert werden. Im letzten Abschnitt wird der Spezialfall untersucht, 
daß P(z) mehrfache Nullstellen besitzt, der in der Theorie des 3. Abschnitts nicht enthalten ist. 

Die wichtigsten Ergebnisse sind die für große Zeiten gültigen asymptotischen Entwick- 
lungen (14), die Potenzreihen (20) und (24) und die Stabilitätskriterien (16). 


2. Allgemeine Lösungsmethoden 


Zur Lösung von (4) reicht die spezielle Laplacetransformation (wie sie etwa von Doetsch 
[3] ausführlich behandelt wird) nicht aus. Wir müssen vielmehr eine allgemeine Integraltrans- 
formation mit Laplace-schem Kern (vgl. z. B. [4]) durchführen, d. h. wir machen den Lösungs- 
ansatz 


TO LE a a 
& 


Dabei werden wir den Integrationsweg & in der komplexen z-Ebene erst später festlegen. 
Durch verschiedene Wahl von & werden wir die verschiedenen Lösungen von (4) erhalten. Damit 
der Ansatz (5) diese Gleichung erfüllt, muß gelten 


0=P "n | t el? G(z) dz + [ O(z) e? G(2) dd = — P Fi il G'(z) e® de + | O(z) e'? G(z) dz. 


& & 6 € 


d RE i 
Dabei wurde unter dem ersten Integral te’ = a6 e'? benutzt und partiell integriert unter 


der Voraussetzung, daß die Endpunkte von & keinen Beitrag liefern. Dies ist erfüllt, wenn 
1. & eine geschlossene Kurve auf der eventuell mehrblättrigen und verzweigten Riemann- 
schen Fläche von G(2) ist, oder 
2. e®G(z) an den Endpunkten von E verschwindet. 


Der zweite Fall tritt insbesondere auf, wenn & so ins Unendliche geführt wird, daß e'* G(z) 
stark genug gegen Null konvergiert. 
Unter diesen Voraussetzungen stellt (5) also sicher dann eine Lösung von (4) dar, wenn 
22) GE) = PA) C'C) 
ist, Hieraus folgt 


a nn 0 ae ar (D), 
 .. EEE 4 NER, 
wobei R(z) ein unbestimmtes Integral von Pe) ist,.so dal: 
) | 0 
Rilzy = eo (6') 


gilt. R’(z) ist also eine meromorphe Funktion und, falls P und Q Polynome sind, sogar eine 
rationale Funktion. Die Pole von R’(z) wollen wir mit}, = 1,2,...) bezeichnen. 
30* 
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Sie seien der „‚Größe‘‘ nach geordnet, d. h. es soll 
IA, | = IAu| für 1 <u 


sein. Die Hauptteile von R'(z) an den Polen A, seien: 


’ a, —__ BE u j 
h,(z) = Se Fr (2 == 1,)? + + ( ER 1,)* 
Nach dem Teilbruchsatz ist dann: | 
Ku), (h;(z) e= 9%(2)) A (6"'), 
vn e 


wobei k’(z) eine ganze Funktion ist und 9,(z) geeignet zu wählende Polynome sind, die die Konver- 
enz der Reihe erzwingen. 
; Ist R’(z) rational, so kann man natürlich die g,(z) = 0 setzen, da die I nur endlich viele 
Glieder enthält; weiterhin ist k’(z) dann ein Polynom. N £ 
Aus (6'’) folgt für R(z) 


ao 


R@ = 3 (hd) 9,9). kKo)n. 0 2 Eee 


.\ =1 
wobei die g,(z) Polynome mit den Ableitungen g,(z) sind; k(z) ist eine ganze Funktion mit der 
Ableitung k’(z); und die h(z) sind durch 


n a®) 1 ar, v 
h,(z) = a) log (z — A,) an RE: Teen a®) log (z—7,) — 4,2) (7b) 


gegeben. Damit wird dann schließlich nach (6) 


G(z) = ek) IT 90,9) _ exe) IT (Be 2, ER AA TEE 
v=1 v=1 

Die Funktion G(z) ist also eventuell mehrdeutig mit den A, (und oo) als möglichen singulären 
Stellen, und ist demnach eventuell als Funktion auf einer verzweigten Riemannschen Fläche 
zu interpretieren, deren Verzweigungspunkte nur an den Stellen A, und oo liegen können. Ob und 
von welcher Art die Stellen A, singulär sind, ist durch die a bestimmt. 

Um nach (5) Lösungen zu finden, kann der Integrationsweg eine oder mehrere der Stellen A, 
umschlingen, um sich entweder wieder zu schließen oder aber in einer der wesentlich singulären 
Stellen (A, oder co) so einzulaufen, daß e'?G(z) (für <0 !) beiAnnäherung von z and, oder © 
so stark gegen Null konvergiert, daß die durchgeführte partielle Integration erlaubt ist. 


3. Die Lösungen der homogenen Gleichungen, wenn P und Q Polynome sind 


Um die Lösungen (5) mit (7) auszuwerten, machen wir folgende Voraussetzungen: 
&) P und Q sind Polynome. 

ß) P hat nur einfache Nullstellen. 

y) a“) sind keine ganzen reellen Zahlen. 

6) Der Grad n von P sei größer als der von Q. 


‚Unter den obigen Voraussetzungen steht in (6’’) eine endliche Summe (wegen (x)), so daß 
die (2) = 0 gewählt werden können. Wegen (6) ist k'(z) = 0 und wegen (ß) sind die q,() = 0. 
Somit geht (7) über in 


n a) 
G(z) = UI (ze —4,) ee EL: 
z-£be u 
ni und (5) in 

f n a) 
e(t) a ee IE (z—A,) der. Ms 
- " ee (9) 
Man erhält n Lösungen &,, &,...,&, von (4) 


(die linear unabhängig sind, wie noch gezeigt wer- 
den ‚soll) wenn man als Integrationsweg Schleifen 
L, b(=1,3,2°,n)um>die Punkte ), nach Bild 1 
ns SE wählt. Wegen lim e?G(@)—0 für t< 0 sind die 
i Re (2)— © 
Bild 1. Schleifenweg zur Definition von e, an den Integrationsweg gestellten Forderungen 
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erfüllt. G(z) besitzt an den Stellen A, Verzweigungspunkte. Zur eindeutigen Festlegung 
der &, sollen die Schleifen €, in dem Blatt der Riemannschen Fläche beginnen, in 


. r . * a. e) 0; 
dem G(z) für große reelle z praktisch gleich Be ist. Dieses Blatt wollen wir das A- 
Blatt nennen. 


Wir wollen zunächst weiter voraussetzen, daß die Imaginärteile der A, alle verschieden sind: 
E).IUA)F It) für vEn. 
Dann numerieren wir so um, daß 
IKEA) SUR). Kür Bl Se ee) 
gilt. Die geraden Stücke der Integrationswege C, nach Bild 1 laufen auf der Geraden J{(z) = Jt(A,). 


a) Asymptotisches Verhalten für große Zeiten 


Um die e, für große Werte von (— !) zu untersuchen, d.h. für die Anfangszeit, suchen wir 
eine asymptotische Entwicklung. Wir führen als neue Integrationsvariable 


Witz 
in (9) ein. Transformieren wir noch für jedese, den Punkt (— {)A, in den Nullpunkt durch 


w=u-+l),, 
so wird 
A & j ae) 0 \u) 
= (—b) (ir) er BE TET Er +4, UN, dw. ee Sy 
6 4 et / 
wobei &, eine zu E, analoge Schleife um den Nullpunkt der komplexen w-Ebene ist. Der Strich 
am Produktzeichen bedeutet, daß der Faktor u = v fehlen soll. 


(+1) 


Wir müssen in (11) allerdings noch überlegen, welche Werte unter (— 1) und 


ae) ®% a(#) 5 

unter w IT’ nn +4, =) -! verstanden werden müssen. Wirsetzen fest, daß wegen— >00 
Te 

+1) ” le u 


rm" 


sein soll, wobei log (— ?) reell und positiv zu wählen ist! Nachdem (11a) festgelegt ist, können 


nella) 


wir den anderen Ausdruck w IT er ) nicht mehr willkürlich wählen, da die Schleife &, 

aus der Schleife €, durch die obigen Transformationen hervorgegangen ist. Geht der Realteil 

von z (bei konstanten Imaginärteil, so wie es der Weg, angibt) gegen + ©, so auch der von w. 
Za) 


N a®) = i 
Da für den Anfang von E, /T(z—4,) ” praktisch den Wert |z|” ' annimmt, muß also 
v—] 


n (v) as.v) v) n ‚(w) 
re Par en —, er ehe 
v=1 K=LN 7 : 
für den Anfang von @, sich dem Wert 


a") a®)| w | zu) 
tm 


nähern, d. h. aber wegen (11) für den Anfang von 6, gilt: 


EZ’ al“) 


w w) ae a 
1 . 


a a ul er 


Wir wollen dieses Blatt auch das A-Blatt der Riemannschen w-Fläche nennen. 
a (a) 


w : ! 
Unter dem Wert, den das Produkt Ir Er +4, ar =! annimmt, wenn man auf der reel- 


ae) 
at 
u 


len Achse dieses A-Blattes in den Nullpunkt wandert, wollen wir 
n aka) 
II(, >” Au) % 
u=1 
verstehen. Damit schreibt sich (11) in der Form 
M a(#) See, ey ee er 
& — Y (A, —4,) — (— t) ( Ei ) ua, er) 1 II 1 


„1 


w 


oe) 


\u() 
do (12), 


u=1 oA 
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wobei In e (1 + bie » 1 wird, wenn w sich auf dem Wege &, dem Punkte w = 0 
w=1 (A —A,) 

nähert. 


Entwickelt man in (12) das Produkt IT’ in eine Taylorreihe mit Restglied nach Potenzen 
von —— 
(— 2) n w y* sh ( ei 
1 a Far N et, 3 
Rt a A 
so kann man in (12) gliedweise integrieren. Weil der Limes 


) 
im (— Gr, erw R,, dw 


(N)>» j 


verschwindet, wie man leicht abschätzt, erhält man auf diese Weise eine asymptotische Ent- 


une 


rg 


wicklung nach Potenzen von > . Benutzt man noch den bekannten Wert des Integrals 


[et wdv ee (er" Fl Eee re ee (13) 
& 
und elementare Eigenschaften der Gamma-Funktion, so erhält man 


ae (a?) +1) er 


N a(#) 2= 
\ 1) Ta) + I a, — 1) = =) 


u=1 
x ı +2] 3 ce +1)(a® + 2)....- (a®) +0) ER SIT 


Da die A, alle verschieden sein sollen, folgt aus (14) die lineare Unabhängigkeit der Lösungen. 
(14) gestattet, den Verlauf der Lösung nach dem Einschalten der Regelung zu berechnen. Be- 
trachten wir als erste Näherung nur das Glied mito = 0: 


nd 


a ia) L «) - (aaa 
a en 
u=1 


Mit}, =u,+io, und a®” = o* +in®) ist der zeitliche Verlauf der Amplitude von e, 
durch 


—c () 7} 
a 
gegeben. Da die Abweichung 3 proportional (— 2) e ist, ist ihre Amplitude proportional 


(—) re 

Die Nullstellen des Polynoms P(z), die die Frequenzen der Eigenschwingungen des nicht 
geregelten Systems sind, machen sich also bei Vorhandensein einer Regelung kurz nach Beginn 
des Regelprozesses bemerkbar, wie man es auch erwarten muß. Allerdings ist die Amplitude der 
Schwingungen zeitabhängig und diese Zeitabhängigkeit kann durch geeignete Wahl der Regelung, 
d.h. des Polynoms Q(z), passend gewählt werden. Im Gegensatz zu «, gehen nämlich in or) 
die Koeffizienten des Polynoms Q(z) ein. Für genügend große —t nimmt (15’) ab, wenn gilt 


iu <oO ERBE RN: 


d. h. aber, wenn das ungeregelte System stabil ist. Ist es dagegen instabil: u, < 0, so kann die 
Regelung zunächst (d. h. kurz nach ihrem Einschalten) keine Stabilität erzwingen. Wichtig ist 
aber der Fall, daß P(z) rein imaginäre Wurzeln hat, die freie Bewegung also rein periodische 
Komponenten besitzt. Hier kann man durch die Regelung eine Stabilisierung erreichen, indem 
man Q(z) nämlich so wählt, daß 


—c®) s HM, 


a) SO 2 2 
ist. 
Einen wesentlichen stärkeren Einfluß muß das Polynom Qim weiteren Verlauf des Prozesses 
bekommen. 


b) Entwicklung der Lösungen für kleine Zeiten 


Daher wollen wir die Lösungen e für kleine Zeiten (— f) untersuchen. Leider läßt sich eine 
Potenzreihenentwicklung nach Potenzen von-(— £) nicht aus (9) gewinnen, da der Integrations- 


weg, ins Unendliche geht, so daß eine Entwicklung von e'? in eine Potenzreihe im Integral 
unmöglich ist. 
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Wir versuchen deshalb andere Lösungen von (4) zu gewinnen. Statt der in Bild 1 dargestellten 
Wege €, suchen wir ganz im Endlichen verlaufende auf der Riemannschen Fläche von G(z) 
geschlossene Wege. Wir betrachten eine Doppelschleife Z, um A, einerseits und alle anderen 
A, =# A, andererseits, wie sie in Bild 2 dargestellt ist. . Diese Schleife ist tatsächlich auf der Rie- 
mannschen Fläche von G(z) geschlossen, denn A, wie die übrigen A, werden sowohl im positiven 
wie im negativen Sinne je einmal umlaufen. Die Schleife 2, ist dadurch definiert, daß abwechselnd 
A, im positiven, dann die übrigen A, im negativen, dann, im negativen und dann dieA, im 
positiven Sinne umlaufen werden. Weiterhin soll &, so liegen, daß es durch die Halbgerade 
Jt(z) = Jt(},) und Ri(z) < R{(A,) genau zwei- 
mal geschnitten wird. Geht man auf dieser Halb- 
graden von dem Schnittpunkt (durch einen Stern 
in Bild 2 bezeichnet) des Teiles &, aus, der die 


Bild 2, Pochhammerscher Weg zur Gewinnung der Bild 3. Integrationsweg &,, der von den Wegen 2, linear 
Potenzreihenentwicklung unabhängig ist 


Halbgrade von unten nach oben durchtritt, und’ geht auf dieser Halbgraden nach rechts zu großen 
Werten von (2), so möge man sich im A-Blatt befinden. Durch diese Angaben ist &, eindeutig - 
festgelegt. 

Weiterhin führen wir noch den Integrationsweg 2, ein, der in Bild 3 dargestellt ist, und 
alleA, umläuft, ähnlich wie @, nur ein einzigesA, umläuft. Der Beginn von &, soll wieder 
im A-Blatt liegen. 

Die Lösungen (9) mit den Integrationswegen 2, wollen wir mit e;, bezeichnen und die 


mit dem Integrationsweg 2, einfach mit &,. 


Bild 4. Zur Deformation von 2, Bild 5. Deformation von 2, 


Wir können nun den Integrationsweg 2, nach Bild 4 deformieren. Die Integrale über die 
gestrichelten Stücke des Weges gehen gegen Null, wenn man diese Teile weiter und weiter nach 
rechts in derkomplexen Ebene verschiebt, weil auf ihnen der Integrand wegen des Faktors e'? 
gegen Null konvergiert (f ist negativ). j 

Der * in Bild 4 gibt wieder den Teil des Integrationsweges, der im A-Blatt liegt. Aus der 
Deformation nach Bild 4 folgt nun 
—2riz/ al) —2ri zul) —2rizale) 

u 


7 
&,=8,—8 Ey Bee € & 


Pn Kr in ee (17). 
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"Wir haben hierbei benutzt, daß & a) = ir ist. I’ bedeutet Summation über alle u + v. 
E 7 Pa 7 


Wir wollen erst einmal den Fall betrachten, daß 7 keine ganze Zahl ist. Aus (17) folgt dann 
; Li 


( ze ah BR. = c r 
z ? > 2 Iria 
ss =\l—e ? &5, —e le 1.) 1. ARTNET 


&, läßt sich also durch die e;, und e, ausdrücken. 
Den Integrationsweg für &, aus Bild3 kann man nach Bild 5 ebenfalls deformieren, so 


daß man 


2 al z a 
2ziall) In (a) ae «2)) 1 


H9=&äte "ge &g+'''+e N oe ie). 
erhält, wobei die Numerierung der A, in (10) festgelegt wurde. 
Setzt man (19) in (17) ein, so sieht man, daß sich auch die ey, durch die e, ausdrücken 


lassen. 
Die ex, lassen sich leicht in eine für alle { konvergente Potenzreihe entwickeln. Da der 


Integrationsweg &, ganz im Endlichen der komplexen Ebene verläuft, kann man im Integral et? 
in eine Reihe entwickeln und gliedweise integrieren. Es folgt 


© (— 1)? [4 a (n) 
=! [»üe- a N 3, 3: 
A) Be Bl 


In die Gleichung für &, 


führen wir als neue Integrationsvariable u = (— t) z, ein. 
Damit geht e, über in 


> Be: . ı) ne DA\uW 
ii n ae ala 
Bel here Een ir (1 a du 2 Se 
B3) u=1 4 
Für den Beginn des Weges 2, ist dann 
Mm—1 Mm—1 
= n 2 alt 
u "rn  |u| Pr und er ( - _ S 1 
#=1 u 
zusetzen. Das Produkt entwickeln wir in eine Potenzreihe nach Es die für Rs Sl 
konvergiert, d.h. un u | 
al > CH Aeneon 


für alle u. Da der Integrationsweg 2, stetig so deformiert werden kann, daß er in dem durch (21') 
festgelegten Gebiet verläuft, darf man gliedweise integrieren. Es sei 


n vw ) 7,\ a9 6) IN 
MT. 0 | 
x „a - PX — me ee 
Damit ist 
en 5 In 
& = Mt | Pn 3 ) S \ec — : “= 
o= ) lie eu dus. cn 
0=0 ” 
wobei das auftretende Integral in (13) angegeben wurde: 
5 M—l Da In—1 #> we Fr ) 
eu % "du= G Pn — ı). Pa Par a ag a 
J HU Nn— 3 
5, N 1% 1 ee +) 
Pr Pn Pn 


Damit ist schließlich 


( 2 ri ie! — | m—1 h ) 
x Pn In—ı Bu = 
I 5 ) a) 


r en De em). 
pP 


TEE 
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Für sehr kleine Werte von (— i) ist nur das erste Glied der Reihe maßgebend: 
. nm—1 


RR e She Pe) Sek pr 


Da die Abweichung 3 proportional (— !) & ist, verschwindet 3 beim Treffen (t = 0) für die 
Lösung ey, nach (19). Für e, gilt dies nur, falls 


An—ı 
n 


ist. Andernfalls wächst 3 in der Nähe des T reffpunktes sogar stark an. So scheint (24b) eine 
notwendige Forderung an die Regelung zu sein, was aber nur nes der Fall ist, wie wir 
später zeigen wollen. 


Wir haben jetzt für die meisten Zwecke ausreichende Formen der Lösungen. Mit Hilfe 


der asymptotischen Lösungen wird man die Anfangswerte beim Einschalten der Regelung 
erfüllen (wozu allerdings noch das partikuläre Integral der inhomogenen Gleichungen notwendig 


ee. 5 1or8,, 


‚ ist) und den Verlauf der Bewegung bei größeren Zeiten nach (14) und (15) diskutieren. Dann 


wird man mit Hilfe der Formeln (17) und (19) zur anderen Darstellung übergehen, um nach (20) 
und (24) die Verhältnisse bei t » 0 zu betrachten. 


e) Ausnahmefälle. 1 ist ganzzahlig 
n 
Wir wollen jetzt einige der gemachten Voraussetzungen rückgängig machen. Zunächst 
nahmen wir unter (e) an, daß die J?(A,) verschieden waren. Ist dies nicht der Fall, so muß man 
die Numerierung, die nun nicht mehr eindeutig ist, unter der Bedingung (10) festlegen. Die 
Integrationswege lege man dann auch bei gleichen Jt(A,) = Jt(A,) so, als ob Jt(A,) > Jt(A,) 
für» < u wäre. Dann bleiben alle obigen Ableitungen erhalten. 


IAn—1 


Betrachten wir jetzt den in (18) ausgeschlossenen Fall —— gleich einerganzen reellen Zahl, 


Pa 
worin der Fall q,_ı = 0 eingeschlossen ist. q,—ı ist gleich Null, wenn der Grad von Q um mehr als eins 
kleiner als der von P ist. Den Integrationsweg 2, in Bild 3 kann man in einem großen Kreis K, 
verwandeln, da beim Umlauf um alleA, der Integrand zum selben Wert zurückkehrt und sich 
daher die Integrale über die Stücke der reellen Achse in Bild 3 gegenseitig aufheben. In der 
Potenzreihenentwicklung (23) ist also 


0 für o< 
q 1 Pn 
n 0 
[e u» den N — 
RK h — 1 5 Nn— 
K ge var BB für o> I 
Re ) 2 
Pn 
zu setzen, so daß (24) lautet 
ie 
78 mt -1 
& = 2mi 3 kart a 


@-+ 
/ o! 
et 0 
= Mar |0, - (= 1) 
Py 


Leider aber lassen sich die e,nicht mehr durch die &, und &y, linear ausdrücken, sobald Yun 
eine ganze Zahl ist. Nach (17) folgt in diesem Falle Pu 


| nn 
&, =\l—e &g- 


Die es, sind also alle bis auf einen belanglosen konstanten Faktor identisch, d. h. &x,, &, enthalten 


nur eine einzige linear unabhängige Lösung. 
Nun gilt aber mit 


27 Slim 
= A et hi — ee Pn je IE IT EEE 20) 
= —1 
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die Limesbeziehung: 


2, A für Pair ganze Zahl ee Ma: 
Pr 
Aus (18) folgt: 3; 
m— ziar 
\ 1 der, Er he 1 En rag 
er di da) er | da®) 
und im Grenzfall (26a) < 
1 Ö&z, 1 22 il) Ey ee 1 (97 
> Tori dab) gen, A) de) Imi® Auen Ab (27) 


Benutzt man die Integraldarstellung für ez, und &,, SO kann man die Differentiationen ausführen 


und anschließend, genau wie im Fall: ,_ı/p, nicht ganzzahlig, Reihenentwicklungen für kleine 
Zeiten angeben. 

Allerdings sind die Rechnungen etwas langwieriger; wir geben nur das Ergebnis an. Für die 
in (27) definierten Größen ey erhält man: 


zia®) 
= Ft) + v ER, =) - log — DH) lt) — A) — Hd) . .-.- - - (28). 
Dabei ist die ee (25) für &,() einzusetzen. Die Funktion F(f) erhält man durch Ab- 
‚leiten von (20) nach den Residuen: 


aan N u 
Ft) = > en te [re leg@—A)d ı:... 82 
= 0: vl 
e=0 2, 
Die Reihenentwicklung für G(f) lautet: 
le 
ET ee DRaı Del en nn, 
e=0 n 
Dabei sind die c, in (22) definiert und die J(s) sind durch folgende Formeln gegeben: 
2nisl für SEN 
m fü . 
HK zul R we Ufer 1 (28 c) 
> Hr ıkem Kae 5] a E al 
en nlai ae En En NE, Tür se 1 


(C = Euler-Mascheronische Konstante = 0,577...). 
Schließlich gilt für H(f) die Entwicklung: 


H() = 2ni 5) sr Kerne mann» (28d), 


1 
I == o! e+1+ 


e = Max BD — 
\ Pn 


wobei die Koeffizienten d, durch 


I ı = —_ 


- eo Ee x s: % 
= ee 2a ‚dh=0. 


definiert sind. 
Ist .—ı= 0, so gilt speziell 


= 20 KH Gy HE Koaans WR 2353 
und = 

a =; er 0, 

N "ode da unn Gau Dhasa Par 


damit verhält sich e, in diesem Falle für Meine 3 i) wie 
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d.h. De, strebt dem endlichen Grenzwert 


(ER ea 1 he (290) 
zu. 
An—ı 


n 


Für negative verschwinden alle 3(t) = (— I) &(t) an der Stelle t = 0: die Regelung 


arbeitet exakt. 


Für positive en verhält sich 3%(t) wie 
n 
In—1 
ria®) I 
(e seen a te. 
Pn 
Aus (30) muß wieder (24b) folgen. Allerdings werden wir auch den Fall 
DE Le En a ee 


als eventuell möglich in Betracht ziehen müssen, nämlich dann, wenn der nach (29d) auftretende 
Fehler klein genug ist. Dies kann man aber erst entscheiden, wenn wir die Lösungen genau den 
Anfangsbedingungen anpassen. Auf dieses Problem werden wir erst in der genannten späteren 
Arbeit eingehen. 


« —1 .. ee . . . . .. 
Auch negative 1 können möglicherweise zu einer funktionierenden Regelung führen, 


n 
wenn ihr Betrag sehr klein ist. Dann wird 3,(l) erst für sehr kleine Zeiten nach Unendlich hin 
„abbiegen“, was unter Umständen keiner physikalischen Änderung der Amplitude mehr entspricht. 


4. Noch nicht behandelte Sonderfälle, insbesondere mehrfache Nullstellen von P(z) 


Bei den bisherigen Rechnungen wurde immer vorausgesetzt, daß die Riemannsche 
Fläche von G(z) wirklich verzweigt ist, daß also die Residuen von O(z)/P(z) nicht ganzzahlig 
sind. Der Fall ganzzahliger negativer Residuen kann leicht mit Hilfe des Residuensatzes erledigt 
werden, weil die Schleifenwege 6, dann in Kreise um Pole endlicher Ordnung an den Stellen}, 
entarten. Die exakte Lösung für eine solche Nullstelle von P(z) ist ein Polynom in £. Für positiv 
ganzzahlige Residuen verschwinden die entsprechenden Lösungen identisch. Man muß wieder 
zu dem Hilfsmittel des Differenzierens nach den Residuen greifen, um zu nichttrivialen Lösungen 
zu gelangen. Wir wollen die nötigen Rechnungen hier nicht aufführen, da die Methode von der 
schon im vorigen Abschnitt benutzten nicht verschieden ist, und außerdem solche singulären 
Ausnahmewerte für die physikalische Anwendung völlig ohne Belang sind. 

Dagegen kann die Möglichkeit, daß P(z) mehrfache Nullstellen hat, praktisch wichtig werden. 
Daher soll dieser Fall ausführlich behandelt werden. 

Schreiben wir (7b) in der Form 


h,(z) = a®) log (2—4,) — 4, (2—A,) 
mit 


) 
j a®) a") 1 d_a, z 
RU re Bes METER see EP TORE AL): 


so erhalten wir aus (7) und (5): 


ie ae ee... 
Pr v=1 
Hierin ist an jeder einfachen Nullstelle A, wie bisher q, = 0, an einer a,-fachen Nullstelle 
z = A, aber hat q, einen Pol von der Ordnung &, — 1. Der Integrand von (32) hat also dort eine 
wesentlich singuläre Stelle, und zwar verhält er sich in deren Umgebung wie 


K, 
es a ne 
e To) o für z>/7; 


mit der Konstanten 
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Um z = A, herum gibt es also eine solche Aufteilung des Winkels 2x in 2 (&— 1) gleiche 
Teile, daß der Integrand in (32) exponentiell gegen Null geht, falls man sich in einem nicht schraf- 
fierten Sektor in Bild 6 dem Punkte A, nähert. 

Daher lassen sich «, Schleifen wie in Abb. 6 finden, die den für erlaubte Integrationswege 
formulierten Bedingungen entsprechen, so daß wir im ganzen wieder 3 ©, —n Integrationswege 


o 
haben. Wie man zeigen kann, liefern sie lauter linear unabhängige Lösungen, und treten somit - 
an die Stelle der e,. 


Bild 6. Integrationswege im Fall einer dreifachen Nullstelle Bild 7. Integrationswege für die Potenzreihenentwicklung 
von P(z) für das asymptotische Verhalten. (Fall der Abb. 6) 


Ebenso lassen sich leicht die den ey, entsprechenden Lösungen finden, die man in Potenz- 
reihen nach f entwickeln kann. Die im Endlichen liegenden, den $%, entsprechenden Schleifen 
zeigt Bild 7 unter der Annahme, daß A, eine dreifache Nullstelle ist. 

Die Lösung &, aber ist mit ihrer Entwicklung (24) unverändert aus Abschnitt 3 zu über- 
nehmen, wenn wir nur statt (22) schreiben: 


=0 


DV zZ 


Die wesentlich neue Schwierigkeit ist die asymptotische Entwicklung der &,, die allgemein 
gar nicht übersichtlich durchführbar ist. Deshalb wollen wir hier nur einige Sonderfälle 
betrachten. Der Fall 

ODd)=-m+tN:: Pe)= Bruni SERIE EHER Bu 


läßt sich auf bekannte Funktionen zurückführen. Die zugehörige Differentialgleichung ist nämı- 
lich vom Besselschen Typ 


deren Lösung sich in der Form 


l dr 


De ze e - ) & 


angeben läßt. (vgl. [5)- (Z; ist eine Zylinderfunktion von Index A). Es lohnt sich nicht, dafür 
unseren Formalismus zu bemühen, obwohl er natürlich auch ‚‚funktioniert“. 
Dagegen bringt der Fall 


OD) = tn +p2, Pe) pP Zap, Deere 
neue Ergebnisse. Man erhält 
(N) 5, G,, Dr 
ec u En ae un (34a) 
mit 
ER u) + M/ Jo p+tHA+g7 
BR _hhrNh4 _ % ea 92 % 
= a, RR ni ee, Ir > = — a, (4b), 


E | ee (Zee "de ae 2 Pe 5, 


ie ee u pe 
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Nehmen wir zunächst @?> 0 an, so geht der Inte i 
; grand gegen Null, wenn man sich von der posi- 
ge: her dem N ullpunkt nähert. Wählen wir daher die Integrationswege wie in Bild 8, 
erhalten wir die drei Lösungen &,, &, und &,; und wählen wir die in Bild 9 dargestellten Wege 
so ergeben sich drei andere, die wir mit 5, , &x, und &, bezeichnen: 


R 
g 
” 


Bild 8. Integrationswege, wenn P(z) eine einfache und Bild 9 j 
gra ge, . Fall Bild 8: W, s 
eine zweifache Nullstelle besitzt > ve 7. 


Cı 


Für die Lösung e, läßt sich leicht die asymptotische Darstellung finden, mit z—/ = w 
lautet sie: 


- BL : © 
En er | ev (mw Iy-ımware "* dwm ei | = Din u \e du « 
2 ( SF ) (ö w e i e Er ‚ 2 b, - ; ig (36), 
& g, A) 
wobei wir den Integranden wie folgt entwickelt haben: 
(A + m)ıe Kiss DAbFur für w|<|]Al, 
= 
u 
um dann w=- I; zu setzen. 
Es folgt also für großes |!| = — ! die Entwicklung: 
f 1 a, rn: } 1 o e; hr 1 a, +1 16) 0 
&, de er 2 Are) j et nat du = (er — 1) E n e'* 3’ b,(a, + o) Ne 
= re er a er 
(37), 
und wegen b, — 4-1 e ° lautet hiervon das Anfangsglied: 
1—g ; 
2min 72 AR 
& m le" r— 1) al Amı.e He -; A a rn ey 


Um die entsprechende Entwicklung von &, und &, zu finden, setzen wir 


[ w ıw® a 
I —- Ze-—u, alsc ee Pen Sen e } 
: ER Vi ee 


und wollen zunächst untersuchen, wie dadurch die Integrationswege 6, und Q, transformiert 


werden. 
Die schlichte z-Ebene wird auf eine zweiblättrige w-Ebene abgebildet, wobei die Punkte 
z —=0 und z= beide nach w — oo übergehen. Verzweigungspunkte sind die beiden Stellen, 


wo m = —4a,tist, also 


w =2Ja,-—d) wd w=—ıw, 


daher gehen die Wege €, und &, gemäß Bild 10 in &| und &; über. 
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Verschieben wir w, durch w — w, = u in den Nullpunkt, so wird: | 
ae tfet[utw+tur(u+2m)?e.[utw+ u? (u +2 mw)? + 212]% 
& 


( Ar +a3+1 a 12 
er 


du. 
u 

u 
Entwickeln wir den Integranden nach Potenzen von mr m Es 


+++ 


‚so erhalten wir auch hier 


1 
eine asymptotische Entwicklung, deren Anfangsglied lautet: 


5 > a | . (a*du i . 
& » wı1 (—A) = e } e ) ” a ENEMIIN, 


Bild 10. Die Wege C, und C, in der w-Ebene 


Da ein verschwindend kleiner Kreis um den Nullpunkt keinen Beitrag liefert, kann man 
den Weg €, wiederum auf die doppelt durchlaufene Halbachse 0 < u <oo zusammenziehen. 
Daher wollen wir sogleich (13) benutzen und finden 


2)" wer? N | ar 
ers ca “ N 


& 


wegen \e 5) | = Yr wird 


ar v wa+uR F TER re 

Er Y2 IT (— 2)" & ar e a I 1 742 Va-2 (—) un Se 

wobei wir zum Schluß den konstanten Faktor weggelassen haben. Aus (39) erhalten wir auch- 

sofort &,, wenn wir ıw, durch w, = — w, ersetzen, diese Lösung verhält sich also für großes |[t|—= — t 
wie 

— et & 
1 FF yore 
a |— =) elaaten 00. 


Für die Lösungen &,, und &y, , deren Integrationswege ganz im Endlichen bleiben, dürfen 
wir in (35) unter dem Integral entwickeln, und erhalten gewöhnliche Potenzreihen, die bei 2— 0 
keine Singularitäten liefern können. Die Entwicklung von &, verhält sich für {— 0 wie oben 
festgestellt, ebenso wie früher, also wie 


sie liefert demnach auch hier die Bedingung (24b). 
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Über den Zusammenhang verschiedener Näherungsverfahren 


zut Berechnung nichtlinearer Schwingungen 
Von K. Magnus in Freiburg/Breisgau 


Es wird gezeigt, daß einige der in der nichtlinearen Mechanik viel gebrauchten Nährungsmethoden (Aequi- 
valente Linearisierung, Harmonische Balance, Beschreibungsfunktion, Energie-Mittelung, Fehlerquadrat- 
Minimum) als Spezialfälle des allgemeinen Ritz-Galer kin schen Variationsverfahren aufgefaßt werden können. 
Unter Verwendung der Ritz-Galerkinschen Vorschriften wird ein verallgemeinertes Verfahren zur Lineari- 
sierung durch nichtharmonische Funktionen angegeben, das für die Berechnung von Relaxationsschwingungen 
Vorteile bringen kann. 


Some approximation methods frequently used in non-linear mechanics (equivalent linearisation, harmonic 
balance, describing function, mean energy, square deviation minimum) are shown to be special cases of the general 
variational method dueto RitzandGalerkin. Following the prescriptions given by Ritz-Galer kin a generalized 
method is derived allowing the linearisation by means of anharmonic functions. It is expected that this method 
will offer some advantages in the treatment of relaxation oscillations. 


On demontre que des methodes approximatives employees sowvent dans la mecanique non-Lineaire (lineari- 
sation equivalente, equilibre harmonique, fonction de decrire, determination de la moyenne de l’Energie, minimum 
de carre d’ erreur) peuvent Etre comprises comme des cas speciaux: du procede general de variations d’apres Ritz- 
Galerkin. En appliquant les instructions de Ritz-Galer kin on indique un procede generalise pour la lineari- 
sation au moyen de fonctions non-harmoniques qui peut produire des avantages pour la caleulation des oscil- 
lations de relaxation. 


Iloka3bIBaeTcA, YUTO HEKOTOPbIE NPHÖJNLFKEHHLbIe METONbI (3KBHUBANIEHTHAA AINHeapnzanma, 
TapMOHHYyecKHÄÜ ÖaslaHc, N300Paskamıman PyHRUHA, BbIYUCHEeHHE CPenHerTo 3HAYEeHHA 3HEePTuUN, 
MHHHMYM KBanpaTa OINHÖRNM), YACTO HCHOJAIB3YEeMbIE B HeJIHHeHHOH MeXAaHHuKe, MOFKHO PaccMa- 
TPHBAaTb KaR YAcTHbIe CHAy4Yau OÖMeTO BAPHAUHOHHOTO MeTona Pnua-Tasepkuma. IIpı 1omomm 
npasna Puma-l'anepkuHa BbIBOANHTCH O000IMEHHBIH MeToN JIMHeapH3allmM C TOMOIMBIO Herap- 
MOHHYECKUX (PYHKIHÜ, KOTOPbIH MOSKET ÖbITb YMOÖHBIM MIA BBIYUCHEHNUA PeNAKCAamMOHHBIX 
Kole6aHni. 


1. Einleitung 


Für die Berechnung nichtlinearer Schwingungen wurden in den letzten 30 Jahren eine 
große Anzahl von Näherungsverfahren entwickelt, die von mehr oder weniger ähnlichen Grund- 
vorstellungen ausgehend letzten Endes zu den gleichen Ergebnissen führen. Bei der Ausarbeitung 
fast aller dieser Verfahren dürfte der Wunsch, die bequemen und weitreichenden Hilfsmittel der 
linearen Mathematik verwenden zu Können, eine ausschlaggebende Rolle gespielt haben, denn 
fast immer ist man dabei von der Methode der „Kleinen Schwingungen‘ ausgegangen, die sich 
vor allem seit Routh als außerordentlich fruchtbar erwiesen hatte. In dem Bestreben, den Be- 
reich der linearen Schwingungen zu verlassen, ohne die hier bewährten Methoden ganz aufzu- 
geben, versuchte man einerseits durch eine andere Art der Linearisierung eines vorgegebenen 
Problems, andererseits durch die Annahme von ‚‚fast harmonischen‘ Schwingungen zum Ziel zu 
gelangen. Das Schrifttum zu diesen Fragen ist derart umfangreich geworden, daß darauf ver- 
zichtet werden muß, die historische Entwicklung zu streifen. In den einschlägigen Büchern und 
dem am Schluß zusammengestellten Schrifttum findet man weitergehende Angaben. 

Für die früher fast ausschließlich behandelten schwingenden Systeme mit einem Freiheits- 
grad ist schon mehrfach darauf hingewiesen worden, daß verschiedene Näherungsverfahren zu 
identischen Ergebnissen führen. So haben bereits Krylov und Bogoljubov (1934) gezeigt, daß 
das von ihnen entwickelte Verfahren der ‚„Harmonischen Balance‘, das auch unter den Bezeich- 
nungen der „Harmonischen Linearisierung‘ oder der „Aequivalenten Linearisierung‘‘ bekannt 
geworden ist, sowie auch das Verfahren der ‚‚Energetischen Balance‘ nur Spezialfälle des allge- 
meineren Verfahrens der Variation der Konstanten sind, das in den klassischen Arbeiten von 
van der Pol zum ersten Male für die Berechnung nichtlinearer Schwingungsprobleme eingesetzt 
wurde. Lur’e [1] hat bemerkt, daß das Poincar&sche Verfahren der Störungsrechnung bei 
Beschränkung der Betrachtung auf die sogenannte erste Näherung Ergebnisse liefert, die mit 
denen der „Harmonischen Balance‘ identisch sind. Er bewies sein Ergebnis ganz allgemein für 
ein System n-ter Ordnung, in dem allerdings nur eine nichtlineare Funktion vorkommt. In sehr 
allgemeiner Form wurde der Nachweis von Bulgakov [2] geführt, der sich mit einer Verallge- 
meinerung der weitreichenden Methode der Variation der Konstanten auf Systeme n-ten Grades 
beschäftigt hat. 2 

Auf die Wesensgleichheit der Verfahren der Energiemittelung und der Harmonischen 
Balance hat Teodorcik [3] hingewiesen. Außerdem ist das in der Regelungstheorie vielfach ver- 
wendete Verfahren der ‚„Beschreibungsfunktion“ (Kochenburger [4]) nicht nur im Ergebnis, 
sondern bereits im Ansatz mit der ‚„Harmonischen Balance“ identisch. 

Neben den bisher genannten Verfahren ist auch die Variationsmethode von Ritz-Galerkin 
zur Berechnung nichtlinearer Schwingungen herangezogen worden (siehe z. B. Klotter [5], [6] so- 
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wie Lur’e und Cekmarev [7]). Daß auch diese unter gewissen Bedingungen ähnliche Ergebnisse 
wie die zuvor genannten Näherungsverfahren liefert, wurde zwar schon vermutet, scheint aber 
bisher noch nicht allgemein bewiesen zu sein. R 

Das Ziel der vorliegenden Arbeit ist nun der Nachweis, daß die genannten Näherungsver- 
fahren ganz allgemein für Systeme n-ter Ordnung als Spezialfälle der Näherungsmethode von 
Ritz-Galerkin aufgefaßt werden können. Diese Erkenntnis wird uns einerseits die Möglichkeit 
geben, gewisse Verallgemeinerungen der bisherigen Verfahren vorzunehmen, andererseits zeigt 
sie, daß keines der genannten Verfahren prinzipielle Vorteile vor den anderen besitzt. Um die 
folgenden Betrachtungen nicht zu sehr auszuweiten, sollen Einschwingungsvorgänge unberück- 
sichtigt bleiben, also nur die praktisch meist interessierenden stationären Schwingungen unter- 
sucht werden. Dabei werden wir uns auf solche Systeme beschränken müssen, bei denen Schwin- 
gungen einer bestimmten Schwingungsperiode, nicht aber Überlagerungen stationärer Schwin- 
gungen verschiedener inkommensurabler Perioden vorkommen. Diese Einschränkung dürfte 
nicht sehr einschneidend sein, da in der Praxis stationäre Zustände auch in Systemen mit be- 
liebig vielen Freiheitsgraden fast ausschließlich ‚‚ein-periodisch“ sind. 


2. Ausgangsgleichungen und Problemstellung 


Die Bewegungsgleichungen dynamischer Systeme können durch geeignete Wahl der Zu- 
standskoordinaten stets auf die allgemeine Form 
A retro) U ed ar 
gebracht werden. Bezüglich der Zeitabhängigkeit der Funktionen /, wollen wir dabei voraus- 
setzen, daß sie periodisch mit der Frequenz w seien: 


2% i 
onen 4) = han) ee N Er © > 


Wenn wir für die /, weiterhin zulassen, daß sie auch unstetig oder nicht analytisch sein 
können, dann schließen wir damit von vornherein die Anwendung der Störungsrechnung aus. 
Denn bei dieser muß für die nichtlinearen Funktionen Entwickelbarkeit in eine Potenzreihe ge- 
fordert werden. Die „Methode des kleinen Parameters‘‘ soll also hier nicht weiter beachtet wer- 
den — sie ist zudem hinreichend ausführlich von anderen Autoren (z. B. [1] und [2]) behandelt 
worden. ’ 

Für die weiteren Untersuchungen ist es nun zweckmäßig, von den Funktionen /, der rechten 
Seiten von (1) das lineare Teilsystem abzuspalten, daß sich nach den Vorschriften der Methode 
der Kleinen Schwingungen ergibt: 


N 


= Ze. m, Aigle Ba a en 
s= 
mit 
(7) 
= ar 
0%, 


wobei wir voraussetzen wollen, daß die Koeffizienten c,, die Zeit £ nicht mehr enthalten. 

Unabhängig von dem allgemeinen System (1) soll ein wichtiger Spezialfall behandelt wer- 
den, bei dem sich die rechte Seite in eine Summe nichtlinearer Funktionen von nur je einer der 
Variablen und der Zeit ! aufspalten läßt: 


= hand). N ee 


Obwohl dieser Fall bei der allgemeinen Prozedur für das System (1) miterfaßt wird, hat es sich 
als zweckmäßig herausgestellt, ihn gesondert zu behandeln, da für ihn die Verhältnisse besonders 
durchsichtig werden. 

Praktisch alle eingangs erwähnten Näherungsverfahren laufen nun darauf hinaus, die in den 
Systemen (1) bzw. (4) möglichen periodischen Lösungen (d.h. stationäre Schwingnngen) mit 
harmonischen Schwingungen zu vergleichen, die in einem hypothetischen linearen System auf- 
treten können. Wir wollen für dieses lineare Ersatzsystem die Form: 


n 
= 2 (0 Fe] KaaN 
s= ; 


wählen — und zwar gleichermaßen für die beiden Ausgangssysteme (1) bzw. (4). In der Art der 
Bestimmung der in (5) vorkommenden Koeffizienten unterscheiden sich nun die verschiedenen 
Näherungsverfahren. Als wesentlichstes Ergebnis der Arbeit wird sich zeigen, daß trotz ver- 


rt re 
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schiedener Wege in allen hier zu untersuchenden Fällen genau die gleichen Bestimmungsglei- 
chungen für die a,, und af, herauskommen. 


Nach der klassischen Methode der Kleinen Schwingungen erhielte man einfach 
Ang = Cn55 ) | 
mit den Beiwerten c,, von (3). Durch diese primitive Art der Linearisierung wird der Zugang zu 
den eigentlich nichtlinearen Erscheinungen versperrt. Der Trick aller nun zu besprechenden 
andersartigen Näherungsverfahren besteht darin, daß zwar zur Erleichterung der Rechnung in 


den Variablen linearisiert wird, daß aber die Nichtlinearität auf dem Umwege über eine Para- 
meterabhängigkeit der Koeffizienten wieder eingeschleust wird. 


3. Das Prinzip der Harmonischen Balance 
Mit dem harmonischen Ansatz 


A, kan) A sn OT nern. 


bekommt man nach Eingehen in (3) auf den rechten Seiten periodische Funktionen g,, deren 
Fourierzerlegung 


+ =bdbo+ D(b,,sinkot+bi,coskwi) 
weil 


sein möge. Subtrahiert man nun nach dem Einsetzen von (6) das ursprüngliche System (3) vom 
Ersatzsystem (5), so folgt: 


3 (a,,A,sin®, + 0% A,0c0s®,—c,,A,sin®)— bo — I (b,,5inkwt-+ b%coskof)=0. . (7). 
=1 k=1 


Um diese Gleichungen befriedigen zu können, wird zunächst b,. = 0 vorausgesetzt. Diese For- 
derung ist zum Beispiel stets für ungerade Funktionen g, erfüllt. Weiterhin werden die in der 
Fourierzerlegung von g, auftretenden Oberschwingungen vernachlässigt, so daß in (7) nur noch 
periodische Funktionen der gleichen Frequenz & vorkommen. Unter Berücksichtigung von (6) 
kann (7) dann in der Form geschrieben werden: 


N 
sin ® 12 a,, A, C0S 9, — qa7, A, © sin 9, — C,, A, C05 9,) — bl 
s=1 
n 


+ 005 © 12 (a,, A,sin 9, + 0, A,® cos 9, — c,, A, sin @,) — 2 US ERRNON: 


s=1 
Die Erfüllung dieser Beziehung ist auf viele Arten möglich, da die noch verfügbaren Un- 
bekannten a,, und a;, in verschiedener Weise gewählt werden können. Wir wollen hier (8) durch 


die Festsetzung 9, = 0 (Festlegung des Zeitnullpunktes durch den Nulldurchgang der Schwingung 
der Koordinate x,) und die Annahmen: 


b,ı 
d,ı = 61 Bi 3 Ang = [ns 
x h fürs 2.0 0 ee), 
il 
eo: al 
A,® 


erfüllen. Physikalisch bedeutet dies, daß die Grundharmonischen der Fourierzerlegung der g, 
mit den aus dem Ersatzsystem kommenden Schwingungen der Koordinate x „ausbalanciert“ 
werden. Man hätte mit dem gleichen Recht natürlich jede andere Koordinate hierzu heranziehen 
können. .. ; ” . . * . 

Setzt man in (9) die bekannten Ausdrücke für die Fourierkoeffizienten b,, und 1 ein, so 
folgt: 


2n 


) 
Te Wen len n [94 sin D, Assin.®,,.2.,D)sinotidod, | 
A) 
ir 1) De 572107 
a } 9,(A, sin D,, Assin ®,,...,1) cos wLtd(w I) | 
n A, w | 


0) 
Damit sind die Koeffizienten des Ersatzsystems (5) bestimmt. * 
Für das speziellere, aber praktisch wichtigere System (4) geht man in entsprechender W eise 
vor: durch Einsetzen des harmonischen Ansatzes (6) werden alle /,, periodisch und können in 


sl 
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Fourierreihen zerlegt werden. Dabei nimmt man die Zerlegung zweckmäßigerweise so vor, daß 
der Phasenwinkel jeder Koordinate mit berücksichtigt wird: 


Bd Hl sn DB RE) enet Pr 


Durch Vernachlässigung der höheren Harmonischen und die Annahme b,,, = 0 bekommt 
man nun wie zuvor durch Differenzbildung anstelle von (8) die Gleichungen: 


x j 
3 [sin ®, (a, A, — b,,1) + c0s Pd, (a, A,0o—bi)]=0 ..... . (2). 
s=1 

Diese Gleichungen werden durch Nullsetzen der einzelnen Faktoren befriedigt, so daß für die 

Koeffizienten die Bestimmungsgleichungen folgen: 


2n 

b 1 2 : 

a En DV | f,,(A, sin ®,, t) sin D,d®, 2 
v 


FAN ER 
Drei 
A,o 


a, 
dy,;, = 


27 
1 - N 
= a | AArsin®,. 0 00,0, 
0 


Physikalisch gesehen werden hier die Grundharmonischen der periodischen f,,-Funktionen je- 
weils mit der Schwingung ausbalanciert, die als Argument in f,, eingeht. Sowohl die Bestimmungs- 
gleichungen (10) als auch (13) sind Integraltransformationen, durch die die Funktionen g, bzw. 
f,, der Variablen x, in die von Parametern A, abhängigen Koeffizienten a,, und a7, umgewandelt 
werden. (13) besitzt vor (10) den Vorteil, daß nur eine Variable bzw. nur ein Parameter eingeht, 
so daß diese Transformation leicht auch tabellenmäßig erfaßt werden kann. 


4. Die Methode der Variation der Konstanten 


Es sei wieder das Ausgangssystem (3) vorgegeben und ein Ersatzsystem der Form (5) ge- 
sucht. Dann läßt sich der Gedankengang der Methode der Variation der Konstanten etwa folgen- 
dermaßen formulieren: 

Die Lösung des Ersatzsystems sei bekannt und — da hier nur stationäre Schwingungen 
betrachtet werden sollen — von der Form 


= Asa rA,, cos ee 


Dieser Ausdruck soll nun als Ansatz für die Lösung des nichtlinearen Ausgangssystems (3) ver- 
wendet werden, jedoch mit Amplitudenfaktoren A,, und A,,, die noch von der Zeit abhängig sind. 
Man hat also: 

= 4A,ıÖsinot -+-A,.(t)cosot, 


A 
z, = A, sinot+A,coswt+A,wcoswti— A ,wsinwt, En SEN 
—=A,sinwot+A,coswi+tit 


Geht man damit in (3) ein, so folgt unter Berücksichtigung von (5): 
A % n n 
Ausmol r A,cosat— 2 Eos %. tr Gellre Kasen a Ans Na ae, i,) . (16% 
s= sl 
Der Ansatz (15) enthält zwei noch zu bestimmende Funktionen. Da die Ausgangsgleichun- 
gen nur von erster Ordnung sind, kann noch eine willkürliche Beziehung zwischen diesen beiden 
Funktionen eingeführt werden; es ist zweckmäßig, 


A, cos ot — A,snot=0. re (17) 
vorauszusetzen. Nun lassen sich (16) und (17) als Bestimmungsgleichungen für die We und A 
auffassen und lösen: > 
x N ) 
AN 2 (Cs — 4,,) (A, sin ot +A,,cos wat, w (A, coswt — A, sin ® i)] sin w t 
z< 2 


+4....)snoL,; 
(18). 


I 


An 2 (Cs — @,,) (A, sin wi + A,, cos mt) — af, w (A,,cos wi — A,.sin @t)] cos w t 


+ Gl...) eoswt 


Der 


a SE en Sn ee 
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Die rechten Seiten sind periodische Funktionen mit der Frequenz o. Von den Amplituden- 
faktoren A,, und A,, wird aber vorausgesetzt, daß sie langsam veränderliche Funktionen der Zeit 
seien. Man wird also annehmen können, daß sie während einer Periode von der Dauer T = ee 

0) 


praktisch konstant sind. Dann aber kann man anstelle der rechten Seiten von (18) die zeitlichen 
Mittelwerte über eine Periode einsetzen. 
Zur Abkürzung der Schreibweise sei jetzt 


ae a ee 


eingeführt und dann die rechte Seite von (18) über y von O bis 2 x integriert: 


2n 
2 je 2 1 : 
A,ı > Ir 2, (C-; zu Ay) Az T +@s 1) Asa r] = 2. [9 F .) sın Y dy, 
d (20). 


EHER ef 
Ars — >, [(e,s ur 4.) As RZ Ars @ Ay ] + 3m | y .) cos y dy 
s=1 ri 
Beschränkt man sich nun auf die Untersuchung stationärer Schwingungen, dann werden die linken 
Seiten Null, und die rechten können als Bestimmungsgleichungen für die Koeffizienten a,, und 
ar, aufgefaßt werden. Es gibt auch in diesem Falle — wie schon bei der Methode der Harmoni- 
schen Balance — viele Möglichkeiten, die Gleichungen (20) zu befriedigen, da nur 2n Gleichungen 
für die 2n® Unbekannten vorhanden sind. Setzen wir jedoch, wie schon bei (9): 


Kaas ern Ai 
ei) 
dann geht (20) über in: 
2 [(&ı — %1) Aı + @ı@ Ay] = u 9(. ..)sinydy, 
J „Bald :(29): 
7 (1 — 9,1) Ale — %ı 8 Aul = — jr g,(. „.) cos y dy 
0 


Nun soll der Zeitnullpunkt — wie früher — so gewählt werden, daß für die Koordinate x, eine 
reine Sinusschwingung herauskommt, also A), = 0 und A,, = A, wird. Dann läßt sich (22) leicht 
auflösen und ergibt genau die Bestimmungsgleichungen (10). 

Für das speziellere System (4) bekommt man mit den gleichen Ansätzen (14) bzw. (15) 
anstelle der früheren Beziehung (18) jetzt: 


[ N 2 . . 
A, = X [f,(A, sin D,, tt) — a,, (A,, siny + A,sC0Sy) — %,; © (A, c0sy— A, sin y)] siny, 
2 (23). 
Äs —= 3 [f,,(A, sin D,,t) — a,, (A, siny + A,5 C0S y) — ar, @ (A,ı c0sy — A, sin y)] cosy 
s=1 


Integration der rechten Seiten über eine volle Periode und Nullsetzen der linken Seiten führt dann 
zu: 


3 (4, sin D, N) sin % dy —Ays Ası ut Ay Asa 2) re] =0, 
FE . (24). 


I 1f,,(A,sin Dt) cosydy— a, Agn — %, A,non) = 
s=1 


Durch gliedweises Nullsetzen bekommt man daraus ein System von 2n? Gleichungen für die 
zu bestimmenden 2n? Koeffizienten. Multipliziert man nun die aus der oberen Beziehung von (24) 
folgenden Gleichungen jeweils mit cos p, (bzw. sin 9,), entsprechend die aus der unteren Beziehung 
folgenden Gleichungen mit sin p, (bzw. — cos 9,) und addiert dann paarweise, so bekommt man 


unter Berücksichtigung der Gleichungen 
sin ycosp,+t cosysing, — sin ®,, sin y sin 9, — cosycosp,—=c0S ®,, 
A,ıc0spg,+A.sing,=Ä,, A, sin 9, — A, 0089, = 0 


. “ EHER » * 
wieder genau die schon vorher erhaltenen Bestimmungsgleichungen (13) für die a,, und a;,;. 
3l* 
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5. Fehlertheoretische Berechnung 


Wenn der harmonische Ansatz (6) in die rechten Seiten des Ausgangssystems (3) und des 
zugehörigen Ersatzsystems (5) eingesetzt wird, dann werden die auf den linken Seiten heraus- 
kommenden zeitlichen Ableitungen der Variablen in beiden Fällen nicht miteinander überein- 
stimmen. Der Fehler ist 


=. =46.)+ 2 TE N: 
De 
Dabei ist der aus dem Ersatzsystem folgende Wert für die Ableitung durch eine darübergesetzte 


Schlange gekennzeichnet worden. Wir verlangen nun, daß der mittlere quadratische Fehler ein 
Minimum werden soll, und errechnen den Mittelwert durch Integration über eine volle Periode: 


27 
— 1 
Pu 2 
tr EI DENE EEE 
0) 
Die Minimalbedingungen ergeben nun: 
er 2x 
. 948 1 n . 
v El Kt ER: 
04? 1 : Sr Far 
dar, = or — ag 22 ers — %)% — Ars 2 dy=0 


N) 


Geht man hier mit dem Ansatz (14) ein und wählt den Zeitpunkt wieder so, daß A,, = A, und 
A, = 0 wird, dann bekommt man durch Ausführen der Integration: - 


27 2r 
Auf LAG Q .) sin 7% dy + Ars) g,(: . :) cos y dy 
0 


n n 
= ng Ayı IT 2 [(Crs = As) Ası ar Qys © Ass] 2 Ay 41 > Ne: Fa A,;) Ay >-* Qys © Ayıl 9 
” SER E (28). 
Ayı et 9,(. . .) cos y dy— Ay f 9,(. . .) sin y dy 
0 


N n 
= — Aı®R 2 [(c,, — 4) Ag — % ® Ay) + Ag er I |, — %) Ay + GW As] 
s= s=1 : 


Multipliziert man die erste dieser Gleichungen mit A,, © (bzw. mit A,, ®), die zweite mit A 

(bzw. — A,s) und addiert dann, so folgen nach Ken durch aut . 
Akı + Ala = Ak 

die Gleichungen: 


N 1,Cer)sinydy= ee, (,— 4) Aı + %@®4,sl E 
e 29) 


[ ...)eosydy=—n > [(e,, — 0,) As — as w Ay] 


0 


Das sind aber genau dieselben Gleichungen, wie sie bei der Methode der Variati { 
s S arlation der Konstanten 
nach Nullsetzen der linken Gleichungen von (20) erhalten wurden d ie jetzi 
auf den früheren Fall zurückgeführt worden ist. a 
Im Spezialfall des Ausgangssystems (4) bekommen wir anstelle von (25) die Fehlergleichung: 


A, == .—T —e 2 Urs, l) — A,; L, — Ars &,] == S 7ER a (30) 
‚= s=1 ! 


Jetzt wird für jeden Teilfehler der quadratische Mittelwert 


- 1 
Äh = [ At, d@, 


a U A 


ET a U 
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ein System von 2 n2 Gleichungen für die a,, und a/, gewonnen. Man verwendet in diesem Fall 
am ar den Ansatz (6), in dem die Phasenwinkel mit berücksichtigt sind und bekommt dann 
aus 


— f frs(A, sin ®,,t) A, sin ®, d®, + a,4An=0, 
0 

a EREHBR): 

— f f,s(A, sin ,,t) A,w cos B,d®d, + at, A! on = 0 
ö 


Daraus folgen aber unmittelbar wieder die früheren Bestimmungsgleichungen (13). 


6. Energie-Mittelungsverfahren 


gi Bei diesem, von verschiedenen Autoren (z. B. Barkhausen [8], Den Hartog [9], Teodor- 
CiK [3]) mit großem Erfolg verwendeten Verfahren geht man von der Forderung aus, daß die 
pro Vollschwingung umgesetzte Energie für das nichtlineare Ausgangssystem denselben Wert 
haben soll wie für das Ersatzsystem. Die Energie wird durch zeitliche Integration der Leistung 
gewonnen, so daß die Ausgangsforderung — unter Fortlassung der hier nicht weiter interessieren- 
den Dimensionfaktoren — mathematisch wie folgt formuliert werden kann: 


Shut Sand. en er 
0 
oder: 
j 9.) 2,dy = j lat... Wh 


0) 0 


Mit dem Ansatz (14) läßt sich die Integration auf der rechten Seite durchführen, und man 
erhält nach gleichzeitigem Aufspalten in eine Sinus- und eine Cosinus-Komponente (was physi- 
kalisch der Trennung von Wirk- und Blind-Komponente der Leistung entspricht) das System: 


Ar] gl. smady = Aı,n 2 Kae), 0,04.) 
Be 


0 


n N ee): 
A]. 9A...) cos:y dy— At > Kar cs: aa A 
0 Se 


Das entspricht vollkommen dem früheren System (20) mit A — A,, —= (0, so daß die Rechnung 
an dieser Stelle abgebrochen werden kann. 

Für den Sonderfall eines Ausgangssystems von der Form (4) bekommt man anstelle von (34) 
die Beziehung: 


2n 


n 2% 
2 / ed — (u ta) 2 an aloh); 
s= 0 0) 

aus der durch Einsetzen des harmonischen Ansatzes, nach Ausintegrieren des rechts stehenden 
Termes sowie nach Trennung in die zwei Komponenten wiederum die früheren Gleichungen (24) 
gewonnen werden. Die eingangs geforderte Gleichheit der Energie führt also ebenfalls zu den- 
selben Bestimmungsgleichungen für die Koeffizienten des Ersatzsystems, wie sie nach den vorher 
besprochenen Verfahren erhalten werden. 


7. Die Variationsmethode nach Ritz-Galerkin 


Geht man nach den Vorschriften von Ritz und Galerkin an die Lösung der Systeme (3) 
bzw. (4) heran, so wird man nicht — wie bei Anwendung der bisher besprochenen Verfahren — 
auf ein lineares Ersatzsystem geführt. Daher soll in diesem Fall die Existenz des Ersatzsystems 
nicht von vornherein gefordert werden. Vielmehr werden wir nachträglich zeigen, daß sich die 
nach Ritz-Galerkin ausgerechneten Ergebnisse unter gewissen Bedingungen vollkommen mit 
den für das lineare Ersatzsystem (5) erhaltenen decken und daß auch jetzt wieder dieselben 
Bestimmungsgleichungen für die Koeffizienten a,, und a/, herauskommen. 


Zur Lösung von (3) gehen wir von einem Ritz-Ansatz 


%, a Pe) EEE LEN. 
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aus, bei dem die y, ein zur Approximation geeignetes Funktionensystem bilden sollen, se 
unserem Falle als periodisch mit der Frequenz » vorausgesetzt werden muß. Es liegt nahe, 
das System 


sinkwt 38) 
Nie kei, ET 
v0 et 


zu wählen. Die Bestimmungsgleichungen für die A,, folgen nun mit (3) und nach der Galerkin-. 


schen Vorschrift aus: 


t, 


[#80 00 Im a=0 an 


= 1, 2,008 
0 


Die Integration ist dabei zwischen zwei festen Zeitpunkten Z, und Z, vorzunehmen, die um den 
2 : n u: 
Betrag einer vollen Schwingungszeit T = = auseinanderliegen; wir wählen , = O0undot, =2nr 


und integrieren dann über@f=y von Obis2 x. 


Die nach der Ritz-Galerkinschen Methode zu erwartenden Ergebnisse können natürlich 
nur dann mit den Ergebnissen der anderen Näherungsverfahren übereinstimmen, wenn wir die 
Approximationsfunktionen (38) wählen und uns dabei auf das erste Funktionenpaar (k = 1) 
beschränken. Dann aber bekommt man mit dem Ansatz 


= A,,8i0y 4 A, COS y 0. or De 
aus der Galerkinschen Vorschrift (39) die Beziehungen: 
27 
| E (A,, cos y — A,, sin y) — 3 c,, (A, siny + Ay, cos y) — 1.» ) sinydy=0, 
| 
0) 
2x (41). 
| E (A,ı 608 y — Ay, sin y) — 3 C,, (A,, siny + A,,cosy) —4,(.. ) cosydy = 0 
i s=H# 
{0 


Die Integration läßt sich teilweise ausführen, so daß aus (41) folgt: 


Ann — Ir Aanr— | Gl. )sinydy—=0; 
RB 0 


S . (42). 
o Au n— 3 Ag — [ 9... cosydy—=0 
: s=1 0 
Durch Einführen des Kroneckersymbols ö,, läßt sich dies umschreiben in: 
N 1 ri e 
2 (ö,, © At As) + En [m ..) sin Ve dy —=0, 
0% ir re AR 


n 1 Y 

2 Are Ar =; 91...) C0sYy dy = 0 

v 
Aus diesem System von 2n Gleichungen können die ? n Amplitudenfaktoren A,, und A,, des 
Ansatzes (40) berechnet werden. 


Um nun zu zeigen, daß die aus (43) zu gewinnende Lösung mit den früheren Ergebnissen 
übereinstimmt, setzen wir zunächst den Ansatz (40) in das lineare Ersatzsystem (5) ein und 
erhalten (wiederum unter Verwendung des Kronecker-ö): 


2 {sin y 0,0 Aa — a, A, + ao Ay] +e0sy[6,,@ Au—a,, As — ah o® Al} =0 dA). 


Durch Nullsetzen der Ausdrücke in den 


eckigen Klammern bekommt man die Bestimmungs- 
gleichungen für die Amplitudenfaktoren: 


Ant Asa — a0, Anal, Are (45). 


Dieses System stimmt aber vollkommen mit (43) überein, sofern dieselb 


en Bezeichnungen 
und der gleiche Zeitnullpunkt eingeführt wird. 


Betrachtet man nämlich die Bestimmungs- 


En Zn nn 


2. angew. Math. Mech. = N 
Bd. 37 Nr. 11/12 Nov./Dez.199z Magnus, Näherungsverfahren zur Berechnung nichtlinearer Schwingungen 479 


gleichungen (10) und (21) als Definitionsgleichungen für neu einzuführende Abkürzungen und 
wählt den Zeitnullpunkt — wie früher — so, daß A, = A,und A, = O wird, dann kann man (43) 
in die Form bringen: 


n m 
2 (Or; @ As 5 Ars A, Zi Qfs @ As) —=0 9 3 1 © Az Fr Ups A; = ds 02) A,,) —=0 (46). 


Durch gliedweises N ullsetzen folgen aber daraus wieder die Bestimmungsgleichungen (45). 
Geht man von dem speziellen System (4) aus, so bekommt man anstelle von (42) das System: 


5 n 2x h f n 2x 
— 0A 2 SIrld, sin DD sinydy=0, @Ann— If IndArsin DD cosydy—0 (47). 
‚= ; erg 


en a früheren Abkürzungen a,, und af, von (13) einführen zu können, wird nun unter dem 
ntegra 3 


siny = sin (d, — 9,) = sin ®, cos 9, — cos ©, sin p,, 
cos y = cos (®d, — 9,) = cos B,cosy, + sin BD, sin , 
eingesetzt. Damit folgt aus (47): 


— a A,n —Y je ps frs(A,sin®,,!)sind, dd, — sin p, I1.&, sin®,t)cosd, «0 =(, 
0 0 


s=1 


(48), 
n 2x 2r 
w An — > (eo of f,,(A, sin ®,, ti) sin d, dd, + sin DS frs(A, sin D,t)cosd, a0 —N 
s= ö ö 
das läßt sich unter Berücksichtigung von (13) umformen in: 
3-9, ® Aa — 4,0059, Qq,r +A,sno, or] =0, 
s=1 : 
. (49). 


n ie 
2 [d,; ® A ı a — A, c0s 9, own — A, sing, a,,] = 0 
s= 


Berücksichtigt man nun noch, daß 
eos = An, und NArsıno7=A,, 


ist, dann geht (49) nach gliedweisem Nullsetzen genau in das System (45) über, das aus dem linearen 
Ersatzsystem erhalten war. 

Zusammenfassend kann also festgestellt werden, daß die nach Ritz-Galerkin erhaltenen 
Ergebnisse identisch mit den Ergebnissen sind, die nach den zuvor betrachteten Verfahren 
erhalten werden — sofern als Approximationsfunktionen Sinus und Cosinus gewählt werden, 
und nur das erste Glied des Ritz-Ansatzes mitgenommen wird. 


8. Eine Verallgemeinerung der Harmonischen Linearisierung 


Nachträglich läßt sich nun leicht ein Verfahren angeben, mit dessen Hilfe das lineare 
Ersatzsystem aus den Galerkin-Bedingungen unmittelbar gewonnen werden kann. Dazu 
braucht nur in den Gleichungen (39) für x, der Ausdruck (5) eingesetzt und dann die Integration 
durchgeführt zu werden. Man übersieht dabei sofort die Gleichwertigkeit der Galerkin-Bedingung 
mit den Forderungen, die den anderen Verfahren zu Grunde gelegt wurden. Aber das Ritz- 
Galerkinsche Verfahren weist zugleich den Weg zu einer Verallgemeinerung. Es können ja 
als approximierende Funktionen y;(f) beliebige periodische Funktionen verwendet werden. 
Mit jeder periodischen Funktion kann demnach eine Linearisierung des nichtlinearen Ausgangs- 
systems vorgenommen werden, die der im Abschnitt 3 besprochenen Harmonischen Linearisierung 
vollkommen entspricht. Wenngleich im allgemeinen die Harmonische Linearisierung für die 
Praxis am wichtigsten ist, so können doch Fälle vorkommen, in denen man zweckmäßigerweise 
eine Linearisierung durch andere Funktionen vornimmt. Zwei Beispiele dieser Art sollen im 
Abschnitt 9 besprochen werden. 

Um die Betrachtungen nicht durch zu große Allgemeinheit undurchsichtig werden zu 
lassen, wollen wir uns hier auf ein Ausgangssystem von der Form (4) beschränken und weiterhin 
voraussetzen, daß die darin vorkommenden Funktionen f,, nicht vom Hysteresetyp sein sollen. 
Dann werden die a*, = 0 und das Ersatzsystem bekommt die einfache Form: 


v2 
= 24,2. 2 ee er 
_= 
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Mit dem eingliedrigen Ritz-Ansatz | 
; ; 2, = A, Hl) ROT IRRUL ERS EEE (51), 


in dem zur Vereinfachung y = o t als Veränderliche eingesetzt wurde (y hat also die Periode 2), 
bekommt man die Galerkin-Bedingung: . 


2x 


[zo |rna=0- ie re 


oder wegen (50): 


27 


S [a,, % — Irs(&s: D] y(y) dy 0: 


0 
Mit der Abkürzung 


ar 
h=f yy) dy. N a es, 
6 
folgt nach Vertauschen von Summation und Integration: 
| 2x 
N 
> a Ah | [Av Z | @r =0 2.222.000. N. 
s=1 5 
Wählt man nun die Koeffizienten a,, nach 
a m, A| von 
TS h A, rs Ss 7) 
0 
an | Fear so ist (54) sogar gliedweise erfüllt, also der 


| Galerkin-Bedingung Genüge getan. Ent- 
sprechende Bestimmungsgleichungen lassen 
sich auch für die allgemeineren Ausgangs- 
systeme (3) und (4) und das allgemeinere Er- 
satzsystem (5) ableiten, nur werden die Berech- 
nungen dann umfangreicher als es bei der Har- 
monischen Linearisierung der Fall ist. 


Die Bestimmungsgleichung (55) geht mit 
y(y) = siny in die früher gewonnene ‚har- 
a monische‘‘ Formel (13) über; es wird dann 
h=n. Der Wert von h für andere y(y) läßt 
sich aus (53) im allgemeinen leicht bestim- 
men. Mit Rücksicht auf die späteren Bei- 
spiele soll hier nur eine trapezförmige Schwingung nach Bild 1 erwähnt werden, die durch 


Bild1. Die Trapezschwingung yyı(y) und ihre erste zeitliche 
Ableitung d,p(y) 


done 0<y<y 
Yo 
1 für YSYSR—YH 
N h 
Ynly)= > FÜR... sg VS Pe SYS Ehe 
0 
— 1 für zn S<Sys2n— 
MTTET mie 
für 27 —y, <y<2n 
Yo 
definiert ist. Man findet nach einfacher Integration aus (53): 
Ya) 8 y0 
hF an. A me 


Die Trapezschwingung enthält als Sonderfälle: 


die Rechteckschwingung: „y=0; h=! . (98) 


ER HNN 


und die Dreieckschwingung: y, = 5 sd 
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u Anwendungsbeispiele 


Zur Demonstration der Verwendbarkeit der verallgemeinerten Linearisierungsformel (55) 
sollen zwei einfache und besonders durchsichtige Beispiele berechnet werden, deren strenge 
Lösung bekannt ist, bzw. leicht ausgerechnet werden kann Auf diese Weise kann der Fehler 
der Näherungsberechnung unmittelbar bestimmt und die Güte der Näherung beurteilt werden. 


a) Eigenschwingungen eines selbsterregten Schwingers von einem Freiheitsgrad 
Der Schwinger genüge der Differentiagleichung: 
BEL O Der. kKsenz 00 etz 4 SERER (00) 


Solange der dimensionslose Dämpfungsfaktor D< 1 ist, kann der Schwinger selbsterregte 
Eigenschwingungen ausführen, deren Amplituden streng berechnet werden können (siehe [10)). 
Man bekommt: 

nD 


AT=K EN U. N in 
2yı — D? 
Die ebenfalls in der Arbeit [10] ermittelte harmonische Näherung für die Amplitude ist: 
2k i 
a (62). 
Ay ıD ( ) 


Für die nun auszurechnende nichtharmonische Näherung soll als approximierende Funktion 
die Trapezfunktion (56) (bzw. Bild 1) gewählt werden. Wir setzen also an: 


An, 1, . (63) 
und berechnen für die in (60) vorkommende nichtlineare Funktion einen Ersatzausdruck 
ksenz=qa& EN) 0 ee 
Nach (55) gilt für den Koeffizienten: 
EN ER | 
Be LESE Zn 
0 


Der Verlauf der zeitlichen Ableitung y7(y) ist in Bild 1 gestrichelt eingezeichnet worden. Man 
findet nun leicht: 


ar a A 
h= | Yry)y—=4 dp — „nupgaiguibait nah dam 
\ Yo Yo 
1) BR 0) 
Ak [ 4k ky 
= ET 
R Farce RA A un 


Setzt man diesen Wert in (64) ein, so folgt aus (60) das lineare Ersatzsystem: 
+ (en -)ite=0- N ee (68). 


Dauerschwingungen in diesem System haben die Frequenz © = 1 und müssen der Bedingung 
genügen: 


En Me ones. 


Daraus folgt die Amplitude der Schwingung zu: 
ky 
A grapes —; aD a ee a gran a ernlen. de 1 E a (70). 
Durch weitere Spezialisierung folgt für die: 


4 7 sck 
Dreieckschwingung n= 2 leo r SE Ysleiktissad: Arab Zi 


Rechteckschwingung Yy =) A rechteck — | 


Der Grund für das Versagen der Rechtecknäherung ist leicht einzusehen: die Geschwindig- 
keit x wäre in diesem Falle ja nur an den Sprungstellen y= nn von Null verschieden; da aber 
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nur das Vorzeichen von & und nicht der absolute Wert in die rechte Seite von (60) eingeht, hat 
das Glied sgn & überhaupt keinen Einfluß auf die Schwingung. r 

Eü In Bild 2 ist der Fehler der beiden Näherungslösungen (62) und (71) gegenüber der as 
Lösung (61) als Funktion der Dämpfung D aufgetragen. Wie zu erwarten ‚ist, stimmt Ei ar- 
monische Näherung gut im Bereich kleiner Dämpfungen, jedoch gibt die Dreiecksnäherung 
für D < 0,6 bessere Werte. 


Bild 2. Prozentualer Fehler der harmonischen Näherung (Kurve 1) und der Bild 3. Periodische Lösung von Gleichung (72) im Falle 7 >1 
Dreieck-Näherung (Kurve2) für die Amplituden des Schwingers nach 
Gleichung (60) 


b) Erzwungene Schwingungen eines aperiodischen Systems 


Es sei ein System betrachtet, dessen Verhalten durch die Differentialgleichung 
T<+xı=kyw(), (ve DL) a FREE TR 


beschrieben wird. Y(y) soll die Rechteckfunktion sein, die aus der Trapezfunktion von Bild 1 
mit y, = 0 folgt. Für alle Werte der Zeitkonstanten T ergeben sich periodische Lösungen, 
also stationäre Schwingungen, deren Amplitude leicht zu berechnen ist. Für0 <y <z hat (72) 
die Lösung: 


t 
ee 
Sie genügt den Bedingungen 
z(0)= —A, ke) 


Periodische Schwingungen treten auf, wenn 
x(z) = — x(0) 


ist. Nach Einsetzen dieser Periodizitätsbedingung in (73) bekommt man eine Bestimmungs- 
gleichung für die Amplitude A mit der Lösung: 


ER 
2T 


Schon durch eine qualitative Betrach- 
tung siehtman, daß eine harmonische Näherung 
y zumindest in den Grenzfällen sehr großer und 

sehr kleiner Zeitkonstante nicht besonders genau 

sein kann. In beiden Fällen weicht die wirk- 

liche Schwingung stark von einer Sinus- 

schwingung ab. Für T>1 ergeben sich prak- 

tisch Dreieckschwingungen (Bild 3), während im 

anderen Grenzfall T’< 1 näherungsweise Recht- 
eckschwingungen resultieren (Bild 4). Folglich werden wir im ersten Fall mit einem Dreiecks- 
ansatz, im zweiten Fall mit einem Rechteckansatz die besten Ergebnisse erwarten. Beide kön- 
nen durch Spezialisierung aus dem allgemeinen Trapezansatz gewonnen werden. 


Zunächst sei die harmonische Näherung betrachtet. Die auf der rechten Seite von (72) 
stehende Rechteck-Erregerfunktion kann in der Form 


A=ktgh 2 TA 


Bild 4. Periodische Lösung von Gleichung (72) im Falle 7 al 


Yaly) =.Sgn (Sin), We . (75) 
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geschrieben werden. Mit dem harmonischen Ansatz 


Ey HD 3 ER le ee ner (76) 
können wir (75) umformen in: 


a 
yr(y) = sgu De a (77). 


Setzt man dies in (72) ein, dann folgt eine nichtlineare Differentialgleichung, deren eines 
Glied ein um den Betrag 7 voreilendes Argument hat. Wie in einer früheren Arbeit [11] gezeigt 
wurde, lassen sich auch derartige ‚Totzeitprobleme‘“ nach der Methode der Harmonischen 
Balance erfassen. Man hat das nichtlineare Glied durch einen linearen Ausdruck zu ersetzen: 
ıy+9) 

A 


k sgn —hsen [EG er ar ara ee (78) 


und findet die Koeffizienten nach der in [11] angegebenen Vorschrift aus: 


2 
157 s : 4kcost 
mr, | k sgn (sin y) sin ydy = > 
(Ü 


a | ae 2 (79) 

h 1 1 y 4 ksin 

ee 3 u 

a sin ty | k sen (sin y) sin y dy Ey 
J ) 
Damit kann nun das lineare Ersatzsystem in die Form gebracht werden: 
4ksint\. 4kcost 

(ae) + A )» Via. ic Muss (80). 


Wenn periodische Lösungen möglich sein sollen, dann müssen notwendigerweise die beiden 
Koeffizienten dieser Gleichung verschwinden. Das ergibt zwei Bedingungen für die Größen A 
und r. Sie führen zu der Lösung: 


A 4 k 
harmonisch En: 
ar a En ie ne (81). 
COST — 7 
rn 


Die Dreieck- und Rechteck-Näherung sollen durch Spezialisierung der allgemeineren 
Trapeznäherung gewonnen werden. Wir setzen an: 


N ck A ee (82), 
und können jetzt für die Erregerfunktion von (72) schreiben: 


ı(y-+r) 
A 


Yrly) = sen [yrly)] = sen = sgen[e (y+ 9]. 


Die Koeffizienten für den zu bildenden linearen Ersatzausdruck (78) werden nun nicht aus 
(79), sondern wegen (55) aus entsprechend modifizierten Formeln gewonnen: 
2n 


1 kcost, 
a nat | Kosn Irlvntnd = a. 
Menge een (83). 
e 1 ksint 
dt.=-sin FR k sgn [yr(y)] Yrly) dy = nA. 2 (Ra — yo) 
0) 
Wegen (57) kommt: | 
es nie Sl ee 


Bm—4y)A ’ Any A 


Damit hat man das lineare Ersatzsystem: 
rs 32 ES IR ME VERRESENERTESEST SETZE (85). 
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Entsprechend den früheren Überlegungen führt die Forderung nach ‚dem ae 
der Koeffizienten dieses Systems zu zwei Gleichungen für A und z, die die Lösung haben: 


A PT: 3ka—yı) | 
ee. 
1 
cosr= yı+ Tr 
Durch Spezialisierung folgt daraus: 2 
Arne ArT RES (87), 
> Bert: 


A preieck = Mr Ar Te 
Durch Vergleich mit der harmonischen Näherung (81) findet man, daß sich die drei ver- 
schiedenen Näherungen wie folgt verhalten: 


2, E er 
A Rechteck : een, 5 A preieck =1: 1,27 3 1,5 = 


Im Grenzfall T — 0 findet man aus (74) für die strenge Lösung A— k. Derselbe Wert 
folgt aus (87) für die Rechtecklösung, während die anderen beiden Näherungen (81) bzw. (88) zu 
große Werte ergeben. Wie zu erwarten war, stimmt also hier die Rechtecknäherung am besten, 
sie gibt im Grenzfall T = 0 sogar genaue Werte. 

Im anderen Grenzfall T— oo folgt durch Entwicklung von (74): 

na k 


A>—- — 


21 


Mit diesem Wert stimmt — wie erwartet — die Dreiecknäherung am besten überein. 


Bild 5 zeigt die Fehler der drei verschiedenen Näherungslösungen (81), (87) und (88) als’ 


Funktion von T in halblogarithmischem Maßstab aufgetragen. Man erkennt daraus, daß man 


Bild 5. Prozentualer Fehler der Rechteck-Näherung (Kurve 1), der harmonischens Näherung (Kurve 2) 
und der Dreieck-Näherung (Kurve 3) für die Amplituden des Schwingers nach Gleichung (72). 


für T < 0,6 am besten mit der Rechtecknäherung und für T > 1,5 am besten mit der Dreieck- 
näherung rechnet. Nur in einem recht kleinen Zwischenbereich um den Wert T = 1 herum gibt 
die harmonische Näherung bessere Werte. 

Natürlich muß bei der Anwendung des hier beschriebenen Verfahrens die Form der Schwin- 
gung, die man berechnen will, schon näherungsweise bekannt sein, damit man die geeignete Appro- 
ximationsfunktion wählen kann. Vielfach lassen sich durch qualitative Überlegungen — oder 
auch durch Versuche — schon gewisse Anhaltspunkte über den Schwingungstyp gewinnen. Das 
gilt zum Beispiel für die technisch viel verwendeten Kippschwingungen und gewisse Typen von 
Relaxationsschwingungen. Hier wird man das vorgeschlagene Verfahren mit Vorteil verwenden 
können. 
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KLEINE MITTEILUNGEN 


Use of French curves for determination of 
tangents and center of curvature. 


If the proper use of French curves is considered to be 
suffieient to interpolate a set of distinct points p by 
a smooth curve /, then the following shows that under 
this assumption the construction of tangents and 
centers of curvature is possible also in a given set of 
points p; wherei = 1,2,3,4,5... Let p, be the point 
in which the tangent and the center of curvature should 
be determined. 

The construction of the tangent can be done in the 
following way: Connect by straight lines point 3 with 
points 1, 2, 4, 5,... Intersect these lines by a steady 


Tangente 


curve l,, say a straight line; then this steady curve l,, 
is an image of the curve l and the intersections 1’, 2’, 
4’,5’,... are the images of the points 1,2,4,5... The 
image 37 of 3 would be the intersection with the sougth 
tangent in point 3. Hence one has to determine a se- 
cond image 1, of the curve. Let the images of points 
Na48.,0..Be1l527, 4'',5’’,. .„. Let these images p” 
have distances d;= pf —p, from the images pj- 
Ltd = (H—P) (pi p,), where f (pi — p;) iS 
steady and positive. E. g. f(pi — p,) = vonstant = Y. 
Then by use of the French curves the points 1, 2’, 


4,5”, ... may be interpolated by the second image In 
of the curve I. The intersection of the pieture l, and I, 
is the image 3’ = 3°’ and therefore lies on the tangent. 
The construction of the center of curvature may be 
done in the following way: 

Construct in the points i = 1, 2,4, 5 the normals n; 
to the secants s; and the normal n, in point i = 3 to 
the found tangent. The normals n,, N, N4, 05. - » inter- 
sect the normal n, and together with the secants 8,, 83; 
Sy 85... form rectangular triangles with the normal 
n,. Determine the mid points 1’, 2°, A BEERDE 
the normals n}, 5, 04, 2,. Interpolate the curve l,through 
these points by use of French curves. The intersec- 
tion 3” with the normal n, to the found tangent in 
point 3 then is the center of curvature of the eurve I 
in the point 3. 


Cincinnati 15, Ohio F. S. Weinig. 


Der Konchoidenlenker für sechs unendlich 
benachbarte, Lagen 

1. K. Hoecken!) hat aus der konchoidischen Be- 
wegung, d.h. getrieblich aus einer zentrischen Winkel- 
schleife (Schubschleife), Bild 1, eine zentrische Kurbel- 


Bahn vonA 
Bild 1 


Bild 2 


schleife entwickelt, bei welcher eine Koppelkurve mit 
einer Geraden sechs endlich benachbarte Punkte ge- 
meinsam hat. Man kann aber auch auf gleichem Wege 
zu einer Koppelkurve kommen, welche mit einer Ge- 
raden sechs unendlich benachbarte Punkte gemein- 
sam hat, d.h. bei welcher die Gerade Tangente an die 
Kurve ist und eine sechspunktige genaue Berührung 
mit sehr guter Geradführung vorliegt. 

3, Bei dem Ausgangsgetriebe gleitet die Gerade 
(Stange) g durch den festen Punkt @ (durch eine um @ 
drehbare Hülse bzw. Kulisse) und der Punkt K wird 
auf einer festen Geraden x geführt. Alle mit der be- 


1) Hoecken, K.: Rechnerische Ermittlung eines Konchoiden- 
lenkers. Arch. f. Getriebetechn. 6 (1938), 8. 421. 
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wegten Ebene verbundenen Punkte beschreiben Kon- 
choiden, und für die Bahn des auf g selbst gelegenen 
Punktes A gilt mit KA=u und @ als Polarzentrum die 
Polargleichung?) 


Gereuor 
c0o8 


U)» 


Denkt man sich nun um den Punkt M auf der Sym- 
metrieachse einen Kreis vom Radius a gezogen, Bild 2, 
so trifft dieser Kreis mit der Polargleichung 


r=dcosp + Ja? +d2sin?p, 


Er ER ER ELDN, 
= dcosp + Ya? — d? + d? cos? 


worin M@=d gesetzt ist, die Konchoide in sechs 
Punkten. Denn Gleichsetzen von (1) und (2) liefert 
eine in cog@ = 2 kubische Gleichung von der Form 


2+4A2+B2+0=0..... (3), 


worin 
1a-_W+t2rd4d—a) „pP o__® 
2 ud x de 2ud 
(4)» 


gesetzt ist. Da sich drei Werte für 2= cos @ ergeben 
hat man also sechs Winkel , von denen je zwei sym- 
metrisch zur y-Achse liegen, Bild 1. 


3. Nach dem Wurzelsatz von Vieta gilt mit den 
Lösungen z2,, 2,, 2, der Gl. (3) doch 

At2,+23=—A (a), 

ar 2a ae Bible). 

Hha=C (od 


Bild 3 


Wenn nun die drei Lösungen für 2 zusammenfallen 
sollen, so muß aus Symmetriegründen A =a=%—], 
d.h. 9=0 sein, und es muß gelten 


(u? +2 pd + d2 —.a?) 
ET 


?) Meyer.zur Capellen, W.: Das Konchoidenpendel. 


2. In- 
strkde. 52 (1932), 123—128. 


Daraus, d.h. aus [6b], [6e], [6a] folgen dann die Maß- 
beziehungen 


d= oder p=3d, 
Hr 94 a 
u=, odeu=-; R Nez). 
aa el 
Aber oder a= h 


Der betrachtete Kreis ist jetzt der Krümmungskreis 
in der Symmetriestellung und hat dort mit der Kon- 
choide sechs unendlich benachbarte Punkte gemeinsam, 
Bild 3, sofern die in-Gl. (7) gegebenen Maßbedingungen 
erfüllt sind. Zur Kontrolle diene 


w=p+td-+a. 

4. Führt man jetzt den Punkt A auf seinem Krüm- 
mungskreis, d.h. unter Lösung der Fü bei K, so 
erhält man eine zentrische Kurbelschleife®), Bild 4, bei 
welcher a—=4A,A die Kurbel darstellt (M ist jetzt 
gleich A, gesetzt) und wobei die Gerade g durch den 
festen Punkt @, jetzt B,, gleitet. Dann stimmt dieses 
Getriebe mit dem Ausgangsgetriebe in sechs unendlich 


Bahn vonE& 


Bild 4 


benachbarten Lagen überein, d.h. der Koppelpunkt K 
beschreibt eine Koppelkurve mit sechspunktig berüh- 
render Tangente, Bild 4a. In der Symmetriestellung, 
d.h. für 9=0, ist K der Wendepol W und gleichzeitig 
der Ballsche Punkt der Systemlage, B, der Momentan- 


pol, also B,K =» gleich dem Wendekreisdurchmesser D 
und A=4° der Mittelpunkt des Rückkehrkreises, 
Bild 4b, wie sich auch aus allgemeineren Zusammen- 
hängen als Sonderfall ergab®). Die Euler-Savary- 
sche Formel 


1 1 —1 1 
—— = oder —— + — = — 
—PA —PA .2 a4. 
ist mit den oben gegebenen Maßen erfüllt. Diese Maße 


der Kurbelschleife sind also — um es zusammenfassend 
nochmals zu sagen, 


d —— 9d 
sn AK=u= —=9a. 
5. Würde man in K eine Schleife, wie in Bild 5a 
gezeigt, oder ein Schubglied für C, — 00 (oder auch eine 


1 


®) Vgl. auch: Meyer zur Capellen, W.: Die Harmonische Ana- 
Iyse bei der Kurbelschleife. Z. angew. Math. u. Mech. 36 (1956), 
Heft 3/4, ferner ders.: Kinematik und Dynamik der Kurbelschleife, 
Erste Mitteilung. Werkstatt und Betrieb 81 (1956), 8. 581585. 

*) Meyer zur Capellen, W.: Über Koppelkurven mit fünf- oder 
sechspunktig berührender Tangente in Sonderlagen des Gelenkvier- 
ecks. Konstruktion, im Druck; ferner ders.: Forschungsbericht 


herausgegeben vom Wirtschafts- und Verkehtrsminister des Landes 
Nordrhein-Westfalen, im Druck. 
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umlaufende Schleife) anbringen, so erhält man im 
Winkel der Schleife (bzw. im Weg des Schubgliedes) 
eine gute Rast, vgl. Bild 5b mit dem Verlauf des Aus- 
schlagwinkels y der Schleife und der zugehörigen 
Winkelgeschwindigkeit dy/di =y. 


Bild 5a 


Bild5b 


‚6. Der Flächeninhalt einer Koppelkurve, die von 
einem Punkt der Geraden g bei einer zentrischen 
Kurbelschleife beschrieben wird, errechnet sich zu?) 

"= — (a? n + uJ,) (8), 


worin w=AK und J, den folgenden Ausdruck be- 
deutet: 


Da Cr Hr 


r 


J„=2d ar 2) u + R(a, li 


en 2 1-4 ; 
Hierin ist & aus cos &= ——— zu bestimmen und 


I+4 
bedeutet A das Verhältnis = 7 ‚d.h. hier = =: und 
Be)‘ : 
sa; ferner ist F das elliptische Integral erster 


und E das elliptische Integral zweiter Gattung. Die 
Rechnung liefert dann F=—.a?rn 3,71, d.h. der 
Flächeninhalt der betrachteten Koppelkurve ist das 
3,7fache vom Inhalt des Kurbelkreises. 


Aachen. W. Meyer zur Capellen. 
5) Meyer zur Capellen, W.: Über den Flächeninhalt von 


Koppelkurven, Forschungsbericht herausgegeben vom Wirtschafts- 
und Verkehrsminister des Landes Nordrhein-Westfalen, im Druck. 


Ein Verfahren zur Berechnung der Stabilitäts- 
bedingungen bei Voraussetzung positiver 
Koeffizienten der charakteristischen Glei- 


» chung 


L. Cremer [1] hat die Frage der Verringerung der 
Zahl der Stabilitätsbedingungen bei Voraussetzung 
positiver Koeffizienten der charakteristischen Glei- 
schung untersucht. Er erhielt das Ergebnis, daß bei 
positiven Koeffizienten der charakteristischen Glei- 
chung die Hurwitzdeterminanten der Nebenfolge 


DH He 
stets positiv sind, wenn dies für die Hauptfolge 
Ha H„-3 Ve, ...)9 EI [#;] 


gilt, und umgekehrt. Bei ungeradem Grad n der cha- 
rakteristischen Gleichung sind stets die in eckigen 
Klammern stehenden Endglieder zu wählen. 

L. Cremer [1] hat dann bemerkt, „daß der Weg 
zu der am mühsamsten zu berechnenden (n—1)-stelligen 
Determinante H,„_ı ganz von selbst an allen H-Deter- 
minanten von niedrigerer Stellenzahl vorbeiführt, mag 
man unmittelbar an ihre Ausrechnung herangehen oder 
den praktisch vorzuziehenden Weg über den Routh- 
schen Algorithmus einschlagen“. 

Im folgenden wird aus einer Arbeit von Ch. Her- 
mite [2] ein einfaches „Diagonalschema‘““ hergeleitet, 
welches die Bildung der Hurwitzdeterminanten der 
Hauptfolge unabhängig von denen der Nebenfolge er- 
möglicht. Bei der Stabilitätsprüfung sind demnach nur 

Br n—1 
die Fr bzw. — GE 
folge zu bestimmen. W. Schmeidler [3] hat bereits 
ein „Randrahmenschema“ zur Bildung der Hermite- 
schen Stabilitätsbedingungen angegeben, das dem hier 
mitgeteilten Verfahren verwandt ist. 

Das Hermitesche Stabilitätskriterium [2] lautet 
wie folgt: Ist 


Hurwitzdeterminanten der Haupt- 


fe) =WH+ 24°‘ An? 


die betrachtete Gleichung mit reellen oder komplexen 
Koeffizienten, so bildet man den Ausdruck 


e * ln % n—l e 
‘ Y i, k=0 


Dabei ist f*(z) = f(—z) und fx) bedeutet das Polynom 
mit konjugiert komplexen Koeffizienten gegenüber 


Kennzeichnend dafür, daß die Wurzeln der Gleichung 
f(z)=0 sämtlich in der linken Halbebene liegen, ist, 
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daß die Hermitesche Form 
n—1 
> (—1)# A;,, U; Ur 
i, k=0 


positiv definit ist. 

Setzen wir nun reelle Koeffizienten voraus, SO Ver- 
schwindet Air, wenn i-+%k ungerade ist. Daher lassen 
sich die Determinanten 


wie folgt zerlegen (4): 


ee Mo 
D; m-ı = (d,m Cr m-ı 
(m == /1,02.0.00: 


er .0* 
Dam = (5m CT, m 


r n a 
(m=1, a 3 . 


ge 
Ch,m = 2", Hzm-ı 


TR I KOMTUR 
Ci,m=2 H,m- 


Hierbei gilt: 


Den Bedingungen 


DV 5.2 So 70) 


gleichwertig sind offenbar die Ungleichungen: 


Co,m = % Hz m-ı > 0 (m= 1, De 


C a >o0 Im=1,2,...., 4° 


| 
'ım 2 


Die Hermiteschen Stabilitätsbedingungen führen 
also bei reellen Koeffizienten zu einer Trennung der 
Hurwitzdeterminanten der Haupt- und Nebenfolge, 
wie sie beim angegebenen Satz von L. Cremer [1] 
auftritt. Hieraus ergibt sich ein einfaches ‚‚Diagonal- 
verfahren‘‘ zur unabhängigen Berechnung der beiden 
Determinantenfolgen: 


1. Man bildet zunächst die Dreiecksmatrix A mit 
den Eckgliedern 


a4,05 no Ind dng: 
41,0 
A= | 
| 
\@n, 0 ! F 5 or 
Dabei gilt: 


Wenn ö + k gerade ist: ,;. = 0. 
Wenn i -+-k ungerade ist: 


—(20% 
Ad; = ik 
ik m a; a7 


i gerade 
‘ ungerade. 


Bei der Gleichung 7. Grades erhält man zum 
Beispiel: 
4% 
07 —a,a, 
A; dg 0 a,@, 
7 
et 0 ° —aa 0 —aa; 
K R x 
a, Ay 0 a, 0 aa 
4 £ 
0 °— 4 0 —., 0 —ad 
K K 
A, Ay 0 aa, 0 ara, 020,0, 
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2. a) Man betrachtet die Zeilen 1, 3, 5, ..., an 
n [n 
in A. Wir erhalten die Zeile m (m= Ir Fr paul 
der Dreiecksmatrix B,, indem wir von den nicht ver- 
schwindenden Gliedern der (?m—1)-ten Zeile von A 
diagonal nach links unten gehen und die Summe aller 
Elemente (einschl. des Ausgangselementes) bilden. Von 
der Dreiecksmatrix B, geht man durch symmetrische 
Ergänzung zur vollständigen Matrix C, über. 
Für die Hauptabschnittsdeterminanten C, m vom 
m-ten Grad der Matrix (, gilt: ; 


n [n—+1 
Co, m = % Hr m-ı (m = 1, ee a 
Bei geradem n stellen die Determinanten C,, „ die 
Hauptfolge dar. 
Bei ungeradem n erhält man die Nebenfolge. 


Beispiel für n=T: j 
A454, Adg 

(030,—Q44, +4,54) (A54,—4gA, +4,05) Aryl, 

(4,0,—454,)+4,9,) (A,4;—QgAz;+470,) a0; 

0,0, Anüı AG] 

b) Man betrachtet die Zeilen 2, 4, 6, ...,n [n—1] 

in A. Wir erhalten die Zeile m (m= ER 3 Bi 


der Dreiecksmatrix B,, indem wir von den nicht ver- 
schwindenden Gliedern der 2 m-ten Zeile nach links 
unten gehen und die Summe aller Elemente (einschl. 
des Ausgangselementes) bilden. Dabei sind alle Sum- 
manden mit entgegengesetzten Vorzeichen zu sum- 
mieren. Von der Dreiecksmatrix B, geht man durch 
symmetrische Ergänzung zur vollständigen Matrix (, 
über. 
Für die Hauptabschnittsdeterminanten C,,„ vom 
m-ten Grad der Matrix (, gilt: 
n [n—1 ) 
BEER 
| @d,d, —4zQy 


| Ag, —Ardg 
C1= 01 — 0500 (440; —AzQ, +40) — 4,5) Agds— rl, | 
|‚% 4ı — 474g AsSTA7Q, 

Bei ungeradem n stellen die ‘C‘,, , die Determinanten 
der Hauptfolge dar. Bei geradem n erhält man die 
Nebenfolge. 

Das Ergebnis von L. Cremer [1] führt damit zu 
dem Rechenvorteil, daß bei geradem n nur die Ma- 
trix C,, bei ungeradem n nur die Matrix C, zu be- 
rechnen ist. An Stelle von je einer Determinante vom 
Grad 1 bis a — 1 bei den Hurwirtz bedingungen ist nur 
n—1 
2 
berechnen, wobei im letzteren Falle die Elemente der 
Determinanten leicht zu berechnende Koeffizienten- 
funktionen sind. Das mitgeteilte Schema läßt sich 
auch benutzen, um die Koeffizientenbedingungen für 
Impedanzfunktionen von elektrischen Netzwerken mit 
zwei Widerstandsarten zu bilden [5], [6]. 


Adzdg 


C,m = Ham m=1,..., 


Beispiel für n =": 


Ad, —A,4g 


%Al; 4704 


je eine Determinante vom Grad 1 bis > bzw. zu 
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 BUCHBESPRECHUNGEN 


C. Carath6odory (weil. Prof. a. d. Universität Mün- 
chen), Variationsreehnung und partielle Dif- 
ferentialgleichungen erster Ordnung. Bd.I. 
Theorie der partiellen Differentialgleichun- 
gen erster Ordnung. Zweite Auflage, herausge- 
geben von Dr. E. Hölder, o. Professor an der Uni- 
versität Leipzig. 171 S. Leipzig 1956. B. G. Teubner 
Verlagsgesellschaft. Preis geb. 14,— DM. 

Reichlich zwanzig Jahre liegen nunmehr seit dem 
Erscheinen der ersten Auflage des vorliegenden Wer- 
kes zurück, und seit dieser Zeit hat es durch die Neu- 
artigkeit seiner Behandlungsweise von Variations- 
problemen — die man heute direkt als Carath6o- 
dorysche Methode bezeichnet — und durch die Ele- 
ganz seiner Darstellung nicht nur die Aufmerksamkeit 
vieler Fachkollegen auf sich gelenkt, sondern darüber 
hinaus einer großen Zahl von Studierenden als unent- 
behrliches Lehrbuch gedient. Da dieses Buch bereits 
seit Jahren vergriffen war, kann man die Herausgabe 
der zweiten Auflage durch Prof. Dr. E. Hölder wärm- 
stens begrüßen und daneben auch die Tatsache, daß 
das ursprünglich einbändige Werk, welches aus den 
beiden Teilen ‚Theorie der partiellen Differential- 
gleichungen erster Ordnung“ und „Variationsrech- 
nung‘“ besteht, jetzt in zwei Bänden erscheint (denen 
sich später noch ein dritter Band als Ergänzungsband 
hinzugesellt); denn dadurch wird dem Leser, der 
unter Verzicht auf ein näheres Studium der Variations- 
rechnung sich lediglich über die Theorie der partiellen 
Differentialgleichungen erster Ordnung informieren 
will, weitgehend entgegengekommen. 

Die Darstellung von Band I stimmt im wesentlichen 
mit der des ersten Teiles der alten Auflage überein, 
getreu der Auffassung des Herausgebers, „daß die 
Diktion Carath6odorys, die seine ganze Persönlich- 
keit widerspiegelt, nicht durch Änderung des Textes 
verwischt werden dürfte“. Aus diesem Grunde sind 
auch Ergänzungen nur in Fußnoten bzw. durch 
kleine Bemerkungen im Anhang von Bd.I gegeben. 
Beachtenswert erscheinen hier besonders die Er- 
weiterungen zum Begriff der stetigen Konvergenz. 
Zu den Kapiteln der Berührungstransformation und 
des Pfaffschen Problems verweist der Herausgeber 
auf weitere Ergänzungen im Bd. III. 


Leipzig. R. Klötzler. 


A. Piluger (Prof. a. d. ETH Zürich), Theorie der 
Riemannschen Flächen. (Grundlehren der mathe- 
matischen Wissenschaften, Bd. LXXXIX. XII + 
2848. m. 5 Abb. Berlin/Göttingen/Heidelberg 1957. 
Springer-Verlag. Preis geh. 36,— DM, geb. 39,20 DM. 

Wenn trotz hervorragender Bücher über Riemann- 
sche Flächen nun ein weiteres Buch über dieses Thema 
erschienen ist, so ist dieser Umstand dadurch gerecht- 
fertigt, daß hier eingehend die diversen Methoden zur 
Führung von Existenzbeweisen untersucht und ver- 
glichen werden. Zuerst wird das Perronsche Ver- 
fahren besprochen nebst den zugehörigen sog. kon- 
struktiven Varianten (Verfahren von Schwarz, 
Poincars, Kellog). Im Gegensatz zu diesem Ver- 
fahren, das nieht für alle Riemannschen Flächen 
brauchbar ist, läßt sich das Dirichletsche Prinzip 
auf beliebige Riemannsche Flächen anwenden. Es 
sei erwähnt, daß sich das klassische alternierende Ver- 
fahren als konstruktive Variante des Dirichletschen 
Prinzips herausstellt. Die Riemannsche Fläche wird 
als ein zusammenhängender Hausdorffscher Raum 
mit einer zweidimensionalen konformen Struktur auf- 
gefaßt; es wird hierbei nicht gefordert, daß auf der 
Fläche eine abzählbare Umgebungsbasis existiert. 
Einige Einzelheiten: Das erste Kapitel behandelt 
außer dem Begriff der Riemannschen Fläche u.a. 
die Kompaktifizierung einer solchen Fläche und 
orientierbare topologische und differenzierbare Flä- 


chen. Die nächsten zwei Kapitel beschäftigen sich mit 
analytischer Fortsetzung, Überlagerungsflächen, Ho- 
mologie und Cohomologie. Das vierte Kapitel ist den 
oben genannten Methoden, Existenzbeweise zu 
führen, gewidmet. Das nächste Kapitel geht auf die 
Uniformisierungstheorie ein; hier wird auch der 
Riemannsche Abbildungssatz bewiesen und die 
Riemannsche Fläche als Fundamentalbereich einer 
Gruppe linearer Substitutionen betrachtet. Durch 
die Stichworte Abelsches Differential, harmonisches 
Differential, Typenproblem ist der Inhalt der letzten 
zwei Kapitel kurz gekennzeichnet. Ein ausführliches 
Literaturverzeichnis befindet sich am Ende des außer- 
ordentlich interessanten und sehr modern geschrie- 
benen Buches, das dem am Gegenstand interessierten 
Leser eine Fülle von Wissen vermittelt. 


Dresden. M. Landsberg. 


Rudolf Sonntag, o. Prof. a. d. Techn. Hochschule 
Karlsruhe), Aufgaben aus der Technischen 
Mechanik. Graphische Statik, Festigkeits- 
lehre, Dynamik. XI, + 2098. m. 324 Abb. 
Berlin / Göttingen / Heidelberg 1955. Springer-Verlag. 
Preis geh. 16,50 DM, geb. 19,50 DM. 

Die vorliegende Sammlung bringt Aufgaben zu den 
im Untertitel angegebenen Gebieten der Technischen 
Mechanik. Der größte Teil der Aufgaben entstammt 
der technischen Praxis und ist erfreulicherweise neu. 
Der Verfasser legt großen Wert auf eine präzise For- 
mulierung der Aufgabenstellung und auf eine gründ- 
liche Darstellung der Lösungen, was auch notwendig 
ist, da die meisten Aufgaben nicht ganz einfach sind. 
Das Fehlen einfacher Routinebeispiele wird dem An- 
fänger Schwierigkeiten bereiten, für einen Fortge- 
schritteneren ist die vorliegende Aufgabensammlung 
eine gute Schule für selbständiges Denken. 


Dresden/Darmstadt. A.Kromm. 


Dr.-Ing. Istvän Szabö (o. Prof. a. d. T. U. Berlin- 
Charlottenburg), Höhere Technische Mechanik. 
XII + 472 8. m. 402 Abb. Berlin/Göttingen/ 
Heidelberg 1956. Springer-Verlag. Preis 31,50 DM. 

Das vorliegende Buch ist eine Fortsetzung der ‚„Ein- 
führung in die Technische Mechanik“ des Verfassers 
und umfaßt die zweite Hälfte seiner Vorlesungen und 
des dazugehörigen Übungsstoffes. Den Inhalt des 
Buches ersieht man aus den Abschnittsüberschriften: 

I. Prinzipien der Mechanik, II. Ausgewählte Pro- 
bleme der höheren Elastizitätstheorie, III. Einblick in 
die Plastizitätstheorie, IV. Theorie der Flüssigkeiten 
und Gase. Da der Verfasser in jedem Abschnitt sehr 
viele Fragen und Probleme behandelt, so konnte der 
für die 472 Buchseiten ungeheuer umfangreiche Stoff 
nur sehr knapp dargestellt werden. So werden z. B. 
die Verträglichkeitsbedingungen der Elastizitäts- 
theorie kurz abgeleitet, auf die Bedeutung dieser Be- 
dingungen wird aber nicht eingegangen. Deswegen 
ist das Buch zwar als Nachschlagewerk durchaus 
wertvoll, zum selbständigen Lernen aber dürfte es zu 
schwer verständlich sein. 

Mit Rücksicht auf die Praxis war der Verfasser be- 
strebt, zur Entwicklung beizutragen, in der die theo- 
retische Mechanik und die technische Mechanik sich 
angleichen. Man muß aber sagen, daß dieses Vor- 
haben dem Verfasser nicht ganz gelungen ist. So 
werden z. B. von den, den Praktiker interessierenden 
Schwingungsproblem kaum einige behandelt. Trotz 
dieser Nachteile ist das vorliegende Buch sehr wert- 
voll, da es ein hohes wissenschaftliches Niveau hat 
und absolut zuverlässig ist. 


Dresden/Darmstadt. A.Kromm. 
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Dr. H. von Sanden (o. Prof. a. d. Techn. Hochsch. 
Hannover), Darstellende Geometrie. 4. bis 7. 
erw. Aufl. Mit einem Beitrag von Priv.-Doz. Dr.- 
Ing. K. Stohler, Hannover. 115 $. m. 125 Abb. im 
Anhang. Stuttgart 1956. B. G. Teubner Verlagsges. 
Preis 7,80 DM. : 

Die Neuauflage des bekannten Buches unterschei- 
det sich von früheren [vgl. ZAMM 11 (1931), 8. 336] 
durch einen neuaufgenommenen Abschnitt: Per- 
spektive mit geneigter Bildebene (nach Dr. Stohler). 
In diesem wird nach einem einführenden Beispiel die 
Theorie der Meßpunkte entwickelt und anschlie- 
ßend auf ein Beispiel mit verschieden gedeuteter 
Aufgabenstellung (Konstruktion des perspektivi- 
schen Bildes eines Bauwerks bzw. Rekonstruktion 
des Bauwerks aus einer photographischen Aufnahme) 
angewendet. 


Dresden. A.Schubert. 


Dr. Wolfgang Haack (o. Prof. a. d. Techn. Univ. 
Berlin), Darstellende Geometrie. Dritter Teil: 
Axonometrie und Perspektive (Sammlung 
Göschen, Bd. 144). 127 S. m. 100 Abb. Berlin 1957. 
Verlag Walter de Gruyter & Co. Preis 2,40 DM. 

Mit dem vorliegenden dritten Band wird. das 1954 
begonnene Werk [vgl. ZAMM 35 (1955), S. 75] ab- 
geschlossen. 

Im ersten Abschnitt wird ausführlich die senk- 
rechte Axonometrie behandelt, und zwar bis zu 
Durchdringungen und Schattenkonstruktionen, wäh- 
rend die schiefe Axonometrie mit dem abschließenden 
Beweis des Satzes von Pohlke etwas kurz weg- 
kommt. Darauf folgt ein Abschnitt: Grundzüge der 
ebenen Perspektive, der im wesentlichen die Theorie 
der Darstellung von Punkten, Geraden und Ebenen 
enthält. Zahlreiche Konstruktionshilfsmittel und 
praktische Hinweise bringt der dritte Abschnitt. In 
zwei weiteren Abschnitten werden noch Perspektive 
von Kreisen und einfachen krummen Flächen sowie 
Schattenkonstruktionen in der Perspektive behan- 
delt. Den Schluß bildet ein Literaturverzeichnis und 
ein Register für alle drei Bände. 


Dresden. A. Schubert. 


R. Kraus, Getriebelehre. Band Ill, Maßbe- 
stimmung. 336 S. m. 200 Abb. Berlin 1956. Verlag 
Technik. Preis 19,— DM. 

Wie der Titel sagt, steht hier die Aufgabe, Ge- 
triebe für bestimmte Zwecke zu entwerfen und ihre 
Maße zu berechnen oder zeichnerisch zu gewinnen, 
im Vordergrund. Unter Wiederholung einiger, in 
früheren Bänden bereits behandelter Verfahren, 
werden zunächst Methoden zur graphischen Ermitt- 
lung von Geschwindigkeiten und Beschleunigungen 
angegeben und werden dann die synthetischen Wege 
rechnerischer und zeichnerischer Art auf Schub- und 
Schraubengetriebe und auf die ebenen viergliedrigen 
Kurbel- und „Kurbel-Schubgetriebe‘‘ angewendet. 
Späterhin beschäftigt sich der Verfasser mit den ent- 
sprechenden sechsgliedrigen Getrieben, schiebt aber 
dazwischen ein Kapitel über den Entwurf von Kurven- 
trieben, während der Abschluß durch ein Kapitel 
über „Maßberechnung mit Hilfe der Mittelpunkt- 
kurven“ gebildet wird. 

Die eigenwillige und persönliche Darstellung des 
Verfassers wirkt sehr anregend, so daß sie jedem 
kinematisch interessierten Leser sehr empfohlen 
werden kann. 


Aachen. W. Meyer zur Capellen. 


Dr.-Ing. habil. Alf Pflüger (Prof. a. d. T.H. Han- 
nover),, Elementare Schalenstatik. VII n- 
112 8. m. 56 Abb. Berlin/Göttingen/Heidelberg 1957. 
Springer-Verlag. Preis geb. 19,50 DM. 

Das vorliegende Büchlein ist unter verändertem 
Titel die zweite, erweiterte Auflage des von der 


Wolfenbütteler Verlagsanstalt 1948 herausgebrachten 
sog. Notdruckes „Einführung in die Schalenstatik a 
Es setzt an Vorkenntnissen einige der in den An- 
fangsvorlesungen vermittelten Grundlagen der Tech- 
nischen Mechanik voraus und ist zum Studium der 
Grundtatsachen der Schalentheorie sehr gut geeig- 
net. 

Nach der Definition der Schalen wird die Membran- 
theorie der Drehschalen an Beispielen des Kuppel- 
und Behälterbaues erläutert und das Meißner- 
Geckelersche Näherungsverfahren zur Abschätzung 
der Randstörungen hergeleitet. Ausführlich wird an 
Hand des Membranspannungszustandes die Trag- 
wirkung der Zylinderschalen dargelegt. Die Pucher- 
sche Differentialgleichungfür den Membranspannungs- 
zustand in Flächen mit rechteckigen Grundrissen, ein 
Abschnitt über die Spannungsberechnung aus den 
bekannten Schnittgrößen und eine Formelzusammen- 
stellung der Längskräfte von 69 Lösungen des Mem- 
branspannungszustandes beschließen das ansprechende 
Buch. 


Freiberg/Sa. D. Rüdiger. 


Die Registres der Berliner Akademie der 
Wissenschaften 1746—1766. Dokumente für 
das Wirken Leonhard Eulersin Berlin. Heraus- 
gegeben und eingeleitet von Eduard Winter. IX + 
393 S. Berlin 1957. Akademie-Verlag. Preis 28,50 DM. 

Das vorliegende Werk wurde von der Deutschen 
Akademie der Wissenschaften zu Ehren des 250. Ge- 
burtstages Leonhard Eulers herausgegeben, der 
von 1741—1766 als Mitglied der Berliner Akademie 
in Berlin weilte. Leider sind die eigentlichen 
Sitzungsprotokolle für die Jahre 1746—1766 nicht 
mehr vorhanden; dagegen existieren noch die von 
dem ständigen Sekretär der Akademie, J. H.S. For- 
mey, geführten Protokollauszüge, die ‚‚Registres‘“. 
Sie bilden trotz ihrer Kürze ein wichtiges Zeugnis für 
das Wirken Eulers in Berlin und für die Rolle, die er 
in der Zeit der deutschen Aufklärung gespielt hat. 

Der Herausgeber, der in dankenswerter Weise die 
Registres der Öffentlichkeit zugängig machte und 
mit vielen Erläuterungen versah, stellte in einer 
fesselnd geschriebenen Einleitung eine Reihe wich- 
tiger biographischer Tatsachen über Euler zusam- 
men Besonders ausführlich werden dort naturgemäß 
die Berliner Zeit und die starken Einflüsse gewürdigt, 
die Euler auf die Entwicklung der Akademie aus- 
übte. Auf diese Weise befriedigt das Buch nicht nur 
den Fachmann durch die Erleichterung des Quellen- 
studiums, sondern auch jeden gebildeten Laien, der 
sich für die geistigen Auseinandersetzungen jener 
Zeit interessiert. 


Dresden. K. Maruhn. 


K. Brücker-Steinkuhl, Anwendung mathe- 
mathisch-statistischer Verfahren in der In- 
dustrie (Forschungsberichte des Wirtschafts- und 
Verkehrsministeriums Nordrhein-Westfalen, Nr. 288). 
103 8. m. 28 Abb. u. 14 Tabellen. Köln und Opla- 
den 1956. Westdeutscher Verlag. Preis 24,20 DM. 

Mit dem Ziel, in der Metallindustrie Wirtschaftlich- 
keit und Qualität durch Einführung geeigneter ra- 
tioneller Prüfverfahren zu verbessern, schildert Verf. 
in allgemein verständlicher Weise spezielle, bei der 
Herstellung von Bandstahl und Ketten durchge- 
führte statistische Untersuchungen. Dabei handelt 
es sich im Wesentlichen um Anwendung bekannter 
statistischer Verfahren wie: Normalitätsprüfung em- 
pirischer Verteilungen mittels x2-Test, Abnahme- 
prüfungen, Kontrollkarten für Mittelwert, Spann- 
weite, Median- und Extremwerte, Ausschußzahl, ein- 
und zweifache Varianzanalyse, Bartlett-Test, Dif- 
ferenzmethode (O.N. Anderson), wobei die unent- 
behrlichen Formeln ohne Begründung den ent- 
sprechenden Kapiteln jeweils anhangsweise 'beige- 
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ehandelt und ver 
e Varianzanalyse ein-, zwei- 

en, lateinischer und griechisch- 
I Juadrate empfohlen. Der Leserkreis, 
e Schrift Nutzen zu bringen bestimmt und ge- 
t ist, dürfte bei den an der Planung, Über- 
wachung und Durchführung der geschilderten Pro- 
 duktionsprozesse unmittelbar Beteiligten zu suchen 
sein. Es ist zweifellos zu begrüßen, wenn durch 
Schriften der vorliegenden Art auch in der deutschen 
längst bewährte statistische 
Methoden Eingang bzw. stärkere Verbreitung finden. 


Bad Nauheim. M. P. Geppert. 


H. Rutishauser (Prof. a. d. ETH. Zürich), 
Der Quotienten - Differenzen - Algorithmus. 
(Mitteilungen aus dem Institut für angewandte 
Mathematik an der Eidgenössischen Technischen 
Hochschule in Zürich, Nr. 7), 748. Basel/Stuttgart 
1957. Birkhäuser Verlag. Preis brosch. 8,50 SFr., 
850DM 

In dem vorliegenden Heft faßt Rutishauser seine 
drei in den Jahren 1954/55 in der ZAMP veröffent- 
lichten Arbeiten, in den er den sogenannten QD-Algo- 
rithmus entwickelt, zusammen und macht sie da- 
durch erfreulicherweise bequemer und allgemeiner 
zugänglich. Die beiden ersten Teile haben dabei 
einige Änderungen erfahren. Es darf hier daran er- 
innert werden, daß es sich bei dem QD-Algorithmus 
kurz gesagt um ein Verfahren handelt, die Pole 
einer in der Form 


fe) = I, zu 


gegebenen Funktion numerisch zu ermitteln, eine 
Aufgabe, auf die auch andere wichtige Aufgaben der 
praktischen Analysis zurückgeführt werden können. 
Sein enger Zusammenhang mit der Kettenbruch- 
theorie wird gleich im ersten seine theoretischen 
Grundlagen und seine Eigenschaften behandelnden 
Kapitel deutlich sichtbar. In diesem Kapitel werden 
außerdem die zugeordneten Polynome eingeführt 
und die Beziehungen zu dem Biorthogonalisierungs- 
verfahren (BO-Algorithmus) von Lanczos und zum 
Verfahren der konjugierten Gradienten (eg-Algo- 
rithmus) von Stiefel-Hestenes besprochen. Es 
enthält auch schon Hinweise auf Besonderheiten, die 
bei der praktischen Durchführung des QD-Algorith- 
mus auftreten können. 
Im zweiten, den Anwendungen des QD-Algorith- 
mus gewidmeten Kapitel behandelt der Verfasser 
die Umwandlung einer Potenzreihe in einen Ketten- 
_ bruch, die Summation schlecht konvergenter Reihen, 
die Interpolation durch Exponentialsummen und 
sehr ausführlich die Bestimmung der Nullstellen 
einer algebraischen Gleichung. Diese letzte Aufgabe 
führt er auf die eingangs genannte Grundaufgabe 
zurück. Um zu einer praktisch zweckmäßigen Lö- 
sungsmethode zu kommen, entwickelt er die ‚pro- 
gressive Form des QD-Algorithmus“. Auch Einzel- 
fragen, wie die Modifizierung des Verfahrens beim 
Auftreten komplexer Wurzelpaare oder von Wurzeln 
gleichen Betrages und die Verbesserung der Konver- 
genz des Verfahrens durch geeignete Maßnahmen 
werden besprochen. 

Das dritte Kapitel beschäftigt sich eingehend mit 
der Bestimmung der Eigenwerte und Eigenvektoren 
einer Matrix mit Hilfe des QD-Algorithmus, und 


sind, falls die s dieSchwarzsch 
Matrix sind. Auch hier empfiehl 
führung des QD-Algorithmus in der prog 
Form. Ziel IR 
In einem ‘Anhang wird noch kurz auf 
Transformation, ein kontinuierliches Analogon zum 
QD-Algorithmus und die QD-Relaxation einge- 
gangen. £ 
Der 
führten Beispielen und von übersichtlichen Rechen- 


schemata. U.a. ist auch ein Flußdiagramm für die 


progressive Form des QD-Algorithmus eingefügt. 
Die Arbeit ist eine sehr wertvolle Bereicherung 
des Schrifttums der numerischen Mathematik. 


Dresden. H. Heinrich 


Dr.-Ing. Bruno Eck, Ventilatoren. Entwurf und 
Betrieb der Schleuder- und Schraubengebläse. 3. ver- 
besserte und erweiterte Aufl. 4938. mit 559 Abb. 
Berlin/Göttingen/Heidelberg 1957. Springer-Verlag. 
Preis geb. 45,— DM. 

Die Bücher von Dr. Eck haben innerhalb der 
technischen Literatur ungewöhnlich: viele Auflagen. 
Der Einfluß des Autors durch diese Bücher auf die 
Anschauungen und die Arbeitsweise der Ingenieure 
in Betrieben für Strömungsmaschinen ist nach meinen 
Erfahrungen sehr groß. Dr. Eck hat selbst die Ent- 
wieklungsrichtung der Ventilatoren maßgebend be- 
einflußt. Die Technik verdankt ihm in den letzten 
Jahren Hochleistungsgebläse radialer Bauart mit 
90% Wirkungsgrad. Obwohl es — stark auf experi- 
mentellen Reihenuntersuchungen aufbauend — zu 
solchen Erfolgen kam, wird bei der Beschreibung 
jeder Vorgang strömungstechnisch analysiert und 
versucht, mit elementaren Rechnungen zu einer 
mathematischen Beschreibung zu kommen. Wenn 
man häufig über mangelnden Kontakt zwischen an- 
gewandtem Mathematiker und dem Ingenieur Kla- 
gen hört, wenn viele Bemühungen zu einer besseren 
Verbindung gemacht werden, dann ist es sehr erfreu- 
lich, daß ein erfolgreicher Praktiker seine Erfahrun- 
gen in solcher Weise niederlegt. Das Material würde 
neben dem an gleichen Gestaltungen interessierten 
Konstrukteur auch dem Theoretiker manche Anre- 
gung zur weiteren vertieften Bearbeitung bieten. 

Allerdings äußert sich der Autor selbst voller 
Pessimismus gegen Versuche, die Strömung im radi- 
alen Laufrad mit Hilfe verfeinerter mathematischer 
Methoden genauer zu erfassen, da der Vergleich mit 
der Messung (wegen starker Ablösungserscheinungen) 
sehr große Unterschiede zeigte. Jedoch ist zu be- 
denken, daß diese Messungen an Maschinen mit ver- 
hältnismäßig geringen hydraulischen Wirkungsgra- 
den gemacht wurden, deren Durchströmen auch 
nicht als reibungsfrei idealisiert? werden kann. 
Macht man aber solche Messungen an Maschinen 
hohen Wirkungsgrades, wie es Acosta (A.S.M. E. 
1954) tat, so zeigte sich bessere Übereinstimmung 
mit der Theorie von Busemann (ZAMM 1928). 
Leider bleiben im vorliegenden Buch beide Arbeiten 
unerwähnt. Natürlich müßte der Berechnung einer 
reibungsfreien Strömung die Grenzschichtrechnung 
folgen. 

Das vorliegende Werk gliedert sich in 2 wesentliche 
Teile: Radialgebläse und Axialgebläse. Die reibungs- 
freie Strömung wird zum Stromfaden idealisiert, 
daraus mit einigen Annahmen die induzierten Ge- 
schwindigkeiten auf den Schaufeln berechnet. Er- 
fahrungswerte zur Vorausbestimmung der Verluste 
sind angeführt. 
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Gestaltungsfragen werden in allen Einzelheiten be- 
antwortet, ein breiter Raum wird dem betrieblichen 
Verhalten und der Regulierung zugebilligt. Über die 
Grundlagen der Lärmbekämpfung wird geschrieben. 


Eingegangene Bücher — Nachrichten 


Z. angew. Math. Mech. 
Bd. 37 Nr. 11/12 Nov./Dez. 1957 


Abschnitte über Festigkeitsberechnung und die ex- 
perimentelle Erprobung schließen dieses lesenswerte 
Buch ab. 


Dresden. W. Albring. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr. J. E. Hofmann (Honorarprof. a. d. Univ. Tü- 
bingen);, Geschichte der Mathematik. Dritter 
Teil. Von der Auseinandersetzung um den Cal- 
culusbiszurFranzösischenRevolution. (Samm- 
lung Göschen, Bd. 882). 107 S. Berlin 1957. Walter de 
Gruyter & Co. Preis brosch. 2,40 DM. 


Modern Computing Methods. (National Physical 
Laboratory). Notes on Applied Science, No. 16. 128 S. 
m. Abb. London 1957. Her Majestys Stationery Office. 
Preis brosch. 10 s. 10d. 


Dr. K. Prachar (Dozent a. d. Univ. Wien), Prim- 
zahlverteilung. (Die Grundlehren der Mathemati- 
schen Wissenschaften in Einzeldarstellungen mit 
besonderer Berücksichtung der Anwendungsgebiete, 
Bd. XCI). X +415 S. Berlin/Göttingen/Heidelberg 
1957. Springer-Verlag. Preis geh. 55,— DM, geb. 
58,— DM. 


Die Technische Hochschule Carolo-Wilhel- 
mina zu Braunschweig. Berichte aus Forschung 
und Hochschulleben 1954 bis 1957. Herausgegeben im 
Auftrage von Rektor und Senat von Hartwig Peter- 
mann. 262 S. mit zahlr. Abb. Braunschweig 1957. 
Verlag E. Appelhans & Co. 


Dr. 6. Hessenberg (weil. o. Prof. a. d. TH. Berlin), 
Ebene und sphärische Trigonometrie. Fünfte 
Aufl. bearb. v. Dr. H. Kneser (o. Prof. a. d. Univ. 
Tübingen). (Sammlung Göschen, Band 99). 172 S. m. 
60 Fig. Berlin 1957. Walter de Gruyter & Co. Preis 
geh. 2,40 DM. 


Dr. H. Hasse (o. Prof. a. d. Univ. Hamburg), 
Höhere Algebra. I: Lineare Gleichungen.Vierte, 
durchgesehene Aufl. (Sammlung Göschen, Band 931). 
1528. Berlin 1957. Walter de Gruyter & Co. Preis 
geh. 2,40 DM. 


Dr.-Ing. E. Hübner, Technische Schwingungs- 
lehre in ihren Grundzügen. XI -+ 3228. m. 
208 Abb. Berlin/Göttingen/Heidelberg 1957. Sprin- 
ger-Verlag. Preis geb. 29,40 DM. 


Dr.-Ing. W. Flügge (Professor of Engineering Me- 
chanies Stanford University), Statik und Dyna- 
mik der Schalen. Zweite, neubearb. Aufl. VIII + 
286 S.m. 121 Abb. Berlin/Göttingen/Heidelberg 1957. 
Springer-Verlag. Preis geb. 28,50 DM. 


6. Springer, Introduction to Riemann Sur- 
faces. VI + 3078. m. Abb. Reading, Mass. 1957. 
en Publishing Company, Inc. Preis 
9,50 $. 


H. v. Mangoldt, Einführung in die Höhere 
Mathematik für Studierende und zum Selbststu- 
dium. Seit der Sechsten Auflage neu herausgegeben 
und erweitert von Prof. Dr. Knopp (Tübingen). 
Zweiter Band: Differentialreehnung, unend- 
liche Reihen, Elemente der Differential- 


geometrieundderFunktionentheorie. 10.vollst. 
u. neubearb. Aufl. XIV + 6248. m. 115 Abb. Leip- 
zig 1957. S. Hirzel Verlag. Preis geb. 22,—.DM. 


R. Bellman, Dynamic Programming. XXV 4 
342 8. m. Abb. Priceton, New Jersey 1957. Prince- 
ton University Press. Preis geb. 6,75 $. 


Dr.-Ing. H. Wittke, Ultragon. Fünfstellige Tafeln 
für Sinus, Cosinüs, Tangens, Cotangens. Hamburg 
1957. Hanseatische Verlagsanstalt GmbH. 


Dr. Fr. A. Willers, Methoden der praktischen 
Analysis. Dritte Aufl. 429S. m. 93 Abb. Berlin 
1957. Walter de Gruyter & Co. Preis geb. 283,— DM. 


Jenaer Jahrbuch 1956. Herausgegeben von 
VEB Carl Zeiß Jena. 263 S. m. 139 Abb. u. 15 Tafeln. 
Jena 1957. Kommissionsverlag VEB Gustav Fischer. 
Preis geb. 20, — DM. 


H. S. M. Coxeter und W. O. J. Moser, Generators 
and Relations for diserete Groups. (Ergebnisse 
der Mathematik und ihrer Grenzgebiete, Neue Folge, 
Heft 14). VIII + 155 8. m. 84 Abb. Berlin/Göttingen/ 
Heidelberg 1957. Springer-Verlag. Preis geh. 32,— DM. 


C. F. Gauss Gedenkband anläßlich des 100. 
Todestagesam 23. Februar 1955. Herausgegeben 
von H. Reichardt, Berlin. VIII + 252 S. m. Abb. 
im Text. Leipzig 1957. B. G. Teubner Verlagsgesell- 
schaft. Preis geb. 25,80 DM. 


H. Ch. Sharaf, Beeinflussung der Lage des 
Verdichtungsstoßes bei Profilen in hoher Ge- 
schwindigkeit. (Mitteilungen aus dem Max- 
Planck-Institut für Strömungsforschung, Nr. 16). 
63 S. m. 42 Abb. Göttingen 1957. Selbstverlag des 
Max-Planck-Instituts für Strömungsforschung. Preis 
geh. 7,— DM. 


H. Reichardt, Über Vielfachverzweigungen 
für bestimmte Verteilungen der Durchfluß- 
menge. (Mitteilungen aus dem Max-Planck-Institut 
für Strömungsforschung und der Aerodynamischen 
Versuchsanstalt, Nr. 17). 498. m. 15 Abb. Göttingen 
1957. Selbstverlag des Max-Planck-Instituts für 
Strömungsforschung und Aerodynamische Versuchs- 
anstalt. Preis geh. 5,— DM. 
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Braunschweig. Professor Dr. H. Schlichting ist 
zum Wissenschaftlichen Mitglied der Max-Planck- 
Gesellschaft zur Förderung der Wissenschaften er- 
nannt worden. 


München. Der bisherige Privatdozent an der TH. 
Stuttgart, Dr.-Ing.E. Truckenbrodt, wurdeam 15.7. 
1957 zum ordentlichen Professor für Technische 
Mechanik in der Fakultät für Maschinenwesen und 
Elektrotechnik ernannt und zum Direktor des Insti- 
tuts für Strömungsmechanik bestellt. 
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A.1I. LURJE 


Einige nichtlineare Probleme 


aus der Theorie der selbsttätigen Regelung 


Übersetzung”aus dem Russischen 


In deutscher Sprache herausgegeben von Prof. Dr. H. KINDLER und Dr. R. REISSIG 


1957. XII, 164 Seiten — 24 Abbildungen — gr. 8° — Ganzleinen DM 15,— 


Das Werk enthält einige Arbeiten des Verfassers über nichtlineare Probleme der Regelungstheorie. Nach 
der Fixierung des Problems entwickelt er eine Methode zur Reduktion der Bewegungsgleichungen eines 
Regelungssystems auf die von ihm als kanonisch bezeichnete Form. Ferner behandelt er die Stabilitäts- 
kriterien für nichtlineare Regelkreise und geht im Anschluß hieran auf Selbstschwingungen in Regelkreisen 


und das Verhalten eines Systems an der Stabilitätsgrenze ein. 


In der deutschsprachigen Literatur sind die von Lurje erläuterten Probleme bisher nur wenig untersucht 
worden. Da sie sehr aktuell sind, kommt diesem Werk eine besondere Bedeutung zu. Mit seiner guten und 
verständlichen Darstellungsweise wendet sich das Buch in erster Linie an Ingenieure, die sich mit der 
Berechnung und dem Bau von selbsttätigen Regelungssystemen befassen, aber auch an Physiker und 


angewandte Mathematiker. 


Bestellungen durch eine Buchhandlung erbeten 
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